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K V Jayakumar, Ph.D. 

Department of Civil Engineering 
Indian Institute of Technology, Kanpur 
May, 1992 

FLOOD ESTIMATION IN SMALL CATCHMENTS - A REGIONAL CASE STUDY 
Summa ry : 

Estimates of flood are required for the design and economic appraisal of a 
variety of engineering works such as culverts, spillways etc. When adequate and 
appropriate data are available for the basin, frequency analysis may be applied or 
else unit hydrograph (UH) and other rainfall-runoff procedures may be applied to 
model the basin system. When this is combined with design rainfall, the design flood 
may be derived. Data available for small watersheds are generally inadequate for 
these approaches. Regional flood frequency analysis is generally not valid for small 
watersheds because of data limitations. Then, other empirical approaches for 
regionalisation are to be used to model either the variation of flood peaks or 
rainfall-runoff relationships. 

From available data, regression relationships are derived to relate peak 
discharge to hydrometeorological features of the storm and physiographic features 
of the watershed. Even when adequate data are available for establishing such 
relationships, they are limited to the regions for which such relationships have been 
derived. Since generally adequate data are available for estimating design storms, it 
seems worthwhile to derive relationships between storm rainfall and streamflow. The 
parameters of the model may then be related to the characteristics of the storm and 
the physiographic features of the watershed to yield a regional relationship between 


storm and floods. 
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A more realistic and alternate approach is to model the catchment runoff 
process according to physical principles. Information regarding physiography, soils, 
vegetation and the drainage pattern of the catchment is used in the physically 
based models. The catchment is decomposed into a network of channels and overland 
flow segments. Using the physical characteristics of each component and appropriate 
hydraulic equations, the hydrologic processes are simulated using the distributed 
parameter model. Necessary values of the coefficients and constants are specified 
from physical characteristics of the watershed or obtained from calibration of the 
model with observed data. 

This study is essentially a regional case study for flood estimation in small 
catchments in a part of the North Brahmaputra basin, designated as subzone 2(a) by 
the Central Water Commission (CWC), to compare in a very limited data environment 
the two approaches for rainfall-runoff relationships, viz., regionalisation using UH; 
and a physically based model . 

Objectives of the Study: The objectives of the present study are: 

i) To develop a regional relationship for UH parameters from available limited 
data for a number of basins to evaluate the validity and applicability of UH 
based regionalisation procedure for flood estimation; 

ii) To apply a physically based model to small watersheds in the region; study 
the regional variation, if any, of the parameters of the model and also to judge 
their validity for flood estimation; and 

iii) To compare these two approaches for regional flood estimation for small 
catchments. 

Scope of the Study: The study is limited to small watersheds in the western 
part of the North Brahmaputra basin, for which data were availble. Furthermore, 
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the UH is represented by the Nash conceptual model and the physically based model 
considered in the study is the USGS Distributed Routing Rainfall Runoff Model 
(DR3M). 

Data Used: Catchment plan, crosssection at the bridge site and hourly rainfall 

is 

and stage data at railway bridges for about 10 catchments in the region were 
available for the study. The catchment areas ranged from about 7.00 Km 2 to 215.00 
Km 2 . The data were available for a period ranging from two to three years. After 
preliminary analysis for data errors, only data of seven basins were found to be 
suitable for analysis. 

Details of the Study: 

In the UH approach for regionalisation, it is necessary to represent the UH in 
terms of a limited number of parameters. Nash model with two parameters N and K 
was adopted in the study. 

Attempts to correlate coefficient of runoff, the infiltration index and base 
flow with storm and runoff parameters were unsuccessful. Variation of the 
infiltration index among the storms was very large. However, a value of 5.00 
mm/hour can be considered to be the general minimum value for the infiltration 
index for all the basins, particularly during intense storms generally met with in the 
design. This value corresponds to the silty loam soil or sandy clay loam soil of the 
region. 

The various steps involved in the regionalisation of the parameters of this 
model are: 

i) Estimate the UH parameters for a storm in a basin; 

ii) Repeat the step for a number of storms in a basin to explain the variation, if 
any, of the parameters in terms of hydrometeorological characteristics of the 
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storm; and 

iii) Repeat the above two steps for a number of basins in the region so that the 
residual variation of parameters can be correlated to the physiographic 
characteristics of the basin. 

The UH parameters N and K of the Nash model vary from storm to storm and 
from basin to basin. A smaller NK generally leads to a quicker flood and a high, 
peak. A quasilinear variation of N and K as a function of peak discharge (Q p< ,» t ) is 
assumed and for each of the basin, the regression equations relating N, K and NK 
with Q peak are derived. Regional regression equations relating N and K corresponding 
to the maximum observed peak in each basin is as follows: 

N - 1.792 + 0.0088 A, and 
K - 2.766 + 0.016A 

where A is in Km 2 , K is in hours and N is dimensionless. 

Baseflow varies from storm to storm and generally increases with the basin 
area. Since the higher values of baseflow lead to larger floods, the enveloping 
curve for baseflow/Km 2 as a function of basin area is recommended for use with 
the design flood. 

A design flood estimated using the regional relationships for a small 
watershed in the region compares with other approaches generally used in India 
satisfactorily. 

A drainage basin is represented in the USGS DR3M as a set of channel and 
overland flow segments in such a way that the essential basin geometry and 
physiography are taken into account for runoff computation. The model has four 
components - a soil moisture or water balance component, the infiltration or 
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rainfall excess component, a routing component and an optimisation component. The 
model has seven physically based parameters - four soil moisture accounting 
parameters viz., (i) DRN, a constant drainage rate for redistribution of soil moisture 
in inches/day, (ii) EVC, a pan coefficient for converting measured pan evaporation 
to potential evapotranspiration; (iii) RR > the proportion of daily rainfall that 
infiltrates into the soil for period of simulation excluding unit rainfall days and 
(iv) Soil-moisture storage at field capacity in inches, and three infiltration 
parameters viz., (i) PSP, the suction at wetted front for soil moisture at field 
capacity in inches of pressure, (ii) KSAT, the effective saturated value of the 
hydraulic conductivity in inches/hour and (iii) RGF, the ratio of suction at wetted 
front for soil moisture at wilting point to that at field capacity. These parameters 
are fitted to the model for a basin by Rosenbrock’s algorithm and the iterative 
optimisation among basins is done heuristically. 

It is realised from the study that the directly contributing impervious area 
percentage to be given as input data from the characteristics of the watershed is 
very important in estimating the flood peak. Such data were not available for the 
region. They were parametrised and estimated by simulation and compared with 
general topography and basin area and are considered to be satisfactory. However, 
with the use of remote sensed data (not available for the study but generally 
available for government organisations), it will be possible to use this information 
for realistic modelling of small watersheds. The results of the study indicated that 
a constant parameter model for th® small catchements in the region is realistic. The 
model parameters viz., PSP — 10 inches; KSAT ■= 0.178 inches/hour; RGF — 18; and 
BMSN - 10 inches agree with the watershed characteristics. The design flood 
estimated for a small watershed using this model compared favourably with other 
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Cone 1 us i ons : 

1. Since design rainfall can be estimated separately by standard procedures, it 
seems preferable to use for estimation of design flood in a small basin, a regional 
rainfall-runoff relationship with design storm rather than using a purely 
empirical approach. 

2. For the design flood for the region under study, a constant inflitration 
parameter during the design storm, a baseflow related to the basin area and a 
regional quasilinear relationship for the UH parameters were found to be 
satisfactory. Regional UH approaches involve a large amount of empiricism. 
When sufficiently large amount of data are available they may be used for 
regionalisation. 

3. DR3M is a relatively simple distributed parameter model which can be used 
for simulation of runoff in small watersheds. Impervious area is a very 
important parameter in the DR3M which has to be estimated with care. 

4. Even in a very limited data environment, DR3M has been found to work 
satisfactorily. 

5. Since field data concerning topography, soil and vegetation may be easily 
collected by remote sensing and other means, DR3M may be used when data are 
very limited. 

6. Because of easy availability of, and competence with computers and softwares 
in field organisations, it seems preferable to adopt a distributed parameter model 
with relevant field data rather than empirical, statistical relationships of 
questionable validity. 

7. Distributed parameter models, may, eventually repalce empirical regression 
models in hydrology not only for large watersheds but also for small watersheds. 



CHAPTER 1 


INTRODUCTION 


1 . 1 Genera 1 

A flood is any large peak flow, but the term is generally used to mean a flow 
that exceeds the capacity of the banks of the stream at a given location, thus 
creating potential damage. When a stream is In flood at one location, it does not 
mean that it is In flood all along its length. The estimation of flood peak is location 
specific. It is very important for design of various water control structures. A 
very major problem in hydrology is the prediction of magnitude of flood peaks at a 
site for design purposes. 

The primary requirement for any such hydrologic analysis and design is the 
availability of data, especially those pertaining to rainfall and runoff in streams. 
Larger catchments, in general, are gauged at various locations. But for smaller 
catchments, the major constraint is the nonavailability of data in the stream at or 
near the site. A small catchment may be defined as one that is so small that its 
sensitiveness to high intensity rainfalls of short durations and to land use are not 
suppressed by the channel storage characteristics (Chow,1964). Numerous attempts 
have been made to delimit the small catchment either by area; or as function or 
type of storage. 

Only a limited number of small watersheds are gauged throughout the 
country and a few studies have been carried out on hydrometeorological aspects of 
small watersheds. Problem specific studies and data collection have been carried 
out by research organisations in India like International Crop Research Institute for 
the Semi- Arid Tropics (ICRISAT), Central Arid Zone Research Institute (CAZRI), 
Central Soil and Water Conservation Research and Training 
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Institute (CSWCRTI) and Research Designs and Standards Organisation (RDSO) of 
the Ministry of Railways. One of the primary reasons for the limited number of 
general studies has been the lack of data. In the absence of data, the engineer has 
to find ways and means to estimate flood peaks. Since small structures are 
associated with small watersheds, it is possible to accept some risk of failure of 
these small structures and so they are generally designed for a small probability of 
failure, say, in terms of recurrence interval of storms or floods for design. 
Estimation of floods for various recurrence intervals for watersheds where no 
observations has been made has been a challenging task for engineers who are 
involved in hydrologic analysis and design. 

Several approaches are available for estimation of floods. When adequate and 
appropriate data are available for the basin, frequency analysis may be used or else 
unit hydrograph (UH) and other rainfall-runoff procedures may be applied to model 
the basin system. When this is combined with design rainfall, the design flood may 
be derived. 

Data available for small watersheds are generally inadequate for these 
approaches. In such a case, empirical approaches are to be used to estimate the 
design flood. This can be through the development of a relationship either for the 
regional variation of flood peaks or for rainfall-runoff relationships. From available 
data, regression relationships are derived to relate peak discharges to 
hydrometeorological features of the storm and physiographic characteristics of the 
basin. Even when adequate data are available for establishing such relationships, 
they are limited to the regions for which such relationships have been derived. 


Regional frequency analysis is generally not possible for small watersheds 
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because of data leng“th limitation. Since generally adequate data are available for 
estimating design storms, it seems worthwhile to derive relationships between storm 
rainfall and stream flow, perhaps in terms of unit hydrographs. The parameters of 
the UH may then be re lated to the characteristics of the storm and physiographic 
features of the watershed to yield a regional relationship between storm and flood. 

A more realistic approach is to simulate the catchment runoff process 
according to physical pr inciples. Data concerning topography, soil and vegetation can 
be obtained from remotes sensing. The physically based distributed parameter models 
can now be implement «d even in a PC- environment. There has been also a greater 
availability of technicaE know-how on the applicability of such models. They may 
also eliminate the n«e«d£Tor regionalisation. Hence it seems reasonable to consider the 
use of physically Imsed distributed parameter models as an alternative to 
regionalisation fox flo-cd estimation in small catchments. Thus the study proposes to 
explore the possibility -of using physically based distributed parameter models for 
flood estimation and to compare them with a regionalisation procedure based on UH 
approaches. 

1.2 Statement of tBie Problem 

Small catchments *®re characterised by: 

i) relatively short da ration of storm, 

ii) a nearly uniform ra infall distribution in space and time, 

iii) steeper slopes end r elatively less permeable surfaces, 

iv) predominance of ov~erland flow and channel flow characteristics in comparison 
to storage characteristics, 

v) a nearly wet condition at the beginning of design storms, and 
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vi) generally a quick response of the catchment in producing runoff. 

Empirical relationships are often region specific and may not be available for 
the region under study; and small watershed data are generally of very short 
length thus precluding regional frequency approaches. Design storms for storm 
precipitation can be fairly well estimated from the data available with India 
Meteorology Department (IMD). Accordingly estimation of design flood in a small 
watershed reduces to establishing the rainfall-runoff relationships for the watershed 
and combining with design rainfall to produce design runoff including peak runoff. 

Rainfall-runoff relationships for small watersheds can be regionalised as 
follows: 

i) Estimate the UH parameters for a storm in a basin, 

ii) Repeat the step for a number of storms in a basin to explain the variation, 
if any, of the parameters in terms of hydrometeorological characteristics 
of the storm; and 

iii) Repeat the above two steps for a number of basins in the region so that 
the residual variation of parameters can be correlated to the physiographic 
characteristics of the basin. 

Alternatively a distributed parameter model may be fitted to the data from 
each watershed. Information regarding physiography, soils, vegetation and the 
drainage pattern of the catchment is used in the physically based models. The 
catchment is decomposed into a network of channel and overland flow segments. 
Using the physical characteristics of each component and appropriate hydraulic 
equations, the hydrologic processes are simulated using the distributed parameter 
model. Necessary values of the coefficients and constants are specified from 
physical characteristics of the watershed or obtained from calibration of the model 
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with observed data. In case the value of any of the parameters vary very 
significantly, then regional variation may be modelled appropriately. While such a 
variation may be large in large basin and in different regions, they are expected to 
be fairly small for small watersheds in a region. A comparison of the performance of 
the two approaches in modelling the rainfall-runoff relationships in small watersheds 
may perhaps lead to identification of better procedures. 

For purposes of hydrometeorological studies, India has been divided into 26 
subzones by the Central Water Commission (CWC), each of which is assumed to be 
hydrologically homogeneous. Research Designs and Standards Organisation (RDSO) of 
the Ministry of Railways has collected and in collaboration with CWC and IMD has 
analysed a large amount of hydrometeorological data for a large number of 
catchments draining into various railway culverts and bridges in different subzones. 
However they seem to have problems in the analysis of data and regionalisation in 
flood estimation for small watersheds because of large errors. Hence the estimation 
of floods in small watersheds, particularly when there are serious data limitations, 
is very important. Thus it seems worthwhile to consider the suitability of the 
models in a data scarce environment. 

1.3 Objectives of the Study 

This study is essentially a regional case study for flood estimation in small 
catchments in a part of the North Brahmaputra basin with a very limited data 
environment to compare the two approaches for rainfall-runoff relationships, viz., 
regionalisation using unit hydrograph (UH) approach and a physically based 
distributed parameter model. The objectives of the study are: 

i) To develop a regional relationship for UH parameters from available limited 
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data for a number of basins to evaluate the validity and applicability of 
UH based regionalisation procedures for flood estimation; 

ii) To apply a physically based model to small catchments in the region, study 
the regional variation, if any, of the parameters and also to judge their 
validity for flood estimation; and 

iii) To compare the two approaches for regional flood estimation for small 
catchments. 

1.4 Significance of the Study 

Development of regional procedures for estimation of design floods in small 
to medium sized drainage basins is an important problem in hydrology. Even though 
many methods have been developed and used widely, in many cases, there is a 
serious doubt about their adequacy and suitability. 

Small watershed hydrology has received limited attention in research and 
design because of a number of factors, viz., 

i) Small watersheds have small hydraulic structures and projects with low 
expenditure or financial investment; 

ii) A routine design carried out in a minimum time is considered to be 
satisfactory; 

iii) Unfortunately the large number of such structures or projects is not 
taken note of, and hence economy and scientific design in this area are 

ignored; and 

iv) This is not considered to be an intellectually or scientifically 
satisfying area of study for researchers (Pilgrim,! 989). 


Yet this area is becoming very important because of the potential saving in 
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the very large number of hydraulic structures and projects in the small watersheds. 
Any significant improvement for rational hydrologic design will lead to technically 
viable and economical system. Such procedures have been developed and applied in 
urban drainage systems which are relatively small. Accordingly, the study 
objectives are considered to be very important for rational scientific design with a 
better insight into the hydrology of small watersheds and hence the study is 

considered to be significant. 

1.5 Scope of the Study 

The scope of the study is limited because of a number of factors: 

i) Because of data availability relevant to the problem under consideration, 

the study is limited to a part of the North Brahmaputra basin; 

ii) While several models are available for establishing rainfall-runoff 
relationship using a black box approach, the study is limited to UH 
approach using Nash Model and a combination of optimisation and 
heuristics in parameter estimation for the model; 

iii) Similarly while several distributed models for a watershed are available, 

this study is limited to the USGS DR3M in case of physically based 

distributed parameter model; and 

iv) Some of the details of the analysis are limited by the inaccessibility and 
remoteness of the region for collection of additional data, and in the 
final stages due to nonavailability of remote sensed data in limited time. 

1.6 Organisation of the Thesis 


The study is reported in the following sequence: 
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i) Flood estimation procedures in small watersheds with particular reference 
to the study are briefly reviewed in Chapter 2; 

ii) Data availability, their limitations, preliminary analysis and selection of 
models used in the study are reported in Chapter 3; 

iii) Chapter 4 contains the unit hydrograph approach using the Nash Model 
used for regionalisation. The analysis of the various basins used in the 
study are also presented therein; 

iv) Chapter 5 contains a brief description of the USGS Distributed Routing 
Rainfall Runoff Model and its application to the basins chosen for the 
study; 

v) The last chapter (Chapter 6) presents the summary of the work carried out, 
the conclusions and suggestions for future work. 



CHAPTER 2 


FLOOD ESTIMATION PROCEDURES IN SMALL 
CATCHMENTS - A REVIEW 

2.1 Introduction 

Lack of adequate runoff information is one of the major problems faced world 
over by hydrologists and engineers who are involved in the design of water control 
structures and in the evaluation of surface water potential. Engineers are required 
to estimate the peak discharges at locations where there are adequate or dearth of 
data. The term dearth may mean nil or it may mean inadequate. When there are no 
data or when the available data are very limited, the hydrologist has to find ways 
and means of somehow rationally estimating the design flood. In such cases, it 
becomes necessary to interpolate or extrapolate the available data for the basin 
under study in time and over events, or use the information from nearby sites in 
the basin or from neighbouring streams and basins. These approaches are used to 
establish spatial correlations or regional relationships, as the case may be, for flood 
estimation. If the regional approach is used, considerable care should be taken to 
select streams as nearly similar in hydrologic characteristics as possible. They may 
have similar vegetal cover, land use, topographic conditions and geologic 
characteristics and the range of extrapolation of drainage area should not be 
excessive. They should also have similar rainfall and evapotranspiration regimes. 
The data pertaining to hydrology, topography, drainage network, soils and vegetation 
from several locations within a region are collected and analysed to derive empirical 
relationships for flood estimation for ungauged sites within each region. The 
procedures for data analysis and flood estimation may depend upon the length of 
record, the nature of variables, their accuracy and reliability and the amount of 


data available. 
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2.2 Factors Controlling Runoff 

The factors that control the magnitude of runoff in a region can be broadly 
classified into permanent factors and transient factors ( Fig. 2.1). The permanent 
factors include the basin characteristics such as the size, shape, slope etc.; drainage 
network characteristics like density, order, length of channel etc.; and channel 
characteristics such as crosssectional area, roughness and slope. The transient 
factors consist of storm characteristics such as its space and time distribution, 
evaporation, infiltration, and land use characteristics which, in turn, include man’s 
influence including afforestation, deforestation, cropping pattern, cultivation 
practices etc. This* classification is by no means exact because many of the factors 
are interdependent to a large extent and some of the factors may be permanent or 
transient. In order to estimate the flood magnitude, it is important to study the 
influence of various factors on the flood hydrograph and its peak. 

2.3 Classification of Reg iona 1 i sat i on Procedures 

Regional flood estimation procedures are used for ungauged sites or when the 
available data at the site are inadequate. A simplistic approach is generally to 
develop graphs or equations which may be used to estimate the peak discharge or 
flood hydrograph for any catchment in the region. Such empirical relationships 
should be used with care because they often depend on particular physical and 
climatic conditions not accounted for in the relationships and so they are applicable 
at best only to those regions for which they are developed . Depending on the 
purpose and availability of data, the procedures for data analysis and flood 
estimation may vary. 
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In many cases, an estimate of peak flow rate is sufficient for engineering 
design purposes, e.g., for design of bridges, culverts, etc. But for others, 
particularly those involving the provision of storage capacity, a complete flood 
hydrograph may be necessary. Hence, regionalisation procedures may broadly be 
classified into: 

a) Flood peak estimation procedures, and 

b) Flood hydrograph estimation procedures. 

They are briefly considered in the following sections. 

2.4 Flood Peak Estimation Procedures 

Many hydraulic design problems require simply an estimation of the peak 
flow rate in a stream, perhaps, for a specified frequency. Then, the general shape 
of the hydrograph and the time of occurrence of the peak flow are of no special 
significance. Approaches for estimating peak flood discharges fall into one of the 
two broad categories, viz., (i) deterministic methods and (ii) probabilistic methods. 
Fig. 2.2 shows the various approaches for peak flow estimation and they are 
considered in greater detail here. 

Peak Flow Estimation 


Deterministic Methods Probabilistic Methods 

i i 

i \ 1 f 1 

Enveloping Curves Rational Complex Rational Flood Frequency 

and Area Based Method Formulae Method Analysis 

formulae 


Fig. 2. 2 Approaches for flood peak estimation 
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2.4.1. Deterministic methods 

(i) Enveloping curves and area based methods 

The catchment area of a basin is hydrologically important because it directly 
affects the magnitude of flood peaks. The larger the size of the basin, the greater 
is the amount of rainfall intercepted and higher is the peak discharge that results. 
The observed maximum peak discharges for a number of catchments in a region are 
plotted on a log-log scale against the catchment area. The plotted points are then 
enveloped by a smooth curve referred to as the enveloping curve for observed 
maximum peak discharges. Two of the well known enveloping curves are those of 
Creager (Creager and Justin, 1950) and Meyer-Jarvis. Enveloping curves have been 
used as a guide for determining flood discharges in some projects in India like the 
Nagarjunasagar in Andhra Pradesh and Matatila in Uttar Pradesh (CWC, 1972). A 
typical enveloping curve as reported in CWC (1972) for Indian rivers is shown in 
Fig.2.3(a) and for different parts of the world (Gray, 1970) in Figs. 2.3(b) and 2.3(c). 
Mutreja(1986) has given details of enveloping curves for Romania and U.S.A. 

This method is based on available records and for maximum observed flood 
peaks only. The enveloping curves have to be reconstructed as more and more data 
become available since the maximum observed flood will be larger, the longer the 
period of record. It can hence be used only as a guide for the order of magnitude 
of the largest possible flood in a basin. 

A functional form may be fitted to the enveloping or other curves relating 
peak discharge to the area. These equations are generally of the form 

Q-CA" ' (2.1) 

where Q is the peak discharge, A is the basin area and C and n are constants. 
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Dickens in 1885 (CWC.1972) made the first attempt in India to derive a general 
formula for determining maximum flood discharge on the basis of studies for 
determining the waterway capacities for bridges. His formula is 

Q — CA 3/H (2.2) 

where Q is the discharge and A is the area. Mutreja(1986) presents a map of India 
which shows the values of Dicken’s constant C to be adopted in different parts of 
India. Chow (1964) and Ghosh (1986) have listed various area based formulae for 
flood peak estimation and Wolf (1966) gives an excellent review of early 
relationships. The manual on "Estimation of Design Flood - Recommended 

Procedures ” (CWC.1972) lists a wide variety of empirical formulae. Table 2.1 gives 
a number of area based formulae used in India and elsewhere in the world 
(Gray, 1970). 

The variation in the empirical formulae may be due to limited data, 

differences in curve fitting procedures adopted and because the data in each case 
are limited to a specific region. Hence such empirical relationships should be applied 
with care and to the specific region for which they may be relevant. 

(ii) Rational method 

The rational method assumes that the intensity of rainfall is an important 
parameter in the estimation of flood peaks and so estimates the peak runoff rates 
from data on rainfall intensity, catchment area and drainage characteristics. The 
method assumes that a rainfall of uniform intensity occurs for a sufficiently long 

time over the whole basin and losses and runoff are constant proportions of the 

rainfall. Furthermore, rainfall over different parts of the basin take different times 
to reach the gauging stations and so runoff will increase as the water from 



Table 2.1 Comparison of formulae for estimating flood discharge 

(Gray, 1970) 
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Table 2.1 contd... 


(1) 

(2) 

(3) 
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(5) 
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Table 2.1 contd... 
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more and more distant parts of the catchment reaches the outlet. When the whole 
area is contributing and the rainfall is uniform, a steady state is reached and 
discharge becomes a constant and maximum. Since the rainfall intensity generally 
decreases with duration, the maximum peak runoff will occur when the duration of 
the storm is around the time when the whole of the basin area starts to contribute 
to runoff at the gauging station. This is, in turn, estimated as the time of 
concentration of the basin given by Kirpich equation (1940), 

t c - 0.00032 L° 77 S ° 385 (2.3) 

where t c - time of concentration in hours, L - maximum length of travel in metres, 
and S - slope -= H/L where H is the difference in elevation excluding rapids and 
falls between the farthest point in the basin and the outlet in metres. 

The rational formula is given by 

Q - 0.278 C I A (2.4) 

where Q is the peak discharge in m 3 /sec, C is the coefficient of runoff depending 
on the characteristics of the drainage basin, I is the intensity of rainfall in 
mm/hour for a duration equal to the time of concentration and A is the drainage 
area in Km 2 . The choice of the value of C is the most critical factor in the use of 
the rational method. This coefficient has to account for various climatic 
conditions and physiographic characteristics of the basin and in particular for the 
characteristics of the design storm. In reality, it is not a constant and it varies 
from storm to storm in the same basin. The judgement required in estimating the 
value of C is, therefore, considerable. However, used with discretion, this method 
can provide a reasonable estimate of peak discharge. The value of C for different 
surfaces have been listed in Dunne and Leopold (1978) and Chow et al (1989). 


Chow (1964) suggests a modified form of rational formula for small basins 
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given by 

Q . ,.oo8 *42 <2.5, 

where Q is the peak runoff in cfs, Z is the fraction of unitgraph peak for a 
duration ’t’ to the equilibrium discharge for a continuous rainfall excess at a rate 
of l’Vhour over an area of A acres and R e is the rainfall excess in inches for a 
given duration of t hours for the basin . Comparing this with the rational formula, 
it is observed that the coefficient C is equal to 1.008 Z. 

(iii) Complex formulae 

The climate of a region often controls the soil and vegetation in a basin. It 
is related to precipitation and temperature or solar radiation. The rational formula 
explicitly considers only the area of the basin and the intensity of rainfall and all 
other factors are implicitly lumped in the factor C. In addition to the area and the 
intensity of rainfall, there are many other factors that control the magnitude of 
flood peak (Sec. 2.2). A large number of basin and storm factors have been 
incorporated into a wide range of flood formulae by means of standard statistical 
techniques producing expressions of the form 

Q " a Xj b X 2 C X 3 d (2 6) 

where Q is the peak discharge; a,b,c,d etc., are the regression coefficients; and 
X lf X 2 ... are the factors controlling the flood peaks. 

These procedures require data pertaining to various factors for several 
basins in the region. Chow(1964) and Wolf (1966) list formulae of this type. 
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2.4.2. Probabilistic approaches 

The various methods considered earlier for estimating flood peaks are 
deterministic. But the engineer or hydrologist is more often interested in the 
number of times a flood of given magnitude is equalled or exceeded in a specified 
period, i.e., he is more interested in estimating the flood magnitude with a given 
recurrence interval particularly for the analysis and design of medium and minor 
structures. In such a case probabilistic approaches are necessary for estimating 
the magnitude of flood peaks. 

(i) Rational Method 

The rational method can be considered as a probabilistic method. It is 
assumed that the rainfall intensity has a duration equal to the time of 
concentration of the catchment. The probability distribution of annual flood peaks 
under constant infiltration differs from that for infiltration which really varies 
from storm to storm (Ramaseshan, 1964). Yet as an approximation it is assumed that 
the average recurrence interval (ARI) of the computed maximum rate of discharge 
for the* catchment under constant infiltration condition is the same as the ARI for 
the corresponding storm. Thus if it is required to calculate the flood peak with a 
specified recurrence interval, a storm of the specified recurrence interval and 
duration t c is used in the calculation. Thus, when used as a probabilistic model, 
the rational formula takes the form (Pilgrim, 1989) 

Q(Y) - 0.278 CCY) I(tc,Y) A (2.7) 

where Q is the discharge and I is the intensity of rainfall for a duration t c with an 
average recurrence interval of Y years. C, being a function of Y, is not a 
constant. This formula converts a rainfall of ARI of Y years derived from 
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intensity-duration-frequency data for the region under study into a peak discharge 
with the same ARI. Both the design rainfall and peak discharge are derived from 
frequency analysis and hence this formula has nothing to do with a particular 
storm. Pilgrim (1989) used this approach for Australia and has found that it fulfills 
the requirments for design methods for small and medium catchments. 

(ii) Regional flood frequency analysis 

Estimates of likelihood of flood of a given magnitude occurring within a 
stated interval of time are commonly required for planning and design of hydraulic 
structures. Where stream gauging records of sufficient length are available at or 
near the point of interest, a frequency analysis of the historical floods is useful 
in estimating the likelihood of future floods with a given probability of 
occurrence. The flood parameters for which estimates of frequency are required 
may be the peak stage, peak discharge or volume of runoff. 

Approaches to regionalisation are generally based on the assumptions that 
the frequency distribution applicable to flood peaks is unique, say Extreme Value, 
Log Pearson Type III etc., and that the parameters of the distribution, say, mean 
and standard deviation are related to the area of the basin and parameters like C v 
( Coefficient of variation) and C s (Coefficient of skew) may vary in space. 
Furthermore the parameters may depend on the length of available data. 

Estimates of flood frequency are often required at ungauged sites for which 
no data are available. Regionalisation is useful in making such estimates. The 
record at any gauged site may not be of sufficient length to provide reasonably 
accurate parameter estimates. In such cases, information from nearby sites in the 
region may be correlated with the record of the site with data of inadequate 
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length to increase its information and to provide a regionalised estimate of flood 
frequency. 

i Le index flood method devised by the United States Geological Survey and 
reported by Dalrymple (1960) ia a method of regionalisation of flood estimation . It 
attempts to extrapolate statistical information of runoff events for flood 
frequency from gauged catchments to ungauged catchments in the vicinity having 
similar catchment and hydrologic characteristics. 

To illustrate the empirical nature of such approaches, the various steps of 
the index flood method are briefly listed: 

i) Select gauged catchments within the region having characteristics similar 
to the ungauged catchments. 

ii) Determine the time base period to be used in the study; 

iii) Establish flood frequency curves for each gauging site using Gumbel 
Extreme Value-Type I distribution; 

iv) Estimate the mean annual flood Q 2.33 for each station; 

v) Test homogeneity of data by a specific procedure; 

vi) Establish relationship of mean annual flood and catchment 
characteristics, usually drainage area, at each station; 

vii) Rank ratios of selected return period floods to the mean annual flood 
at each station; and 

viii) Compute median flood ratio for each of the selected return periods; 
multiply by the estimated mean annual flood of the ungauged catchment; 
and plot against the recurrence interval on Gumbel’s probability paper 
to get the flood frequency curve for an ungauged catchment. 

The regional flood frequency analysis for British Isles reported in the Flood 
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ood Studies Report (NERC.1975) employed a modification of Dalrymple’s approach 
inparable to that of Benson(1962) (Hall, 1981). 

The usefulness of frequency analysis is restricted because observed flow 
.ta for small catchments are usually very limited (Raudkivi,1979). Apart from the 
gional flood frequency analysis a large number of studies using statistical 
iproaches have been carried out. Condie et al (1987) used the analysis of variance 
ichnique to compare the regional flood frequency methods in Southern Ontario, 
heir study allows an assessment of the influence of drainage area and return 
sriod and they found that the drainage basin size has a significant effect on the 
ssult of the study. Stedinger and Tasker (1985) proposed a method of estimating 
le parameters of regional regression models that take into account their length of 
icord available at gauged sites and the between-site correlation among concurrent 
ows. This method uses the generalised least square estimator of the regression 
odel parameters. Several investigators like Shane and Garver (1970); Vicens, 
odriguez and Shaake(1975); and Kuczera(1982) have addressed the issue of how 
est to combine at-site and regional regression estimates. Mimikou and Gordio(1989) 
sed the multiple regression techniques in the development of regional relationships 
etween the basin characteristics and the mean annual floods. They successfully 
sed the model in predicting mean annual floods needed in the hydrologic design for 
ngauged catchments for the study region in Greece. Many more studies have been 
arried out for the regional flood frequency analysis ( Singh,1987). 

1.5 Flood Hydrograph Estimation Procedures 


The time distribution of runoff during a storm is often needed for many 
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planning purposes. A hydrograph depicts the time distribution of runoff. It is used 
for engineering planning such as reservoir design, for assessing the flood damage 
potential and the influence of flood detention structures in reducing the flood 
peaks. The whole hydrograph is necessary when hydrographs from dissimilar 
tributary areas are added and routed downstream to a channel reach of interest. 
The duration of flooding, a critical factor in many planning problems, can also be 
studied from a knowledge of the hydrograph. 

2.5.1 Ra i nf a 1 1 - runof f relationships 

From the hydrologic point of view, runoff from a drainage basin may be 
considered as an output from the watershed due to effective rainfall as input. The 
runoff is influenced by both physiographic and climatic factors. It is hence often 
necessary in hydrologic analysis and design to develop relationships between 
rainfall and runoff possibly using some of the factors affecting runoff. These 
relationships differ with the type of precipitation, consideration of the volume or 
peak runoff, or the time distribution of runoff. The need for such cause-effect 
relationships has lead to the development of hydrologic models. Such models help 
to predict the hydrologic output, for example, the runoff hydrograph from a given 
rainfall. 

From available data, regression relationships are derived to relate peak 
discharge to hydrometeorological features of the storm and physiographic features 
of the watershed. Even when adequate data are available for establishing such 
relationships, they are limited to regions for which such relationships have been 
derived. Since generally adequate data are available for estimating storms, it 
seems worthwhile to derive relationships between storm rainfall and streamflow. 
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without bothering about the processes governing the water movement. Such models 
where no consideration is given to the internal structure of the system end where 
the model parameters and variables are lumped in space are termed as lumped 
parameter models (Eagleson,1967; Holton, 1970). The development of lumped models 
can be stated to have started with the advent of the rational formula. Many 
models of the rainfall-runoff process treat the precipitation input as uniform over 
the watershed and ignore the internal spatial variation of the watershed flow. 
During 1960’s and 70’s, a large number of conceptual lumped models have been 
developed (Fleming, 1975). 

A unit hydrograph representation is a classical example of lumped rainfall 
runoff system and since it was first introduced by Sherman (1932), in many 
situations it continues to be the most practical tool available to hydrologists for 
making flood estimates. 

The unit hydrograph of a catchment is defined as the direct surface runoff 
(DSRO) resulting from unit excess rainfall over the drainage area generated at a 
constant rate during a specified duration. In the estimation of flood discharges 
from rainfall intensity and duration, the first concern is to find the amount of 
rainfall excess and the second is to obtain the flood discharge from the rainfall 
excess. With sufficient rainfall-runoff data, UH can be derived by various methods 
(Clarke, 1945; O’ Kelly, 1955; Nash.l 957; Dooge,1959 ). 

a) Synthetic unit hydrograph 

To derive a UH for sites where data are not available, it is important to 
have procedures by which UH can be constructed. For ungauged basins, the UH is 
synthesised from the physical characteristics of the basin. 
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Estimates of UH for an ungauged basin can be made if the information 
derived from the analysis of gauging records on other streams is regionalised. UHs 
are computed for the gauged streams in an area and the lag, peak and duration of 
these hydrographs are related to the hydrometeorological parameters of the 
catchment such as drainage area, channel gradient and drainage density. These 
physiographic and geomorphologic characteristics of the catchment represent the 
various factors that affect the storage and transmission of the temporally and 
spatially distributed runoff generated by a rainstorm and thereby control the time 
distribution of runoff. A set of regression equations is derived between the 
characteristics of the UH and geomorphological variables. Knowing the values of 
the geomorphological features of the basin from physical measurements and 
meteorological characteristics of the storm, flood peaks from ungauged basins can 
be estimated using regression equations. Several methods (Snyder, 1938; 
Commons, 1942; William, 1945; Taylor and Schwarz, 1952; SCS, 1957; Bender and 
Roberson, 1961; Gray, 1961) have been developed for synthesis of UH. 

i) Soil Conservation Service triangular hydrographs . From a 
study of UH for a large number of small drainage basins, SCS (1972) has 
developed modification of the synthetic unit hydrograph technique. The simple 
method is an approximation of the UH by a triangular unit hydrograph as shown in 
Fig .2.4. 
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From the area of the triangle, the volume of runoff V is equal to 

V - iQpt. 

— ^ (tp + t«) 

- Qp " £ru 


( 2 - 8 ) 

(2.9) 

( 2 . 10 ) 


Examination of a large number of hydrographs from small agricultural watersheds 
throughout the United States of America lead to the empirical generalisation that 


t. - 1.67 tp 

which gives peak discharge, Q p 


( 211 ) 


— , and time to peak t P - T/2 + t,_. 
tp 


ii) SCS dimensionless hydrograph. The SCS has also used a 
dimensionless hydrograph approach. The dimensionless hydrograph is a synthetic 
unit hydrograph in which the discharge is expressed by the ratio of discharge Q 
and peak discharge Q p and the time by the ratio of time T to the time of rise of 
UH t p . Given the peak discharge and the lag time for the duration of excess 
rainfall, the UH can be estimated from the synthetic dimensionless hydrograph of 


the basin. 
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The SCS triangular and dimensionless UH’s employ a standard hydrograph 
form to all the catchments. If one does not want to rely on such a generalised 
procedure, he may adopt the SCS approach but derive dimensionless or triangular 
UH for his own region. 

iii) Snyder’s synthetic unit hydrograph. Snyder introduced the 
procedure for developing the UH for basin with no data. He correlated the timing 
and peak rates of hydrographs derived from basins in the Appalachian mountains 
with measures of their physiographic characteristics. He used three parameters to 
describe the hydrograph, viz., the time lag to peak, t p ; peak discharge, Q p; and the 
baselength t B . Snyder’s equation in metric units for a rainfall excess of standard 
duration t r - t p /5.5 are as follows: 


tp 

0.75 Q (LLc)° 3 

hours 

(2.12) 

Op - 

2.75 Cp A 
tp 

m 3 /sec, and 

(2.13) 

te - 

3 + 3 (tp/24) 

days 

(2.14) 


where L is the length of the main stream in Km. from the outlet to the stream 
divide; L c is the distance in Km. from the outlet to a point on the stream nearest 
to the centroid of the catchment; and C t and C p are constant coefficients derived 
from the gauged watersheds in the same region. For a different duration t, the lag 
time was adjusted by 

tpR “ tp + Ct - tr)/4 (2.15) 

Q pr and t BR are calculated using this value of t PR . With these values, the shape of 
the hydrograph may be sketched. The values of Snyder’s coefficients C t and C p 
have been found to vary considerably depending on the topography, geology and 
climate. Although Snyder indicated that the coefficient C t is affected by basin 
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2.272 (L Lc/S) ° 408 

(2.17) 

tp 


2.167 (q P )° 940 

(2.18) 

W 60 


2.084 (qp)' 1 065 

(2.19) 

W R 60 


0.856 (q P )-° 865 

(2.20) 

^R76 


0.440 (q P )-° 918 

(2.21) 



5.428 (t p ) 

(2.22) 

tm 


tp + t r /2 

(2.23) 

Qp 


“ q P X A 

(2.24) 

where q p 

- 

Peak discharge of UH/unit area in m 3 /sec; 


tp 

- 

Time to peak of UH from centre of unit rainfall duration, 

in 



hrs; 


W 60 

- 

Width of UH measured at 50% peak discharge ordinate, in 

hrs; 

w 75 

- 

Width of UH measured at 75% peak discharge ordinate, in 

hrs; 

Wrso 

- 

Width of rising side of UH measured at 50% peak discharge 
ordinate, in hrs; 

w Rr5 

- 

Width of rising side of UH measured at 75% peak discharge 
ordinate, in hrs; 


- 

Base width of UH , in hrs; 


tm 

- 

Time to start of rise to the peak of UH, in hrs; 


Qp 

- 

Peak discharge of UH in m 3 /sec; and 


S 

- 

Equivalent stream slope in m/km. 



Similar relationships have also been derived for other subzones in the country. 

b) Use of conceptual models 

Various methods are available for deriving unit hydrographs (subec.2.5.2). In 
order that the UH can be regionalised, it is necessary to represent the UH in 



terms of a limited number of parameters, perhaps in terms of conceptual model 
very useful and simple model for UH is the one proposed by Nash (1957). Sti 
have been carried out (Singh, 1976; NIH,1985) to compare the various ccmcep 
models and it has been found that the agreement between the observed 
predicted hydrographs was quite good for the Nash model vis a vis the oi 
models. Nash model has two parameters N and K. Nash (1960) applied his model 
British catchments. He showed that the two parameters of his model were rela 
to the first and second moments of the IUH about the origin as 


n»i 

NK ; and 

(2. 

*** 

1/N 

(2.1 

These two moments were than correlated 

characteristics as 

empirically with the watersh 

“i 

27.6 A 0 3 S,, 0 - 2 ; and 

(2-2 

m 2 — 

0.41 V 01 

(2-2! 


where S 0 is the overland slope in parts per 10,000 calculated as mean of a gri 
sample of slopes, A is the watershed area in square miles and L is the length o 
the channel in miles. 

Two approaches are possible for using the UH method. One can study th< 
storm in a given catchment and the resulting flood hydrographs, isolate or derive 
the UH for the catchment and then use it for estimation of design floods on the 
same catchment. In the second approach, UH from various catchments are taken and 
an attempt is made to correlate the hydrograph parameters with the' storm and 
catchment characteristics. 

In the present study, though many methods are available for the UH 
approach, it was decided to use the Nash model because it is simple and a number 



of computer programs have already been implemented at IIT Kanpur- This mode) 
also well established and has been found to be efficient. Hence this model is uf 
for the derivation of unit hydrographs, for regionalisation of the model paramet) 
and hence for flood estimation for small catchments in the region. A bri 
description of the Nash model is presented in Chapter 4. 

2.5.2 Distributed parameter models 

A hydrologic system can be called distributed if its parameters and variabh 
are distributed in space (Eagleson,1967; Nash, 1969; and Holton, 1970). Distribute 
parameter models are physically based and the hydrologic processes ar 
represented by a set of partial differential equations, interrelated by the concept 
of continuity of mass and momentum. Freeze and Harlan (1969) suggested ( 
blueprint for physically based models. This blueprint suggests the types of input) 
and outputs required for such physically based hydrologic response models. Th« 
output from such a response model would provide a total picture of the hydrologic 
system. This output would include the streamflow hydrographs at any point within 
the basin, the groundwater flow pattern and the soil moisture regime. The output 
would also incorporate surface water, soil water and groundwater zones as 
components in a single system, not as discrete elements. The blueprint suggested by 
freeze and Harlan forms the basis for many physically based distributed parameter 
models. 

In the distributed parameter models, the catchment is divided into a number 
of segments and each segment is simulated separately and the results combined to 
obtain the catchment response. Increasing the number of segments allows a more 
detailed description of the basin, but at the same time increases the computer 



processing time in direct proportion to the number of segments. Moreover witl 
multiple segments, determination of some parameters may become somewhat raor 
difficult. Computation with details of the distributed parameter model is usually 
time consuming. However such a model is useful for research on the consequent 
of heterogeneity in a catchment. It is appropriate to ask under what conditions anc 
for what type of problems is it advantageous to use a distributed parameter model 
The choice of a model for a particular problem is never a simple one. It is 
inevitably based on economic constraints and personal experience and constraints, 
as well as purely hydrological consideration and scientific vigour. Data availability 
is often a crucial factor in such decision making (Beven, 1985). 

Beven and O’Connel (1982) have considered the role of distributed parameter 
models in hydrology in some depth. Four major areas which offer the greatest 
potential for the application of distributed models have been identified by them as: 

i) forecasting the effect of land use changes; 

ii) forecasting effects of spatially variable inputs and outputs; 

iii) forecasting the movements of pollutants and sediments; and 

iv) forecasting the hydrologic response of ungauged catchments where no data 
are available for calibration of lumped models. 

Thus, in particular, the following problems may need distributed parameter models: 

i) Hydrology of urban and small watersheds for runoff estimation; 

ii) Groundwater modelling; 

iii) Estimation of soil erosion, transport, deposition and sediment behaviour; 

iv) Surface water - groundwater interaction; 

v) Water quality simulation; and 

vi) Human impact on environmental processes dealing with the above. 



All of these problems allow distributed models to exploit their majo 
advantages over other models. The parameters of the model are physicall: 
based and they can generally be either measured in the field or derived froi 
field measurements. In addition, the parameters as well as input data an< 
derived variables are distributed over space and time. The physical basis ol 
the model parameters also give some hope that the measured parameter values 
may be extrapolated to other locations or other time periods. 

The development of distributed parameter modelling of catchment 
hydrology has been a slow faltering process (Beven, 1985). There have been 
many studies on modelling individual processes, especially groundwater flows, 
unsaturated soil water flow and channel routing; but models involving 
interactive processes and the application of catchment scale models to real 
world problems are limited. Examples are those of Betson (1979), Dagan (1979), 
Freeze (1971,1980), Jayawardhane and White (1977,1979), Lane and Woolhiser 
(1977), Morris and Woolhiser (1980), Narasimhan and Witherspoon (1977) etc. 
Fleming (1975) lists some of the distributed parameter models that were 
developed in the sixties and early seventies. Dawdy et al (1978) developed the 
USGS Distributed Routing Rainfall Runoff model (DR3M). The computer program 
for the model for routing urban flood discharges through the branched system 
of pipes or natural channels using rainfall as input has been well documented. 
The model combines the soil moisture accounting and rainfall excess components 
developed by Dawdy et al (1973) with the kinematic wave routing method 
presented by Leclerc and Schaake (1973). Other models that have been 
developed include System Hydrologique European (SHE) developed in Denmark, 
Sweden and Britain (Beven,1980), USD A Agricultural Research Service Small 



Watershed Model and the Institute of Hydrology Distributed Model (IHDh 
Morris, 1980). There are in addition a number of models available that treat onl 
interacting infiltration and surface flow process (e.g., Huggins and Monke,1961 
Smith and Woolhiser, 1971; Engman and Rogowski,1974; Ross et a 
1979;Kutchment, 1980). Beven and Kirby (1979) developed a physically base 
model which had certain parameters that are measurable in the field. 

Predicting the response of ungauged catchments is one area o! 
application where a physically based model can play an important role, sinc< 
parameter values may be measured or estimated on the basis of more readily 
available information. There have been relatively few studies where such 
models have been applied to catchments, treated as if they were ungauged, with 
the prediction later compared to available measurements. Beven (1985) reports 
on the Finite Element Storm Hydrograph Model (FESAM; Li et al, 1977; Ross et 
al, 1979; and Shanholz, 1981) which is based on a Hortonian infiltration excess 
concept of runoff generation. 

2.6 Models Selected for the Study 

The present study is essentially a regional case study for flood 
estimation in small catchments. The data available for the study (Chapter 3) 
are very limited and hence regional flood frequency studies are not possible. 
Though the regression and empirical models of CWC and RDSO are available, 
they are not used in the present study because these relationships have 
essentially been developed for large catchments which have adequate and 
appropriate date. Since generally adequate data are available for estimating 
design storms, it seems worthwhile to derive relationships between storm 
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rainfall and streamflow, perhaps in terms of UH and then regionalise the model 
parameters for flood estimation for small catchments in the region. A well 
tested and established UH model, the Nash conceptual model is used in the 
present study ( Chapter 4 ). 

A more realistic approach would be to model the catchment runoff 
process according to physical principles. Information regarding physiography, 
soils, vegetation and the drainage pattern of the catchment is used in the 
physically based models. Though many other models are also available 
(Yen, 1989) which are physically based (Subsec. 2.5.2), a well tested model, the 
USGS DR3M is used in the present work to study the regional variation, if any, 
of the parameters of the model and also to judge their validity for flood 
estimation (Chapter 5). 



CHAPTER 3 


MODEL SELECTION AND PRELIMINARY ANALYSIS OF DA 


3.1 Introduction 

Since there are a large number of models generally available for hydrologic 
studies (Chapter 2), it is necessary to choose one or more of them for detailed 
analysis. This involves both the choice between competing models and procedure for 
estimation of optimum value of the parameters of the models chosen. The choice of 
a particular model depends on many factors like the purpose, data availability and 
the ease with which one can use them. 

3.1.1 Criteria for model selection 

World Meteorological Organisation Commission on Hydrology (WMO,1975) has 
recommended the following criteria for selection of models in various hydrological 
forecasting situations: 

i) General and specific purpose and benefits of forecast (eg., continuous 
hydrograph of flood discharge, water quality etc;) 

ii) Climate and physiographic characteristics of the basin; 

iii) Length and records of various types of data; 

iv) Quality of data both in space and time; 

v) The availability and size of the computer, both for development and 
operation of the model as well as possible use of model by relatively non- 
expert hydrological forecasting personnel. 
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vi) The ability of the model to be conveniently updated on the basis of 
current hydrometeorological conditions. 

The models chosen for the present study are the Nash model and the USGS 
DR3M. These models were chosen with a specific purpose. They are both easy to 
use and the computer programs for these models were available in the computer 
system used for the study. The parameters of these models are used for 
regionalisation for flood estimation. Since the data available are for a limited period, 
the regional flood frequency analysis or any other statistical approaches cannot be 
used. With the available data, the two models chosen can be applied to the various 
basins in the study region. Though many computer systems were available for the 
study, the entire computational work was carried out on the microVAX II system of 
the Centre for Water Resources Engineering and Management, IIT Kanpur. The 
model parameters can be updated with the availability of additional data. 

3.1.2 Steps in model selection 

Dooge (1973) gave a rational methodology for the selection, calibration and 
use of mathematical model for the runoff process that might follow a procedure as 
shown in Fig. 3.1. The first step is to define a problem (Chapter 1). The next step 
involves choosing class of models and selecting a model. The present study uses a 
lumped parameter model and a physically based distributed parameter model . The 
models chosen are the Nash model and the USGS DR3M. The process is explained in 
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Chapter 2 and 3. The models are used for the various basins with the available set 
of data and their performance is evaluated and used in flood estimation (Chapters 4 
and 5). The results are summarised in Chapter 6 along with the discussions on the 
suitability of the selected models for the given problem. Thus, in general, the 
procedure followed in this study follows the method suggested by Dooge. 

3.2 Study Area 

The various basins used in a regional study have to be generally similar in 
hydrometeorological characteristics. However, it is usually very difficult to specify 
the criteria for regional homogeneity even though some empirical tests have been 
suggested (Dalrymple,1960). The region selected for the study may be based on 
empirical criteria. For purposes of hydrometeorological studies, India has been 
divided by the Central Water Commission (CWC) into 26 subzones, each of which is 
assumed to be homogeneous and these are shown in Fig. 3.2. The present work is 
essentially a regional case study for flood estimation, especially in small catchments. 
RDSO of the Ministry of Railways has collected and in collaboration with CWC has 
analysed a large amount of hydrometeorological data for floods for a large number 
of catchments draining into various railway culverts and bridges in different 
subzones. After discussion with the officers of the RDSO, a part of the North 
Brahmaputra basin, designated as subzone 2(a) was identified for this study. It may 
be noted that the data available for small catchments in this region are yet to be 
analysed and interpreted. 

The subzone lies approximately between longitudes 88°E and 97° 20’ E and 
latitudes 26°N and 29°25’N. Actually the various small basins for which data were 
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available and are used in the present study lie in a part of this subzone between 
longitudes 89°E and 93°E and latitudes 26°N and 27°N (Fig. 3.3). 

The Brahmaputra river with a total catchment area of 0.94 million Km 2 drains 
one of the largest river basins in the world. The total length of the river in India 
is 885 km. The drainage area of the river in India is 1,95,000 Km 2 . The drainage 
area of the North Brahmaputra subzone 2(a) is 1,21,444 Km 2 and the South 
Brahmaputra subzone 2(b) is 73,556 Km 2 . The topography of the area is generally 
hilly with rivers having moderate steep slope, meandering channels, beds and bank of 
alluvial soils and comparatively low silt charge. The south-west monsoon and 
cyclonic storms cause rainfall in the subzone from May to October. Normal annual 
rainfall in the basin varies from 2000 mm to 5000 mm. The subzone is covered with 
a variety of soils. Broadly they can be classified as red loamy soil, brown hill soil, 
Terai soil and alluvial soil of recent origin. The subzone has a considerable area 
under forests which may have undergone changes in the recent times because of 
human activities. 

3.3 Data Availability 
3.3.1 Available data 

The data available for various river basins in the subzone include: 

i) Catchment plan; 

ii) Crosssection at bridge site; 

iii) Longitudinal section maps for some of the catchment selected; 

iv) Hourly rainfall data from raingauges at one or more self recording 
raingauge stations for the period, generally from middle of May to middle 
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of October; 

v) Hourly gauge data at bridge site; and 

vi) Data for rating curves in terms of discharge at gauging sites through 
velocity observations, using generally floats, two to three times a day; 

3.3.2 Source of data 

North-East Frontier Railway under the supervision and guidance of RDSO had 
observed and collected hydrometeorological data for the subzone 2(a). Flood 
estimation procedures for the larger catchments have been developed using the unit 
hydrograph principle (CWC.1991). Since data and procedures for the large basins are 
generally available, RDSO is interested in the development of a procedure or 
procedures for estimation of design flood for railway bridges and culverts, 
particularly for small catchments. Some information was also obtained from the 
Small Catchment Directorate of CWC. Meteorological data were gathered from the 
IMD and other publications, e.g.. Irrigation Atlas of India; Report of the National 
Commission of Floods (RBA, 1980) etc. 

3.3.3 River basins selected for the study 

The entire subzone, being in a very difficult terrain, had disadvantages of 
inadequate accessibility and communication. In general the data for the small 
catchments in the region pertain to the period 1961 to 1966. During this period, 
there were two major wars and there were no observations. Hence, this area 
suffered from disturbed conditions during the period of observation. Due to these 
reasons, collection of data and supervision of the same could not be of adequate 
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standard. 

Based on the available information and after initial scrutiny of data, the 
following river basins were identified for use in the present study: 


Table 3.1 River basins identified for the study 


Sl.No 

River Basin 

Bridge No. 

Catchment 
Area (Km 2 ) 

No.of Rain Theissen 
gauges Weights 

1 . 

Batjhora 

173 

7.33 

1 

1.00 

2. 

Sankhini 

112 

7.76 

1 

1.00 

3. 

Kurti 

139 

21.91 

1 

1.00 

4. 

Mujnai 

201 

38.50 

2 

0.24; 0.76 

5. 

Talma 

24 

42.12 

1 

1.00 

6. 

Khaga 

114 

66.00 

2 

0.46; 0.54 

7. 

Dharsi 

285 

91.40 

2 

0.23; 0.77 

8. 

Mansari 

22 

213.00 

2 

0.24; 0.76 

9. 

Gaghra 

506 

13.02 

1 

1.00 

10. 

Balason 

8 

350.20 

4 

0.11;0.27;0.38tf).24 


3.4 Preliminary Analysis 

For both the models used in the study, rainfall that causes the flood and the 
corresponding streamflow values are required in addition to the catchment details. A 
close scrutiny of the available data indicated inconsistencies and errors and these 
have to be taken care of before using the data in the two models. The following 
sections briefly explain the various preliminary analyses carried out and the 
decisions taken. 


3.4.1 Precipitation data 


Precipitation data were available at one or more locations in each basin. The 
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data were, in general, measured using self recording raingauges and tabulated hourly 
rainfall values were made available for the study. For 4 basins, there was only one 
raingauge station. In other cases, where there were more than one raingauge in a 
basin, the Theissen polygon method was used to estimate the average rainfall over 
the basin. The average precipitation over the catchment is calculated as: 

Pmv “ X 1 Pl + ^ Pl X Pn <3 * 1) 

where P av - Average precipitation over the catchment; 

A^Aa A n - Area within each Theissen polygon with corresponding 

rainfalls Pn P 2 P n ; and 

A - Total area of the catchment. 

The Theissen weights for all the basins are given in Table 3.1. 

3.4.2 Stage data 

A detailed study of the available data indicated large errors and 
inconsistencies. The errors were mainly in the stage - discharge observations. The 
flood events were identified based on the rise and fall in the stage of the river at 
the gauging site. On many days, night observations were not available and even if 
they were available, they showed significant inconsistencies. For the computation of 
discharges, velocities were observed using floats and hence were widely varying. In 
a single day, there were widely diff erring values of velocities for the same value of 
stage and vice versa. Further in a number of cases, on days with heavy rainfall, 
the observation sites, as per remarks in the record books, were inaccessible and 
hence stage data were not available. There were also remarks in the observation 
book, such as the instruments being tramped by animals like rhinoceros and 
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elephants. Thus very valuable data were lost for many days till the instruments got 
replaced. 

After detailed scrutiny of records, including comparison of rainfall and 
runoff data, fairly consistent flood events were identified for the study. 

3.4.4 Rating curve 

It is difficult to measure continuously streamflows in rivers, but relatively 
easier to measure the stages or water levels continuously. Thus the establishment 
of a relationship between the rate of flow and water levels (and slopes if necessary) 
in channel, river etc., helps in conversion of stage data to discharge data. Water 
level in the river and corresponding discharge may be hydraulically related for 
uniform flow through Manning’s equation or empirically related by plotting the 
measured discharge against the corresponding stage and establishing a regression 
relationship between stage and discharge. This relationship between stage and 
discharge is referred to as a rating curve and the relationship may be of the form 

Q - K ( H - HoT (3.2) 

where K and n are the constants, and H is the stage while H 0 is the hypothetical 
datum corresponding to zero discharge. 

While establishing the rating curve, it is to be noted that the discharge is 
generally expressed as a function of stage alone and corrections for slope are made 
where necessary. Logarithmic transformation of Eqn.3.2 leads to a linear relationship 

log Q - log K + n log (H - Ho) (33) 

In the present study, the electronic spreadsheet software Lotus 1-2-3, available for 
developing the rating curve (Tirupati,1989), is used. A worksheet is prepared using 
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stage and discharge data. For an assumed value of H 0 , the built-in regression option 
is used for the estimation of regression parameters. The regression output in the 
worksheet gives the additive and multiplicative constants, the standard error of 
estimate etc. A sensitivity analysis is carried out to determine by the method of 
least squares, estimate of H 0 and the corresponding values of K and n., i.e., when 
the standard error of the Q e8 timate is minimum. The corresponding level H 0 is chosen 
as the datum level. A 95% confidence band is also computed. The logarithmic values 
may be transformed back to natural discharge values. Using the graph utility of 
Lotus 1-2-3 , rating curve for the basin with 95% confidence band are plotted in 
logarithmic and discharge domain. 

The rating curves are prepared for all the basins selected for the study. A 
typical worksheet is shown in Table 3.2. Titles and column headings are given in the 
top of the worksheet. The first two columns in the worksheet have the gauge and 
discharge data. The third column contains the difference between the gauge and 
datum level. The datum level selected is given in column 13. Columns 4 and 5 give 
the logarithm of the (stage - datum) and discharge values. A regression is carried 
out between the values of columns 4 and 5 with stage as the independent variable 
and the discharge as the dependent variable. The regression outputs are presented 
starting from column 12. Column 6 gives the values of the discharge calculated 
using the results of the regression. Column 7 and 8 have the 95% confidence band 
values calculated as equal to Q #atlmat# rt:l*96 * standard error 



TABLE 3.2 RATING CURVE FOR HUJNAI BASIN FOR THE YEAR 1964 
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CHARGE STAGE- LN(H-HO) LN(Q) Y coapute 95X confidence Q compute 95Xconf idence 
SECS DATUM X Y level for Y level for Q 

2 3 4 5 6 7 8 9 10 11 12 13 


476 

475 

454 

521 

494 

540 

510 

460 

448 

482 

496 

518 

503 

613 

588 

543 

578 

586 

607 

548 

652 

660 

656 

684 

833 

709 

696 

677 

619 

636 

681 

666 

707 

692 

762 


1.2 0.182321 6.165417 6.201486 6.355659 6.047313 493.4821 575.7421 422.9751 Datum level for which lininua 
1.2 0.182321 6.163314 6.201486 6.355659 6.047313 493.4821 575.7421 422.9751 Error in Y Estimate occurs: 296.5 

1.2 0.182321 6.118097 6.201486 6.355659 6.047313 493.4821 575.7421 422.9751 

1.2 0.182321 6.255750 6.201486 6.355659 6.047313 493.4821 575.7421 422.9751 Regression Output: 


1.2 0.182321 6.202535 6.201486 6.355659 6.047313 493.4821 575.7421 422.9751 Constant 5.976263 

1.3 0.262364 6.291569 6.300363 6.454537 6.146190 544.7700 635.5793 466.9353 Std Err of Y Est 0.078659 

1.3 0.262364 6.234410 6.300363 6.454537 6.146190 544.7700 635.5793 466.9353 R Squared 0.969516 

1.3 0.262364 6.131226 6.300363 6.454537 6.146190 544.7700 635.5793 466.9353 No. of Observations 113 

1.3 0.262364 6.104793 6.300363 6.454537 6.146190 544.7700 635.5793 466.9353 Degrees of Freedom 111 

1.3 0.282364 6.177944 6.300363 6.454537 6.146190 544.7700 635.5793 466.9353 


1.3 0.262364 6.206575 6.300363 6.454537 6.146190 544.7700 635.5793 466.9353 X Coefficient(s) 1.235305 

1.3 0.262364 6.249975 6.300363 6.454537 6.146190 544.7700 635.5793 466.9353 Std Err of Coef. 0.020790 


1.3 0.262364 6.220590 6.300363 6.454537 6.146190 544.7700 635.5793 466.9353 

1.4 0.336472 6.418364 6.391909 6.546082 6.237736 596.9956 696.5105 511.6990 

1.4 0.336472 6.376726 6.391909 6.546082 6.237736 596.9956 696.5105 511.6990 

1.4 0.336472 6.297109 6.391909 6.546082 6.237736 596.9958 696.5105 511.6990 HO Std Err Of 

1.4 0.336472 6.359573 6.391909 6.546082 6.237736 596.9956 696.5105 511.6990 Y estiaate 

1.4 0.336472 6.373319 6.391909 6.546082 6.237736 598.9956 696.5105 511.6990 296.00 0.082661 

1.4 0.336472 6.408528 6.391909 6.546082 6.237736 596.9956 696.5105 511.6990 296.50 0.078659 

1.4 0.336472 6.306275 6.391909 6.546082 6.237736 596.9956 696.5105 511.6990 297.00 0.084321 


1.5 0.405465 6.480044 6.477137 6.631310 6.322963 650.1070 758.4753 557.2221 

1.5 0.405465 6.492239 6.477137 6.631310 6.322963 650.1070 758.4753 557.2221 

1.5 0.405465 6.486160 6.477137 6.631310 6.322963 650.1070 758.4753 557.2221 

1.5 0.405465 6.527957 6.477137 6.631310 6.322963 850.1070 758.4753 557.2221 

1.5 0.405465 6.450470 6.477137 6.631310 6.322963 650.1070 758.4753 557.2221 

1.5 0.405465 6.563855 6.477137 6.631310 6.322963 650.1070 758.4753 557.2221 

1.5 0.405465 6.545349 8.477137 6.631310 6.322963 650.1070 758.4753 557.2221 

1.5 0.405465 6.517671 6.477137 6.631310 6.322963 650.1070 758.4753 557.2221 

1.5 0.405465 6.428105 6.477137 6.631310 6.322963 650.1070 758.4753 557.2221 

1.5 0.405465 6.455198 6.477137 6.631310 6.322963 650.1070 758.4753 557.2221 

1.5 0.405465 8.523562 6.477137 6.631310 6.322963 650.1070 758.4753 557.2221 

1.5 0.405465 6.501289 6.477137 6.631310 6.322963 650.1070 758.4753 557.2221 

1.5 0.405465 6.561030 6.477137 6.631310 6.322963 650.1070 758.4753 557.2221 

1.6 0.470003 6.539585 6.556861 6.711035 6.402688 704.0587 821.4203 603.4654 

1.6 0.470003 6.635946 6.556861 6.711035 6.402688 704.0587 821.4203 603.4654 
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of estimate. Columns 9,10 and 11 have the retransformed values of the discharge 
and the 95% confidence level values. A few iterations are carried out by changing 
the datum level and the value of H 0 corresponding to the lowest value of the 
standard error of Q es timate is chosen for preparing the rating curve. A plot of 
(stage - datum) against Q e9t imate and the 95% confidence band gives the rating curve 
for the basin with 95% confidence. A typical rating curve for Mujnai basin is 
shown in Fig. 3.4. (a). It is observed that a single rating curve does not fit the 
entire range of discharge and hence for this basin, two separate rating curves, one 
for the low flows and the other for the high flows are developed as shown in Fig. 
3.4(b). 

Furthermore, for some basins a single rating curve was sufficient for the 
entire period for which data were available. For others, for eg., Mujnai, different 
rating curves had to be used for different years because of the variation in the 
stage-discharge relationships over the years. 

3.4.5 Evaporation data 

No evaporation data were available. Hence using Christiansen’s equation, pan 
evaporation values were calculated for Guwahati using monthly normal values of 
the meteorological characteristics obtained from IMD. Using a pan coefficient of 0.7 
and assuming Guwahati data to be representative for the entire region, the daily 
potential evaporation values for the basins were estimated for use in the various 
analyses in the study. 

3.5 Number of Floods in Each Basin 

After detailed scrutiny of records, including comparison of rainfall and 
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Fig. 3. 4b TWIN RATING CURVE-MUJNAI RIVER 
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runoff records, consistent flood events were identified for the study. The number 
of floods in each basin is given in Table 3.3. 


Table 3.3 Number of Floods in Each Basin 


SI .No. 

Name of Basin 

No. of floods 

1 . 

Bat jhora 

11 

2. 

Sankhini 

5 

3. 

Kurt i 

5 

4. 

Mu jna i 

13 

5. 

Talma 

10 

6. 

Khaga 

5 

7. 

Dhars i 

18 

8. 

Mansar i 

5 

9. 

Gaghra 

No suitable flood 

10. 

Ba 1 ason 

5 


3.6 Conclusion 

After the preliminary analysis, it was observed that Gaghra basin did not 
have any floods suitable for analysis because rainfall and runoff data were very 
erratic. The data for river Balason was highly inconsistent. In this case, on certain 
days of heavy rainfall, the stage data did not indicate any floods, while on days 
when the stage was on the rise, there were hardly any noticeable rainfall. Further 
the stage and discharge data also showed wide discrepancies and hence this basin 
was not taken up for the study. Thus after the preliminary analysis, the two basins 
Gaghra and Balason were rejected for the purpose of this study. Further analysis 
was carried out only for the remaining eight basins. 






CHAPTER 4 


REGIONAL I SAT I ON USING UNIT HYDRQGRAPH APPROACH 


4.1 Introduction 

Despite a large number of streams that are gauged, planning problenis quite 
often involve ungauged basins. The estimation of p eak runoff from a catchment is 

needed for comprehensive design of water control structures. The climatic and 
physical characteristics of the catchment are the major factors ^ influence the 

runoff. The climatic factors include precipitation, its duration, areal and temporal 
distribution, and intensity, evapotranspiration, interception etc. The physical 
characteristics which influence runoff are the area of catchment, its length, slope, 
shape, elevation, soil type, drainage system, vegetal cover etc. 

Unit hydrograph (UH) is one of the most popular techniques for the 
computation of runoff from a catchment. A unit hydrogr ap h may by defined as the 
hydrograph resulting from unit depth of surface runoff produced by a storm of 
unit effective rainfall and specified duration. It is considered to represent the 
integrated effect of the various physical factors 0n the routing of the effective 
rainfall through the catchment system. The UH approach assumes that the basin 
response is linear and hence hydrographs produced by different storms of uniform 
intensity and fixed duration but varying volumes of effective rainfall are 
essentially proportional to the volume* 

For gauged catchments, the UH can be derived by analysing the available 
rainfall-runoff records. But when the stream is unga Uged) the UH for the catchment 

has to be derived on the basis of regional relationships using their physical and 
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storm characteristics. This, in turn, necessitates the development of suitable 
regional relationships for estimation of unit hydrographs and hence estimation of 
flood. This is generally done by derivation of UH for gauged catchments for a given 
duration of effective rainfall and then by development of regional relationships 
between UH parameters with the physiographic and storm characteristics of the 
catchments. Thus, for development of regional relationships based on UH approach, 
the various steps involved would be: 

i) Estimate the various UH parameters for a storm in a basin; 

ii) Repeat the step for a number of storms in the basin to explain the 
variation, if any, of the parameters in terms of hydrometeorological 
characteristics; and 

iii) Repeat the above steps for a number of basins in the region so that the 
parameters can be correlated to the physiographic characteristics of the 
basin. 

4.2 Unit Hydrograph Analysis 

The UH approach for determining the hydrograph of design flood assumes 
the catchment as a linear system which transforms the rainfall input into direct 
surface runoff (DSRO) as output. Application of this method involves the prior 
selection of storm and the estimation of loss or rainfall excess and the baseflow 
contribution. 

4.2.1 Direct surface runoff separation 


The total flow due to storm, in practice, can be considered to be consisting of 
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two parts, viz., storm or direct surface runoff (DSRO) and the baseflow. DSRO is 
the runoff resulting at the catchment outlet due to the rainfall excess and the 
baseflow is that contribution to the streamflow which results from release of water 
from surbsurface storage. 

To get DSRO, the baseflow has to be separated from the total runoff. Many 
methods have been suggested (Chow, 1964, Gray ,1970). Using a suitable method, the 
baseflow is separated and the DSRO hydrograph is obtained. In the present study, 
baseflow is considered as constant for a particular storm. 

4.2.2 Effective rainfall 

The part of the rainfall which appears over the surface as runoff and hence 
represents the rainfall contribution to the streamflow is called the excess rainfall. 
The rainfall excess can be determined by either of the following methods; 

(i) The average rainfall over the basin for each time interval can be determined 
by arithmetic average, Theissen polygon method, isohyetal method etc. A 
constant abstraction rate is assumed at all time intervals. The excess of rainfall 
over abstraction rate is the effective rainfall and is otherwise zero. From 
principle of conservation of mass, cumulative effective rainfall during the storm 
should be equal to the total DSRO. The abstraction rate for the storm can be 
calculated by a trial and error procedure and this procedure is referred to as 
the infiltration index method. 

(ii) If instead of taking the average rainfall over the basin, the effective rainfall 
for an abstraction rate is calculated for each of the station and the basin 
average of effective station rainfalls is calculated by any of the methods 
mentioned above, the effective rainfall will be more realistic. It is adopted 
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sometimes by the US Army Corps of Engineers and is referred to as the Corps 
of Engineers procedure. 

The latter method is used in the present study if the number of raingauge 
stations in the basin is more than one. But when there is only one raingauge station 
in the basin, the two methods are identical. 

4.2.3 Use of Nash model 

Various models and methods are available for deriving the unit hydrograph 
for a basin. In order that the UH can be regionalised, it is necessary to represent 
the UH in terms of a limited number of parameters, perhaps in terms of conceptual 
models. A very useful and suitable model for UH is the Nash model. In the model 
proposed by Nash(1957), the watershed is represented by a series of N identical 
linear reservoirs, each having the same storage coefficient K. By routing unit 
impulse effective rainfall excess through N linear reservoirs, a mathematical 
equation for instantaneous unit hydrograph (IUH) for the model can be derived. The 
instantaneous unit hydrograph u(t) for a basin is its UH for an instantaneous unit 
effective rainfall. Nash derived the following equation for the IUH, which is 
mathematically a gamma function 

u(t) - e (_t/K1 (t/K) N1 (4.1) 

where T(N) - (N-l)! and is otherwise a constant which makes the integral Ju(t)dt- 1. 
The two parameters N and K determine the shape of the unit hydrograph. Nash 
generalised his model for noninteger values of N. T(N) can be interpolated from 
tables of complete gamma function when N is not an integer (Abramovitz and 
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Stegun, 1965). A noninteger value of N is considered valid in this study. 

4.2.4 Estimation of UH parameters for the Nash model 

The parameters of the Nash model N and K are to be derived from data for 

defining the UH. Two methods, viz., (i) method of moments and (ii) method of least 

squares are generally used for deriving these parameters. 

(a) Method of Moments: For an initially relaxed (i.e., Q(0) - 0.0) linear system, 

the output is related to the input through the UH by the convolution equation 
t 

Q(t) - Jl(t) U(t-r) dr 0 <; t <oo (4.2) 

0 

where Q(t) represents the output direct runoff hydrograph (DRH); I(t), the input 
effective rainfall hyetograph (ERH) and U(t) is the IUH. Let X(t) be defined over 
0 <; t< oo. The k th moment of a function X(t) about the origin is given by 

oo 

M k . | X(t) t k dt (4.3) 

o 

Since different series may be given in terms of different measurement units, to 
avoid the need for standardisation of inputs and outputs in the same units the 
standard form of the moment equation is given by 

t 

JxCt) Atf 

ul ’ £*)* 

J n-0 

0 



Mktx)- 


for sampled data 


(4.4) 



60 


M k (x) 


]T>«) a - §)* 

n~0 

oo 


for quantised data 


(4-5) 



n-0 

Rainfall is generally given as quantised data and the DSRO as sampled data. By 
taking moments of the convolution equation (4.2) and equating them, it is possible to 
derive various equations relating the moments of Q(t) in terms of those of I(t) and 
U(t). The first and second moments of U(t) for Nash model can be obtained from 
integration as M x (u) - N K and M 2 (u) - N (N+l) K 2 . Method of moments provides a 
simple approach to the calculation of the parameters of the IUH since the moments 
themselves can be expressed in terms of the parameters. For example, for Nash 
model 


M^Q) - MjG) + Mj(u) - M,(I) + NK (4.6) 

MjXQ) - M^I) + M^u) + 2Mj(I) Mj(u) (4.7) 

with M 2 (u) - N (N+l) K 2 . For the two parameters of Nash model, only two moment 
equations are needed to uniquely define the parameters and can be solved to yield 
the values of N and K. This method is known as the method of moments approach to 
the estimation of the model parameters. Further details of estimation of these 
parameters for the sampled data Q and quantised data I are available in reference 
(Chow et al ,1988). 

(b) Method of Least Squares: In this procedure, the parameters N and K are 
estimated by minimising the sum of the squares of the differences between observed 
and computed DSRO hydrographs. This procedure requires an initial estimate of the 
parameters and it estimates the parameters by searching in the parameter space for 
the minimum value of the objective function. 


The method proposed by Nash is easily programmable. However, in the 
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present study, programs developed at the National Institute of Hydrology, Roorkee 
(NIH,1985) have been used with certain modifications. Two computer programs for 
estimating the parameters of the Nash model are available one each for (i) method 
of moments and (ii) method of least squares. Both the programs use the phi-index 
method for rainfall separation. In the present study, the Corps of Engineers 
approach (Subsec. 4.2.2) has been programmed and used whenever there is more than 
one raingauge station in a basin. Furthermore, the microVAX II computer system of 
the Centre for Water Resources Engineering and Management, IIT Kanpur has been 
used for the study and some modifications in the software to suit the microVMS 
FORTRAN compiler were also made. 

4.3. Hydrograph Analysis for a Storm 

The input data for the program for estimating the UH parameters and the 
flood hydrograph by the method of moments for a storm flood event are: 

(i) Catchment Area (Km 2 ); 

(ii) Time interval between successive observations of rainfall and runoff 

(Hours); 

(iii) Number of storms to be analysed; 

(iv) Number of raingauge stations and the corresponding Theissen weights; 

(v) Number of rainfall records for the storm; 

(vi) Rainfall depths at specified time interval at each raingauge station 

for the storm in mm; 

(vii) Number of discharge hydrograph ordinates for the storm; and 

(viii) Ordinates of the discharge hydrograph for the storm (m 3 /sec). 

The discharge hydrograph ordinates for each flood event are obtained from 
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the rating curve that has been prepared for the basin. The rainfall and discharge 
(stage) data are available at an interval of one hour. The baseflow, assumed to be 
constant within a storm equal to the initial value, is separated from the observed 
hydrograph values and the DSRO is obtained. The ERH alongwith the infiltration 
index is obtained after separating the infiltration losses. These operations are 
performed by the various subroutines in the two programs. The results include the 
effective rainfall, the infiltration index and the coefficient of runoff . 

For the second computer program based on the method of least squares, 
initial estimates of the parameters N and K are also required. The program gives the 
values of the parameters N and K and the UH ordinates for the catchment 
corresponding to the minimum values of the objective function estimated from the 
data of the storm precipitation. For a storm-flood event, the Nash model parameters 
N and K are estimated by the method of moments. Using the values of N and K 
obtained by the method of moments as initial values, the UH parameters are 
obtained by the method of least squares. Both the programs compute the ordinates 
of the discharge hydrograph. The observed and the computed hydrographs are also 
plotted to see if the flood event is simulated properly. This procedure is repeated 
for all the flood events identified for the basin. The results are tabulated for 
further analysis. 

To illustrate the results from the analysis, the data and results for a typical 
storm in Dharsi river basin (Fig. 4.1.) are presented. Dharsi river basin has a 
catchment area of 91.40 Km 2 and has two raingauge stations. Table 4.1(a) shows the 
input data for the storm dated 12 th to 14 th September, 1970. The input data consist 
of details of the basin and storm, and in particular, rainfall at each station and 
runoff values. The results of the rainfall separation and DSRO 
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also presented in the Table. They include the baseflow, DSRO, weighted rainfall 
values, infiltration capacity, total rainfall excess and the runoff coefficient as well 
as the separated rainfall excess values. Table 4.1(b) shows the results of estimation 
of Nash model parameters by the method of moments and includes the first and 
second moments of DSRO and ERH respectively, Nash parameters N and K, the IUH 
and the one-hour UH. 

The comparison of observed and computed hydrographs (in Fig.4.2) along with 
a consolidated list of input data and results of separation analysis are presented in 
Table 4.1 (c). It also indicates several parameters of the model and the process 
including observed and computed peaks and times to peak respectively and several 
measures of error. The results indicate that the average absolute error is around 
2%, standard error around 2.6% and the absolute error in peak of about 10.8% which 
is a little high. 

For the analysis using the method of least squares, the input data essentially 
consist of DSRO and the separated rainfall values derived in a similar manner as 
that for the method of moments. The initial parameter estimates are the values 
obtained from the method of moments. The program uses a modified Newton method 
for minimising the sum of squared errors with a relative tolerance limit of 0.001. 
The results of analysis are presented in Table 4.2. which is self explanatory. A 
comparison of the errors indicate that the average standard error and average 
absolute error are of the order of 1%. The average percentage absolute error is 
5.26% and the percentage absolute error in peak is 0.07% and the observed and the 
computed time to peak are the same. Hence the method of least squares seems to 
predict the hydrograph better than the method of moments (Fig.4.2). 
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TABLE 4.1 (a) INPUT DATA AND RESULTS OF RAINFALL SEPARATION 
UNIT HYDROGRAPH ANALYSIS - NASH MODEL 


CATCHMENT DHARSI 

DATE 12-14, SEPT, 1970 

CATCHMENT AREA (Ktn 2 ) 91.40 
NO.OF RAINGAUGES 2 

THEISSEN WEIGHTS OF EACH RAINGAUGE STATION 

0.7700 0.2300 

NO OF RAINFALL VALUES= 21 

RAINFALL AT EACH STATION (MM) 

STATION NO: 1 


3.05 

29.97 

0.00 

0.00 

7.72 

5.08 4.57 

0.00 0.00 

2.54 

1.27 

4.57 

1.02 

10.67 

0.00 

25.91 

0.51 

50.80 

0.51 

27.94 

0.00 

STATION NO: 2 
6.10 0.00 
6.10 22.60 

2.03 

25.40 5.59 

0.00 0.00 

0.00 

0.00 

7.37 

0.50 

0.13 

0.00 

31.00 

1.27 

24.90 

3.05 

6.35 

2.29 

NO. OF 

RUNOFF 

VALUES - 37 








OBSERVED DISCHARGE HYDROGRAPH (M 3 /SEC) 


11.75 

11.75 

11.75 

11.75 

11.75 

11.75 

11.75 

13.71 

15.23 

16.99 

18.80 

20.76 

24.30 

27.41 

29.05 

30.75 

30.75 

30.75 

31.63 

31.63 

30.75 

29.05 

27.41 

25.83 

24.30 

22.85 

22.85 

22.14 

22.14 

22.14 

21.44 

20.76 

19.45 

18.80 

18.80 

18.18 

18.18 





WEIGHTED RAINFALL VALUES (MM) 


3.75 

0.00 

9.75 

4.80 

1.96 

5.21 

10.09 

27.08 

44.84 

22.97 

23.48 

11.14 

0.00 

0.00 

0.98 

0.90 

0.00 

0.68 

1.09 

0.53 

0.47 










NO. OF 

DSRO ORDINATES 

- 29 











BASEFLOW (M 3 /SEC) 




11.75 

11.75 

11.75 

11.75 

11.75 

11.75 

11.75 

11.75 

11.75 

11.75 

11.75 

11.75 

11.75 

11.75 

11.75 

11.75 

11.75 

11.75 

11.75 

11.75 

11.75 

11.75 

11.75 

11.75 

11.75 

11.75 

11.75 

11.75 

11.75 
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Table 4.1 (a) Co ntd.. 

DIRECT SURFACE RUNOFF (M 3 /SEC) 


3.48 

5.24 

7.05 

9.01 

12.55 

15.66 

17.30 

19.00 

19.00 

19.00 

19.88 

19.88 

19.00 

17.30 

15.66 

14.07 

12.55 

11.10 

11.10 

10.39 

10.39 

10.39 

9.69 

9.01 

7.70 

7.05 

7.05 

6.43 

6.43 



TOTAL RAINFALL (MM) 170.74 
INFILTRATION CAPACITY (MM/HR) 32.78 
TOTAL RAINFALL EXCESS 13.88 
RUNOFF COEFFICIENT 0.0813 

SEPARATED RAINFALL VALUES (MM) 

13.88 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 0.00 
0.00 0.00 0.00 


TABLE 4.1 (b) RESULTS OF ESTIMATION OF NASH MODEL PARAMETERS BY 

METHOD OF MOMENTS 


FIRST MOMENT OF DSRO 13.030 

SECOND MOMENT OF DSRO 218.720 

FIRST MOMENT OF ERH 0.500 

SECOND MOMENT OF ERH 0.250 

VALUE OF N 3.200 

VALUE OF K 3.915 


I.U.H. ORDINATES 


0.103 

0.365 

0.691 

1.007 

1.275 

1.474 

1.603 

1.666 

1.672 

1.633 

1.560 

1.463 

1.351 

1.232 

1.111 

0.991 

0.878 

0.771 

0.672 

0.583 

0.503 

0.431 

0.368 

0.313 

0.266 

0.224 

0.189 

0.158 

0.132 



SUM OF IUH 0.975 

IUH PEAK 1.672 

IUH TIME TO PEAK 9 HOURS 


AREA OF UH 0.972 


UNIT HYDROGRAPH ORDINATES (M 3 /SEC) 


0.296 

0.324 

0.515 

0.850 

1.146 

1.380 

1.544 

1.639 

1.673 

1.655 

1.598 

1.513 

1.408 

1.292 

1.171 

1.051 

0.934 

0.823 

0.721 

0.627 

0.542 

0.466 

0.399 

0.340 

0.289 

0.244 

0.206 

0.173 

0.145 
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TABLE 4.1(c) COMPARISON OF OBSERVED AND COMPUTED HYDROGRAPHS 


ORD.NO 

RFSTN1 

RFSTN2 

WTRFALL EFFRF BASE 
FLOW 

DSRO 

OBSD. 

DISCH. 

COMPUTED 

DISCHARGE 


MM 

MM 

NM 

HH 





1 

3.05 

6.10 

3.75 

0.00 

11.75 

0.00 

11.75 

11.75 

2 

0.00 

0.00 

0.00 

0.00 

11.75 

0.00 

11.75 

11.75 

3 

5.08 

25.40 

9.75 

0.00 

11.75 

0.00 

11.75 

11.75 

4 

4.57 

5.59 

4.80 

0.00 

11.75 

0.00 

11.75 

11.75 

5 

2.54 

0.00 

1.96 

0.00 

11.75 

0.00 

11.75 

11.75 

6 

4.57 

7.37 

5.21 

0.00 

11.75 

0.00 

11.75 

11.75 

7 

10.67 

8.13 

10.09 

0.00 

11.75 

0.00 

11.75 

11.75 

8 

25.91 

31.00 

27.08 

0.00 

11.75 

0.00 

13.71 

11.75 

9 

50.80 

24.90 

44.84 

13.88 

11.75 

0.00 

15.23 

11.75 

10 

27.94 

6.35 

22.97 

0.00 

11.75 

4.12 

16.99 

15.87 

11 

29.97 

6.10 

24.48 

0.00 

11.75 

5.33 

18.80 

17.08 

12 

7.72 

22.60 

11.14 

0.00 

11.75 

7.16 

20.76 

18.91 

13 

0.00 

0.00 

0.00 

0.00 

11.75 

11.84 

24.30 

23.59 

14 

0.00 

0.00 

0.00 

0.00 

11.75 

15.95 

27.41 

27.70 

15 

1.27 

0.00 

0.98 

0.00 

11.75 

19.21 

29.05 

30.96 

16 

1.02 

0.50 

0.90 

0.00 

11.75 

21.50 

30.75 

33.25 

17 

0.00 

0.00 

0.00 

0.00 

11.75 

22.82 

30.75 

34.57 

18 

0.51 

1.27 

0.68 

0.00 

11.75 

23.28 

30.75 

35.04 

19 

0.51 

3.05 

1.09 

0.00 

11.75 

23.04 

31.63 

34.79 

20 

0.00 

2.29 

0.53 

0.00 

11.75 

22.25 

31.63 

34.00 

21 

0.00 

2.03 

0.47 

0.00 

11.75 

21.06 

30.75 

32.81 

22 

0.00 

0.00 

0.00 

0.00 

11.75 

19.60 

29.05 

31.35 

23 

0.00 

0.00 

0.00 

0.00 

11.75 

17.98 

27.41 

29.73 

24 

0.00 

0.00 

0.00 

0.00 

11.75 

16.30 

25.83 

28.05 

25 

0.00 

0.00 

0.00 

0.00 

11.75 

14.62 

24.30 

26.38 

26 

0.00 

0.00 

0.00 

0.00 

11.75 

13.00 

22.85 

24.75 

27 

0.00 

0.00 

0.00 

0.00 

11.75 

11.46 

22.85 

23.21 

28 

0.00 

0.00 

0.00 

0.00 

11.75 

10.03 

22.14 

21.78 

29 

0.00 

0.00 

0.00 

0.00 

11.75 

8.73 

22.14 

20.48 

30 

0.00 

0.00 

0.00 

0.00 

11.75 

7255 

22.14 

19.30 

31 

0.00 

0.00 

0.00 

0.00 

11.75 

6.49 

21.44 

18.24 

32 

0.00 

0.00 

0.00 

0.00 

11.75 

5.56 

20.76 

17.31 

33 

0.00 

0.00 

0.00 

0.00 

11.75 

4.74 

19.45 

16.49 

34 

0.00 

0.00 

0.00 

0.00 

11.75 

4.02 

18.80 

15.77 

35 

0.00 

0.00 

0.00 

0.00 

11.75 

3.40 

18.80 

15.15 

36 

0.00 

0.00 

0.00 

0.00 

11.75 

2.87 

18.18 

14.62 

37 

0.00 

0.00 

0.00 

0.00 

11.75 

2.41 

18.18 

14.16 


UH PEAK (M 3 /SEC) - 1.67 

UH TIME TO PEAK (HRS) - 9.00 

EFFICIENCY OF THE MODEL (%) - 84.10 

OBSERVED PEAK (M 3 /SEC) - 31.63 

OBSERVED TIME TO PEAK (HRS) - 19.00 

COMPUTED PEAK (M 3 /SEO - 35.04 

COMPUTED TIME TO PEAK (HOURS) - 18.00 

AVERAGE STANDARD ERROR - 2.63 

AVERAGE ABSOLUTE ERROR - 2.08 

AV. PERCENTAGE ABSOLUTE ERROR - 9.47 

PERCENTAGE ABSOLUTE ERROR IN PEAK - 10.77 

PERCENTAGE ABSOLUTE ERROR IN TIME TO PEAK- 5.26 




TABLE 4.2 RESULTS OF ANALYSIS BY NASH MODEL USING METHOD OF LEAST SQUARI 


FINAL SOLUTION FOR THE PARAMETERS 


N 

- 

3.3339 








K 

- 

4.0663 








AREA 

OF UH 

0.98417 











UNIT 

HYDROGRAPH ORDINATES 




0.269 

0.465 

0.566 

0.664 

0.936 

1.170 

1.351 

1.476 

1.546 

1.567 

1.549 

1.498 

1.423 

1.332 

1.231 

1.125 

1.018 

0.913 

0.813 

0.719 

0.632 

0.552 

0.480 

0.416 

0.358 

0.308 

0.263 

0.224 

0.191 

0.162 

0.137 

0.115 

0.097 

0.081 











ESTIMATED DIRECT SURFACE RUNOFF (M 3 /SEC) 


0.111 

0.152 

0.213 

0.386 

0.544 

2.432 

3.638 

4.919 

7.053 

9.741 

13.761 

16.698 

18.943 

20.455 

21.263 

21.445 

21.104 

20.350 

19.292 

18.026 

16.634 

15.183 

13.727 

12.306 

10.950 

9.678 

8.502 

7.427 

6.457 

5.587 

4.814 

4.132 

3.534 

3.013 








COMPARISON OF OBSERVED AND COMPUTED HYDROGRAPHS 


ORD.NO. 

OBSERVED DISCHARGE 
mM-TEc. 

COMPUTED DISCHARGE 

Cl) 

(2) 

(3) 

1 

11.75 

11.75 

2 

11.75 

11.75 

3 

11.75 

11.75 

4 

11.75 

11.75 

5 

11.75 

11.96 

6 

11.75 

12.06 

7 

11.75 

12.14 

8 

13.71 

12.30 

9 

15.23 

15.39 

10 

16.99 

16.18 

11 

18.80 

17.67 

12 

20.76 

18.80 

13 

24.30 

21.49 

14 

27.41 

25.51 

15 

29.05 

28.45 

16 

30.75 

28.99 

17 

30.75 

30.30 

18 

30.75 

31.26 

19 

31.63 

31.58 

20 

31.63 

31.14 

21 

30.75 

30.42 

22 

29.05 

29.48 

23 

27.41 

28.40 
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TABLE 4.2 Contd... 


(1) 

(2) 

(3) 

24 

25.83 

27.24 

25 

24.30 

26.13 

26 

22.85 

25.48 

27 

22.85 

24.06 

28 

22.14 

22.70 

29 

22.14 

21.43 

30 

22.14 

20.25 

31 

21.44 

19.18 

32 

20.76 

18.21 

33 

19.45 

17.34 

34 

18.80 

16.57 

35 

18.80 

15.88 

36 

18.18 

15.29 

37 

18.18 

14.76 


UH PEAK (M 3 /SEC) - 1.567 

UH TIME TO PEAK (HRS) - 10.00 

EFFICIENCY OF THE MODEL - 94.91 

OBSERVED PEAK (M 3 /SEC) - 31.63 

OBSERVED TIME TO PEAK (HRS) - 19.00 

COMPUTED PEAK (M 3 /SEC) - 31.58 

COMPUTED TIME TO PEAK (HRS) - 19.00 

AVERAGE STANDAR ERROR - 1.542 

AVERAGE ABSOLUTE ERROR =1.006 

AV. PERCENTAGE ABSOLUTE ERROR = 5.260 
PERCENTAGE ABSOLUTE ERROR IN PEAK » 0.070 
PERCENTAGE ABS. ERROR IN TIME TO PEAK 0.000 




<- LV I mi\oi urAO 1 1 \i — IVI r N 

OF OBSERVED AND COMPUTED HYDROGRAPHS 
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4.4 Rainfall Runoff Analysis for a Basin 

The data for Dharsi basin were available for the years 1968 to 1970 and 18 
flood events were identified for analysis. Table 4.3 shows the consolidated results 
for all the 18 floods in the basin. The first two columns in the Table give the the 
flood number and the date of occurrence of the event. The next two 
columns indicate respectively the observed peak discharge value for the 
flood (m 3 /sec) and its time of occurrence. Column 5 shows the total rainfall (mm) for 
the flood. The next four columns give the results of the separation analysis which 
include the effective rainfall (mm), baseflow (m 3 /sec), infiltration index (mm/hour) 
and the runoff coefficient. 

From the Table, it is observed that the baseflow in this basin varies from 
11.75 m 3 /sec to 28.20 m 3 /sec. There is a wide variation in the values of the 
infiltration index and the runoff coefficient. A very high runoff coefficient is 
associated with the highest observed flood peak and it has the lowest infiltration 
index (2.414 mm/hour) and relatively low rainfall. A very low runoff coefficient is 
associated with one of the lowest flood peaks and it has the highest infiltration 
index (33.29 mm/hour) and the highest rainfall (218.4 mm/hour). It is found that the 
behaviour of the infiltration index and the runoff coefficient with respect to the 
peak discharge is erratic. This indicates that the floods in this basin are caused by 
highly variable precipitation and highly variable abstraction. 

Columns 10 to 14 represent the results of the Nash model using the method 
of moments. It is observed that in most cases, the computed peak values are within 
±10% and generally around 7% or less, except for flood number 11 where there is 


a large variation. 



TABLE 4.3 RESULTS OF ANALYSIS USING NASH MODEL 


DHARSI RIYER BASIN CATCHMENT AREA 91.40 SQ.KM 


DATE 

QPOBD. 

M**3/SEC 

OB. TIME 
HRS 

RFALL 

MM 

BFFRF 

MM 

BASEFLO 

M**3/SEC 

PHIINDEX R-0 COEFFT 
MM/HOUR 

2 

3 

4 

5 

6 

7 

8 

9 

260668 

39.220 

13.000 

79.250 

12.884 

18.180 

12.556 

0.162 

300668 

106.760 

11.000 

88.590 

64.009 

20.110 

2.414 

0.757 

12130768 

44.540 

22.000 

188.740 

20.885 

11.750 

31.720 

0.111 

140768 

83.170 

10.000 

66.030 

32.497 

28.200 

7.117 

0.492 

170768 

32.510 

11.000 

55.370 

6.997 

21.440 

18.913 

0.126 

20210768 

57.850 

14.000 

119.120 

30.892 

20.110 

14.882 

0.259 

200968 

49.130 

10.000 

63.000 

23.279 

18.180 

7.227 

0.370 

290669 

52.760 

13.000 

123.410 

18.186 

21.096 

18.818 

0.147 

15160869 

47.940 

12.000 

218.440 

19.775 

14.730 

33.292 

0.091 

20969 

75.100 

11.000 

74.920 

32.642 

25.060 

4.460 

0.436 

120969 

45.680 

14.000 

71.880 

24.502 

15.260 

9.432 

0.341 

340670 

69.010 

13.000 

179.830 

24.290 

20,110 

27.460 

0.135 

450670 

102.880 

13.000 

148.160 

47.641 

19.454 

19.540 

0.322 

150770 

72.010 

11.000 

124.460 

25.335 

25.061 

23.728 

0.204 

230770 

91.780 

14.000 

106.940 

48.098 

20.110 

6.087 

0.450 

50870 

37.240 

9.000 

125.740 

12.066 

20.110 

16.416 

0.096 

80870 

44.660 

6.000 

91.800 

13.780 

25.060 

14.457 

0.150 

12140970 

31.630 

19.000 

170.740 

13.880 

11.750 

32.778 

0.081 


IOM.TIME 

METHOD OF MOMENTS 

N K 

METHOD OF LEAST SQUARES 

MK QPCOMP. COMP. TIME N K MI 

N* 

K* 

ESTIMATED 

MK* QPEAK 

T1 

HOURS 


M**3/SEC HOURS 



M**3/SEC 

HOI 


11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

2 

13.000 

3.774 

3.724 

14.054 

37.900 

14.000 

4.000 

3.690 

14.760 

3.327 

3.996 

13.297 

39.210 

13.06 

10.000 

1.490 

7.953 

11.850 

108.730 

10.000 

1.385 

9.265 

12.832 

1.301 

7.779 

10.121 

106.770 

11.00 

25.000 

3.296 

4.534 

14.944 

39.930 

25.000 

2.993 

5.080 

15.204 

3.168 

4.294 

13.603 

44.540 

23.00 

11.000 

2.233 

4.449 

9.935 

82.550 

10.000 

1.840 

5.910 

10.874 

2.009 

6.458 

12.973 

83.180 

10. 00 

10.000 

2.992 

3.564 

10.663 

33.450 

11.000 

3.262 

3.776 

12.317 

3.529 

3.621 

12.776 

32.510 

11.00 

14.000 

2.691 

5.541 

14.911 

60.190 

14.000 

2.691 

5.541 

14.911 

2.768 

5.040 

13.952 

57.840 

14.00 

13.000 

2.593 

5.325 

13.808 

50.010 

13.000 

2.593 

5.325 

13.808 

3.030 

4,551 

13.791 

49.120 

1.00 

17.000 

3.883 

3.444 

13.373 

49.890 

14.000 

3.383 

3.938 

13.322 

2.921 

4.755 

13.889 

52.770 

14.00 

11.000 

2.462 

4.308 

10.606 

48.800 

11.000 

2.302 

4.734 

10.898 

3.066 

4.485 

13.749 

50.000 

13.00 

10.000 

1.706 

5.824 

9.936 

76.190 

11.000 

1.650 

6.628 

10.936 

2.251 

6.006 

13.519 

75.800 

11.00 

16.000 

3.023 

4.136 

12.503 

49.750 

17.000 

2.908 

4.925 

14.322 

3.134 

4.358 

13.656 

46.000 

5.00 

12.000 

2.583 

4.269 

11.027 

61.440 

13.000 

2.669 

4.294 

11.461 

2.434 

5.665 

13.786 

68.300 

3J6 

11.000 

1.452 

6.498 

9.435 

101.620 

11.000 

1.605 

6.416 

10.298 

1.418 

7.561 

10.719 

102.580 

2.06 

9.000 

2.079 

4.386 

9.118 

71.520 

10.000 

2.287 

4.528 

10.356 

2.344 

5.833 

13.670 

72.100 

10.00 

12.000 

1.664 

7.131 

11.866 

82.980 

13.000 

1.736 

7.720 

13.402 

1.751 

6.940 

12.149 

90.870 

13.06 

12.000 

3.411 

3.791 

12.931 

39.640 

12.000 

3.411 

3.791 

12.931 

3.387 

3.885 

13.159 

38.130 

11.06 

14.000 

4.294 

3.028 

13.002 

47.030 

14.000 

4.306 

3.062 

13.185 

3.164 

4.301 

13.609 

44.560 

10.06 

18.000 

3.200 

3.915 

12.528 

31.650 

19.000 

3.334 

4.066 

13.556 

3.555 

3.571 

12.696 

31.690 

19.0! 
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Columns 15 to 18 show the results of the Nash model by the method of least 
squares. The computed peaks are generally closer to the observed data than those 
obtained by the method of moments and the error in the peak discharge values are 
of the order of 3%. However for three floods, viz., 3, 12 and 15, the error in the 
prediction of peak discharge is larger and perhaps a weighted least squares may 
lead to a better result. The computed time to peak are often nearly the same as 
that for observed data except in a few cases where the difference is around two to 
three hours. However, the deviation in the case of method of least squares is 
smaller compared to that in the case of method of moments. Hence, for considering 
the variation of UH parameters and basin characteristics, it is proposed to use the 
results of the method of least squares. 

4.4.1 Parameter variation in a basin 

The values of N, K and NK are shown in columns 17 to 19. The values of N, 
K and NK are plotted against the observed peaks in Fig. 4.3. NK essentially 
represents the time lag for the storm-flood event and varies from about 10.9 to 15 
hours, with a decreasing trend for NK as Q p increases. This is logical as larger 
floods will run off quicker. The figure also indicates clearly that the value of N 
decreases and K increases respectively with the flood peak. A regression line is 
fitted to these variation by the method of least squares, viz., 


N - 4.504 - 0.030 Qpj and 

(4.8) 

K - 1.800 + 0.056 Qp 

(4.9) 


The variation of NK as a function of Q p is not obvious from the figure. Data errors 
also have been noted for some storms. The derived equation for the total 
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time lag in terms of the fitted equation for N and K is given by 

NK - 8.1072 - 0.1982 Qp - 0.00168 Qp 2 (4.10) 

This is also shown in the figure and it seems to fit the observations satisfactorily. 
The decrease in NK values for smaller Q p s may be because of a storm occurring 
near the gauging station producing a larger flood peak (Ramaseshan,1964). The 
parameter estimates (Eqns. 4.8 to 4.10) are shown in columns 20 to 22 respectively. 
The estimated Q p and time to peak for those values are shown in columns 23 and 24 
respectively; and for the storm-flood event of Sec. 4.3 in Fig. 4.2. The fit is 
considered satisfactory. 

Other parameters in the process include coefficient of runoff, infiltration 
index and baseflow. Attempts were made to correlate these parameters with the 
storm and flood characteristics. While some general trends were obvious, it is not 
possible to fit definitive relationships for their variation in a basin especially since 
the variability about the fitted relationships was very large. 

4.5 Comparative Results of Analysis for Different Basins 

For each of the basins, rainfall-runoff analysis for all the storms in the 
basin were carried out. The results are tabulated in Tables 4.4 through 4.10. The 
plots of observed maximum discharge for each of the floods against the model 
parameters are also shown in Figs. 4.4 through 4.10. 

4.5.1 Routing parameters (Nash Model) 

Batjhora basin has a catchment area of 7.33 Km 2 and the results of the 
analysis of the basin are presented in Table 4.4 and Fig. 4.4. They indicate that the 



TABLE 4.4 RESULTS OF ANALYSIS USING NASH KODEL 


BATJHORA RIVER BASIN CATCHMENT AREA 7.33 SQ.KM 



DATE 

QPOBD . 

OB. TINE 

RFALL 

EFFRF 

BASEFLO 

PHIINDEXR-0 COEF 



MM3/SEC 

HOURS 

NN 

NN 

N**3/SEC 

NN/HOUR 


1 

2 

3 

4 

5 

6 

7 

8 

9 

1 

180764 

6.800 

6.000 

29.718 

9.042 

2.270 

3.342 

0.304 

2 

190764 

14.100 

3.000 

88.138 

33.609 

2.800 

15.708 

0.381 

3 

200764 

12.430 

3.000 

58.928 

21.565 

3.650 

4.449 

0.366 

4 

260764 

4.590 

5.000 

36.068 

11.876 

1.950 

4.842 

0.326 

5 

20964 

9.400 

4.000 

59.944 

13.172 

2.800 

10.677 

0.220 

6 

30964 

8.040 

9.000 

79.248 

23.578 

2.800 

6.670 

0.296 

7 

110964 

33.730 

5.000 

75.946 

71.633 

3.650 

0.479 

0.954 

8 

120964 

12.430 

7.000 

87.630 

36.948 

3.650 

13.337 

0.422 

9 

80865 

28.660 

6.000 

84.836 

54.476 

5.650 

2.915 

0.642 

10 

80965 

5.640 

12.000 

92.710 

8.963 

1.610 

13.884 

0.116 

11 

180965 

8.040 

3.000 

53.594 

11.020 

2.270 

10.000 

0.206 




NETHOD OF NONENTS 




NETHOD OF LEAST SQUARES 


ESTIMATED 


• 

CON. TINE 

N 

K 

NK 

QPCONP . 

CONP.TINE 

N 

1 

NK 

N» 

IT 

NET 

BC 

HOURS 




N**3/SEC 

HOURS 







10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

(0 

7.000 

1.549 

2.184 

3.383 

6.080 

6.000 

1.408 

2.184 

3.075 

1.705 

2.776 

4.734 

70 

4.000 

1.783 

1.720 

3.067 

15.150 

4.000 

1.401 

2.707 

3.793 

1.508 

2.747 

4.143 

70 

3.000 

1.346 

1.885 

2.537 

13.170 

4.000 

2.685 

0.867 

2.328 

1.553 

2.753 

4.277 

10 

8.000 

2.107 

2.715 

5.721 

4.170 

7.000 

1.461 

4.740 

6.925 

1.765 

2.785 

4.915 

10 

4.000 

2.029 

1.263 

2.563 

8.420 

4.000 

1.480 

1.919 

2.840 

1.635 

2.765 

4.522 

>0 

8.000 

1.486 

2.827 

4.201 

8.020 

10.000 

1.898 

3.324 

6.309 

1.672 

2.771 

4.633 

[0 

5.000 

0.967 

2.799 

2.707 

33.730 

5.000 

0.665 

3.982 

2.648 

0.978 

2.668 

2.610 

70 

7.000 

2.070 

2.428 

5.026 

11.570 

7.000 

1.701 

3.732 

6.348 

1.553 

2.753 

4.277 

70 

7.000 

1.197 

1.572 

1.882 

25.760 

7.000 

1.222 

1.645 

2.010 

1.115 

2.688 

2.998 

10 

11.000 

1.992 

2.866 

5.709 

4.280 

11.000 

1.896 

3.288 

6.234 

1.737 

2.780 

4.829 

10 

3.000 

1.370 

1.175 

1.610 

8.230 

3.000 

1.073 

1.840 

1.974 

1.672 

2.771 

4.633 
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lag time for large floods is smaller. However for small floods, the data are erratic. 
The regression equation fitted for N, K and NK respectively are considered to be 
reasonable and are accepted in this study. 

Sankhini basin has a catchment area of 7.76 Km 2 and the results for the basin 
are given Table 4.5 and Fig.4.5. A comparison of results for storms 2 and 5 indicate 
large inconsistency of data in that a large flood takes a longer time to travel. The 
fitted relationships indicate an increase of K with Q p and a slight increase of the 
lag time with Q p . This is not logical. Hence a uniform value of K is adopted for the 
basin. 

Kurti basin has a catchment area of 21.91 Km 2 and the results for this basin 
are shown in Table 4.6 and Fig.4.6. The results are quite contradictory in that two 
storms with different effective rainfall each have the same Q p . Storm 3 has the 
highest flood peak but also has the largest time lag. Detailed analysis of storm 
precipitation data indicates two observations of the order of 78mm and 55mm 
which seem to dominate the process and contribute to effective rainfall while 
precipitation of the order of 40 mm/hour does not produce any runoff. The high 
value of phi-index of 41.82 mm/hour for the storm is also unduly large. The results 
indicate that there are data errors which cause the problem. Similarly the results 
of the storm indicating the lowest lag time for one of the lowest observed Q p does 
not seem to be reliable. Since the data for the basin seem to be inconsistent and 
erratic, it is proposed not to consider this basin for further analysis. 

Mujnai basin has a catchment area of 38.50 Km 2 and has 13 storm flood 
events. The results of the analysis for the basin are shown in Table 4.7 and Fig .4 .7 
respectively. Even though there are variations in the results, they are generally 



Fig. 4. 4 BATJHORA BASIN- QPEAK VS N,K,NK 



QPEAK M**3/SEC 

□ N + K o NK v NKFTTTED 


Fig. 4. 5 SANKHINI BASIN- QPEAK VS N,K,NK 





TABLE 4.5 RESULTS OF ANALYSIS USING NASA NOBEL 


SANKHINI RIVER BASIN CATCHKBNT AREA 7.76 SQ.KM 



DATE 

QPOBD . 

OB. TINE 

RFALL 

EFFRF 

BASEFLO PHIINDEX 

R-0 COEF 



N**3/SBC 

HOURS 

NN 

NN 

M**3/SEC NN/HOUR 


1 

2 

3 

4 

5 

6 

7 

8 

9 

1 

290764 

16.200 

6.000 

125.730 

42.760 

6.630 

19.528 

0.340 

2 

240765 

10.150 

14.000 

112.776 

44.884 

2.690 

8.140 

0.399 

3 

120865 

11.260 

12.000 

159.766 

22.078 

6.610 

26.680 

0.138 

4 

140865 

11.840 

15.000 

98.552 

23.163 

7.560 

5.216 

0.235 

5 

18190965 

4.480 

6.500 

67.310 

16.511 

0.730 

10.795 

0.245 




NETHOD OF 

NONENTS 




NETHOD OF 

LEAST SQUARES 


ESTINATED 


1 

IEC 

CON. TINE 
HOURS 

N 

E 

NK 

QPCONP. 

N**3/SEC 

CONP.TINE 

HOURS 

N 

1 

NK 

N* 

1* 

Ni* 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

30 

10.000 

3.874 

2.296 

8.895 

14.900 

9.000 

2.994 

2.956 

8.850 

2.694 

3.321 

8.948 

60 

11.000 

3.020 

2.850 

8.607 

10.020 

13.000 

3.278 

3.251 

10.657 

2.658 

3.164 

8.409 

90 

10.000 

2.330 

2.616 

6.095 

10.850 

10.000 

1.850 

4.390 

8.122 

2.665 

3.193 

8.507 

40 

15.000 

2.564 

2.334 

5.984 

12.000 

16.000 

2.470 

2.623 

6.479 

2.668 

3.208 

8.559 

10 

7.000 

3.687 

1.780 

6.563 

4.250 

7.000 

2.723 

2.686 

7.314 

2.624 

3.016 

7.915 



TABLE 4.6 RESULTS OP ANALYSIS USING NASH MODEL 


KURTI RIVER BASIN CATCHNENT AREA 21.91 SQ.KM. 



DATE 

QPOBD. 

OB. TINE 

RPALL 

EPPRP 

BASEFLO 

PHIINDEI 

R-0 CGEP 



N**3/SEC 

HOURS 

NN 

NN 

N**3/SEC 

NN/HOUR 


1 

2 

3 

4 

5 

6 

7 

8 

9 

1 

3010766 

30.230 

15.000 

337.058 

60.327 

4.930 

18.046 

0.179 

2 

100766 

30.230 

10.000 

218.948 

44.814 

5.870 

19.608 

0.205 

3 

20210766 

63.550 

7.000 

312.420 

45.916 

21.300 

41.812 

0.147 

4 

70967 

41.810 

10.500 

114.046 

62.483 

4.930 

6.459 

0.548 

5 

13140768 

41.810 

14.000 

113.538 

61.872 

8.770 

4.446 

0.545 




NETHOD OP 

NONENTS 




NETHOD OP 

LEAST SQUARES 


ESTINATED 


QPCONP. 

CON. TINE 

N 

K 

NK 

QPCONP. 

CONP.TINE 

N 

K 

NR 

N* 

R* 

NR* 

N**3/SEC 

HOURS 




N**3/SEC 

HOURS 







10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

37.000 

15.000 

1.688 

3.782 

6.384 

32.300 

15.000 

1.306 

6.107 

7.976 

1.227 

5.396 

6.621 

34.750 

10.000 

1.405 

3.446 

4.842 

34.500 

10.000 

1.163 

4.515 

5.251 

1.227 

5.396 

6.621 

59.380 

9.000 

3.394 

2.583 

8.767 

58.500 

8.000 

2.840 

3.072 

8.724 

2.826 

3.163 

8.940 

51.980 

11.000 

1.310 

3.735 

4.893 

44.400 

11.000 

0.946 

7.222 

6.832 

1.783 

4.620 

8.236 

44.-240 

11.000 

2.545 

2.067 

5.261 

42.440 

11.000 

2.631 

2.270 

5.973 

1.783 

4.620 

8.236 



Fig. 4.6 KURTI BASIN - QPEAK VS N,K,NK 
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T A B L B 4.7 RESULTS OF ANALYSIS USING NASH MODEL 


MUJNAI RIVBR BASIN CATCHMENT AREA 38.50 SQ.KM 



DATE 

QPOBD. 

OB. TIME 

RFALL 

EFFRF 

BASEFLO 

PHIINDEX R-0 COEFFT 



M**3/SEC 

HOURS 

MM 

MM 

K**3/SEC 

MM/HOUR 


1 

2 

3 

4 

5 

6 

7 

8 

9 

1 

260664 

52.440 

9.000 

81.681 

23.396 

16.905 

7.114 

0.286 

2 

140764 

52.440 

9.000 

135.037 

30,318 

15.433 

12.928 

0.225 

3 

200764 

38.060 

9.000 

122.480 

20.767 

13.960 

12.312 

0.170 

4 

280764 

69.550 

12.000 

162.027 

45.578 

17.132 

11.579 

0.281 

5 

110964 

54.280 

9.000 

96.036 

24,096 

15.433 

6.712 

0.251 

6 

120964 

45.110 

8.000 

72.353 

13.687 

18.406 

8.078 

0.189 

7 

140964 

104.150 

9.000 

125.401 

51.980 

19.935 

8.334 

0.415 

8 

60765 

13.900 

11.000 

85.501 

2.380 

10.392 

10.095 

0.028 

9 

200765 

25.170 

7.000 

46.688 

6.085 

14.640 

7.248 

0.130 

10 

190865 

32.000 

6.000 

38.797 

5.126 

24.749 

7.856 

0.132 

11 

130965 

30.270 

6.000 

83.151 

15.945 

20.615 

7.177 

0.192 

12 

140965 

33.670 

11.000 

85.993 

13.668 

22.711 

7.277 

0.159 

13 

150965 

54.820 

12.000 

88.685 

27.379 

24.749 

6.341 

0.309 


METHOD OF MOMENTS METHOD OF LEAST SQUARES 

QPCOMP. COM. TIME N I NK QPCOMP. COMP. TIME N l NK It 
M**3/SEC HOURS M**3/SEC HOURS 

10 11 12 13 14 15 16 17 18 19 20 


51.150 

9.000 

1.826 

2.220 

4.054 

49.090 

52.580 

10.000 

2.515 

2.137 

5.375 

51.740 

39.910 

10.000 

2.630 

2.105 

5.536 

38.240 

77.140 

11.000 

3.279 

1.874 

5.489 

72.450 

52.760 

10.000 

2.114 

1.879 

3.972 

53.580 

47.950 

8.000 

1.418 

1.930 

2.737 

45.380 

89.790 

9.000 

1.675 

2.166 

3.628 

92.150 

16.170 

12.000 

1.883 

1.678 

3.160 

14.090 

27.650 

4.000 

1.391 

1.740 

2.420 

25.300 

32.970 

5.000 

2.382 

2.065 

4.919 

31.040 

33.950 

14.000 

1.552 

3.613 

5.607 

31.770 

36.340 

13.000 

2.803 

1.406 

3.941 

34.530 

51.390 

13.000 

1.509 

3.284 

4.956 

51.490 


9.000 1.681 2.629 4.419 1.825 U91 , m 

9.000 2.312 2.351 5.436 1.825 2.691 im 

9.000 2.217 2.620 5.809 1.709 2.993 5 « 

11.000 2.986 2.005 5.987 1.961 2 331 ' 

9.000 1.592 2.370 3.773 1.839 2.652 « o 

8.000 1.325 2.399 3.179 1.766 2.845 5 fl ! 

9.000 1.873 1.863 3.489 2,238 1.605 5 ? 

12.000 1.332 3.731 4.970 1.516 3.500 ini 

5.000 1.344 2.391 3.214 1.606 3.263 5 ?< 

5.000 1.768 3.498 6.184 1.661 3.120 iw 

15.000 1.772 4.272 7.570 1.847 3.156 51 

13.000 1.745 3.062 5.343 1.674 3.085 j’iJ! 

13.000 1.382 3.624 5.008 1.844 2.641 4 g6 
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consistent and the fitted regression equations appear to be satisfactory. 

Talma river basin has a catchment area of 42.10 Km 2 and has 10 rainfall- 
runoff events (Table 4.8 and Fig.4.8). The largest flood has a large lag time and 
yet the fitted regression equations are not much affected by it. The relationships 
are considered satisfactory. 

Khaga river basin has a catchment area of 66.00 Km 2 and has 5 storm flood 
events. The results are shown in Table 4.9 and Fig .4 .9. The results are quite 
consistent and are considered satisfactory. 

Mansari river basin is the largest basin considered in the study with a 
catchment area of 213.00 Km 2 . It has 5 storm flood events and the results are 
shown in Table 4.10 and Fig.4.10. Flood number 1 and 3 have the lowest peaks and 
have a lag time of 26 hours. Flood number 5 has the highest flood peak and has the 
lowest NK value of 20.25 hours. These indicate the general trend normally expected 
for the parameters. Flood number 4 having the second largest flood peak, however, 
indicates a very large NK value nearly twice that for flood 5 with which it differs 
in flood peak magnitude by about 3%. This is presumably due to errors, and the 
parameters of this storm were not used in considering the parameter variation 
within the basin. 

The consolidated list of the relationships between the model parameters and 
the flood peaks is given in Table 4.11 




QPEAK M**0'SEC 






TABLE 

4.8 RESULTS OF ANALYSIS USING NASH MODEL 




TALMA RIVER BASIN 


CATCHMENT AREA 42.12 SQ.KM 



DATE 

QPOBD. 

OB. TIME 

RFALL 

EFFRF 

BASEFLO 

PHIINDEI R-0 COEFFT 



M**3/SEC 

HOURS 

MM 

MM 

M**3/SEC 

MM/HOUR 


1 

2 

3 

4 

5 

6 

7 

8 

9 

1 

18190765 

16.340 

10.000 

32.004 

13.389 

5.125 

5.189 

0.418 

2 

220765 

27.240 

13.000 

85.344 

39.194 

6.315 

5,468 

0,459 

3 

23240765 

143.030 

13.000 

317.754 

212.835 

10.959 

7.876 

0.670 

4 

300765 

14.300 

17.000 

55.118 

11.564 

6.966 

6.390 

0.210 

5 

140865 

31.570 

11.000 

67.056 

17.620 

17.387 

7.754 

0.263 

6 

170767 

56.320 

12.000 

161.290 

60.139 

7.646 

11.442 

0.373 

7 

180767 

107.410 

11.000 

176.276 

65.512 

51.000 

14.227 

0.372 

8 

180968 

32.560 

14.000 

95.504 

35.946 

11.242 

25.715 

0.376 

9 20210968 

87.250 

11.000 

188.214 

64.513 

10.704 

23.792 

0.343 

10 

12130969 

56.320 

7.000 

55.000 

17.651 

24.126 

11.749 

0.321 


QPCOMP. 

COM. TIME 

METHOD OF MOMENTS 

N K 

NK 

QPCOMP. 

COMP. TIME 

METHOD OF LEAST SQUARES 

N K NK 

N* 

K* 

NK* 

M**3/SEC 

10 

HOURS 

11 

12 

13 

14 

M**3/SEC 

15 

HOURS 

16 

17 

18 

19 

20 

21 

22 

17.800 

12.000 

3.623 

2.905 

10.525 

16.370 

13.000 

3.661 

3.262 

11.942 

2.485 

5.211 

12.959 

27.900 

14.000 

2.869 

5.793 

16.620 

27.910 

14.000 

2.865 

5.793 

16.597 

2.485 

5.124 

12.733 

140.660 

16.000 

2.525 

5.385 

13.597 

135.690 

15.000 

2.151 

6.386 

13.736 

2.485 

4.198 

10.431 

15.640 

16.000 

1.481 

7.127 

10.555 

14.300 

16.000 

1.431 

8.959 

12.820 

2.485 

5.228 

12.991 

34.450 

11.000 

1.945 

4.293 

8.350 

31.820 

11.000 

1.641 

6.000 

9.846 

2.485 

5.089 

12.647 

58.920 

12.000 

1.983 

4.981 

9.877 

56.130 

12.000 

1.843 

5.634 

10.383 

2.485 

4.891 

12.155 

119.140 

11.000 

3.297 

2.602 

8.579 

107.920 

11.000 

2.688 

3.701 

9.948 

2.485 

4.483 

11.140 

40.510 

12.000 

3.347 

3.584 

11.996 

39.090 

12.000 

3.206 

3.884 

12.452 

2.485 

5.082 

12.628 

69.240 

10.000 

2.300 

4.129 

9.497 

78.380 

11.000 

3.172 

2.803 

8.891 

2.485 

4.644 

11.540 

63.730 

7.000 

2.445 

1.614 

3.946 

55.780 

7.000 

1.941 

2.438 

4.732 

2.485 

4.891 

12.155 



TABLE- 4.9 RESULTS OF ANALYSIS USING NASE KODBL 


KHA6A RIVER BASIN CATCHKENT AREA 66.00 SQ.KM. 



DATE 

QPOBD. 

OB. TIME 

RFALL 

EFFRF 

BASEFLO PHIINDEX 

R-0 COEF 



M**3/SEC 

HOURS 

MM 

MM 

M**3/SEC MM/HOUR 


1 

2 

3 

4 

5 

6 

7 

8 

9 

1 

150781 

241.020 

8.000 

132.250 

97.857 

16.590 

3.558 

0.724 

2 

240881 

129.500 

9.000 

149.750 

39.815 

12.000 

17.931 

0.266 

3 

280784 

45.930 

8.000 

58.400 

15.327 

4.430 

7.550 

0,262 

4 

140984 

62.130 

14.000 

161.500 

39.225 

5.510 

12.983 

0.243 

5 

160984 

53.900 

9.500 

55.100 

16.740 

16.590 

5.749 

0.304 




METHOD OF MOMENTS 




METHOD OF LEAST SQUARES 


ESTIMATED 


QPCOMP. 

M**3/SEC 

COM. TIME 
HOURS 

N 

E 

NE 

QPCOMP. 

M**3/SEC 

COM. TIME 
HOURS 

N 

E 

NE 

N* 

K* 

NR* 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

21 

242.650 

8.000 

2.922 

1.603 

4.684 

234.800 

8.000 

2.637 

1.820 

4.799 

2.847 

1.817 

5.173 

126.710 

9.000 

3.829 

1.643 

6.291 

118.400 

8.000 

3.910 

1.768 

6.913 

3.404 

1.928 

6 .566 

34.440 

10.000 

3.488 

2.172 

7.576 

35.600 

9.000 

3.093 

2.269 

7.018 

3.822 

2.012 

7.691 

64.000 

14.000 

3.192 

2.568 

8.197 

64.850 

14.000 

3.421 

2.140 

7.321 

3.741 

1.996 

7.467 

57.500 

9.000 

5.580 

1.248 

6.964 

52.500 

9.000 

4,557 

1.640 

7.473 

3.782 

2.004 

7.581 


T A B L E 4.10 RESULTS OF ANALYSIS USING NASH MODEL 




MANSARI ! 

RIVER BASIN 


CATCHMENT AREA 213,00 SQ.KM 



DATE 

QPOBD. 

OB. TIKE 

RFALL 

EFFRF 

BASEFLO 

PHIINDEX R-0 COEF 



K**3/SEC 

HOURS 

HM 

MM 

M**3/SEC 

MM/HOUR 


1 

2 

3 

4 

5 

6 

7 

8 

9 

1 

80670 

25.510 

32.000 

81.502 

12.760 

2.095 

16.330 

0.157 

r 

L 

26290970 

35.930 

23.000 

90.921 

13.231 

11.214 

17.761 

0.146 

*\ 

300771 

20.730 

22.000 

64.510 

7.595 

5.918 

11.597 

0.118 

4 

300971 

59.210 

23.000 

32.600 

15.993 

13.706 

4.488 

0.491 

5 

15190871 

57.400 

25,000 

158.671 

29.476 

17.359 

10.112 

0,186 



METHOD OF MOMENTS 




METHOD OF LEAST SQUARES 



ESTIMATED 


QPCOMP, 

COM. TIME 

N 

K 

NK 

QPCOMP. 

COM. TIME 

N 

K 

NK 

N* 

K* 

NK* 

M**3/SEC 

HOURS 




K**3/SEC 

HOURS 







10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

27.220 

30.000 

5.491 

5.514 

30.277 

27.600 

25.000 

4.011 

6.575 

26.372 

3.998 

6.826 

27.291 

40.570 

26.000 

5.558 

4.792 

26.634 

34.140 

26.000 

4.002 

7.550 

30.215 

3.998 

6.399 

25.583 

20.990 

28.000 

4.763 

5.993 

28.545 

21.240 

25.000 

4.004 

6.563 

26.278 

3.998 

7.022 

28.074 

58.120 

23.000 

4.440 

4.243 

18.839 

56.960 

24.000 

4.018 

5.041 

20.255 

3.998 

5.444 

21.767 

57.320 

32.000 

3.197 

10.696 

34.195 

58.730 

30.000 

2.999 

10.737 

32.200 

3.563 

4.225 

15.054 


Fig. 4. 10 MAN SARI BASIN- QPEAK VS N,K,NK 



QPEAK M**3/SEC 

a N + K « NK v NKF1TIED 
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Table 4.11 Regression relationships between model parameters and Q, 


Basin 

Regression Relationsips 

Batjhora 

N - 1.889 - 0.027 Q p 

K - 2.803 - 0.004 Qp 

Sankhini 

N - 2.597 + 0.006 Qp 

K - 2.900 + 0.026 Qp 

Mujnai 

N - 1.405 + 0.008 Qp 

K ■» 3.792 - 0.021 Qp 

Talma 

N - 2.485 

K - 5.342 - 0.008 Qp 

Khaga 

N - 4.502 - 0.005 Q p 

K - 2.053 - 0.001 Qp 

Dharsi 

N - 4.504 - 0.030 Qp 

K - 1.800 + 0.056 Qp 

Mansari 

N - 3.998 

K - 7.872 - 0.042 Qp 


4.S.2 Storm runoff parameters 

Perusal of Tables 4.3 to 4.10 show that the infiltration index, the coefficient 
of runoff and the baseflow vary from storm to storm in each basin and is often 
erratic. Attempts were made to correlate them with storm and runoff parameters; 
but they were unsuccessful. This indicates the inherent variability of the process. 
However, the Tables indicate that the basef low/unit area may be a better parameter 
than the baseflow itself. 

4.6 Reg iona 1 i sat ion of Rainfall-Runoff Parameters 

From the analysis of the storm flood events in the various basins, parameters 
are to be identified as representative of the region under study for use in flood 
estimation for the catchments in the region. Table 4.12 show the ranges of these 
parameters for the various basins. 
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Table 4.12 Ranges for storm runoff parameters 


Sl.No. 

Basin 


Ranges 

For 


Infiltration 

index 

mm/hour 

Baseflow 
per unit area 
m 3 /s©c/Km 2 

Coeff. 

of Runoff 

1. 

Batjhora 

0.50 - 15.70 

0.22 - 0.77 

0.12 

- 0.95 

2. 

Sankhini 

5.00 - 26.70 

0.94 - 0.97 

0.14 

- 0.40 

3. 

Mujnai 

7.00 - 13.00 

0.27 - 0.64 

0.13 - 

0.42 

4. 

Talma 

5.00 - 26.00 

0.12 - 1.20 

0.21 - 

0.67 

5. 

Khaga 

3.50 - 18.00 

0.07 - 0.25 

0.24 - 

0.72 

6. 

Dharsi 

2.40 - 33.00 

0.13 - 0.31 

0.10 - 

0.76 

7. 

Mansari 

4.50 - 17.80 

0.01 - 0.08 

0.12 - 

0.49 


As seen from the Table, the coefficient of runoff varies from 0.1 to 0.95. Since 
variation is too large and no trend is observed for its variation, this parameter 
cannot be used for regionalisation. The infiltration rate in the basin is another 
parameter that could be considered for regionalisation and for the design flood a 
minimum value of the infiltration rate which will lead to a large flood can be taken 
as a design parameter. 5.00 mm/hour seems to be the minimum value for the basins 
in general. This value corresponds to a silty loam soil or sandy clay loam soil 
which is considered realistic. So this value can be taken for computation of design 
flood for the catchments in the region. This value can be considered to be the 
general minimum value for the phi-index for all the basins, particularly during the 
intensive storms generally met with in the basin. Accordingly an abstraction index 
of 5mm/hour is adopted for all the basins in the study. . 

The baseflow generally increases with basin area. The range of variation of 
baseflow/unit area for each of the basins is also shown in Table 4.12. An enveloping 










93 


iii) Estimate the effective rainfall for the storm assuming the minimum phi-index 
of 5.00 mm/hour; 

iv) Estimate the Nash model parameter for the basin using Eqns.4.10 and 4.11; 

v)Calculate the IUH ordinates in mm/hour from Eqn. 4.1; 

vi) Estimate the 1-hour UH in mm/hour; 

vii) Convert the UH values to m 3 /sec; 

viii) Calculate the DSRO ordinates by convolving the 1 hour effective rainfall 
values with the UH values; 

ix) Estimate the baseflow from Fig. 4.12; and 

x) Calculate the flood discharge values by adding the baseflow to the DSRO 
ordinates. 

To demonstrate the use of the model developed in this Chapter, the Nash 
model is applied to Kurti basin in the region, the data for which were not used in 
the study. 

4.7.1 Design storm 

For Kurti basin, Kirpich equation (Eqn .2.3) is used to estimate the time of 
concentration (t c ) with L - 14.81 Km. and S - 0.03166 as equal to 118.63 minutes, 
say, 2 hours. 

Depth -duration -frequency (DDF) maps developed by IMD (1974) for 50 year 
return period are used to estimate the design storm. The basin is located in the 

maps and the rainfall for various durations of 30 minutes, 1 hour, 3 hours, etc., 

upto 24 hours were noted from the respective DDF maps. Point rainfall depths for 
50 years frequency were plotted as a function of duration on a log-log paper 
(Fig.4.13). The 2 hour precipitation value is interpolated as 105 mm. The point 
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rainfall values are converted to areal rainfall values by multiplication with area 
reduction factors (CWC, 1973). The hourly storm values are then derived. 
CWC/IMD adopts a sequencing procedure for the design storm precipitation, viz., 
the highest hourly precipitation is kept in the middle, the next highest value 
follows it, the third precedes it and the others alternate. Since the design 
hyetograph has a duration of 2 hours only, it has hourly precipitations of 67.22 mm 
and 19.57 mm (Table 4.13). 


Table 4.13 Calculation of Design Storm Hyetograph 


Duration 

Point 

Areal to 

Areal 

Hourly 

Abstrac- 

Effective 


Rainfall 

Point 

Rainfall 

Rainfall 

-tion 

Rainfall 


from DDF Rainfall 


of Design 




Curve 

Ratio 


Storm 




mm 


mm 

mm 

mm/hour 

mm 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(7) 

30 min. 

65 

68.33 

46.46 

- 

_ 

— 

1 hour 

88 

76.39 

67.22 

67.22 

5.00 

62.22 

2 hours 

105 

82.66 

86.79 

19.57 

5.00 

14.57 


4.7.2 Design unit hydrograph 

The Kurti basin has a catchment area of 21.91 Km 2 . Using Eqns.4.10 and 4.11, 
the values of the Nash model parameters for the basin are found to be N — 1.9855; 
and K — 3.1170. Table 4.14 shows the estimation of design flood for the basin. The 
IUH ordinates at hourly intervals are calculated from Eqn.4.1. The 1 hour UH 
ordinates are estimated by integration of 1 minute IUH values. Since the area of the 
basin is 21.91 Km 2 , the corresponding values in m 3 /sec are shown in column 3. The 
effective rainfall values for the design storm are convolved with the 1 hour UH 
values to yield DSRO ordinates (column 4). The baseflow of 21.91 m 3 /sec is 






Fig.4-13 Depth duration curve tor 50 year return period 
tor the region 
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ESTIMATION OF DESIGN FLODD FOR KURTI BASIN 
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estimated for the basin from Fig. 4.12. Using this value, the flood discharge values 
are computed and are shown in column 5. The flood peak of 76.36 m 3 /sec occurs at 4 
hours for this basin. 

For small basins, it is desirable to derive the design storm at shorter 
intervals, say 15 to 30 minutes and use it to derive the design flood. The design 
flood for 15 minute interval has been derived and is shown in Table 4.14a. The 
results indicate that the hydrograph peak is changed marginally (from 76.36 m 3 /sec 
to 77.21 m 3 /sec) but its time of occurrence is earlier by 45 minutes. 

4.7.3 Comparison of results with Snyder’s approach 

Snyder’s approach (Eqns.2.12 to 2.15) requires the values of the coefficients 
C t and C p . C t and C p are computed for a gauged watershed nearby the ungauged 
one. In the present study, the coefficients are computed for the Mujnai basin. From 
Eqns. 4.10 and 4.11, the Nash model parameters for this basin are calculated as 
N- 2.131 and K = 3.382 . The 1 hour UH ordinates are obtained (Table 4.15). From the 
catchment map of this basin, the values of L, L c and S are obtained as 21.33 Km, 
11.75 Km and 0.03163 respectively. The peak discharge of 1 hour UH for the basin 
is 1.0863 m 3 /sec and it occurs at t = 4 hours. Using these values, the Snyder’s 
coefficients C t and C 3 are calculated. For the Mujnai basin, t R = l hour; 
t pR = 3.5 hours. Hence t r = 13/21 hours; t p - 71.5/21 hours; q pR - 0.02822 m 3 /sec/Km 2 . 
Thus C t - 0.888 and C p = 0.0359. Using the above values of C t and C p for the Kurti 
basin, t R » 1 hour; t pR - 2.98 hours; q pR = 0.0331 m 3 /sec/Km 2 and Q PR = 0.726 m 3 /sec 
compares well with t pR - 3.5 hors and Q PR = 0.709 m 3 /sec derived using the UH 
approach in the study. 
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4.7.4 Comparison with the method proposed by CWC 

CWC has derived a set of equations for UH parameters for this region for 
larger catchements (CWC,1991) and these equations (Eqns. 2.17 to 2.24) have been 
given in Subsec 2.5.1. The values of peak discharge/unit area (q p ) , time to peak 
(t p ) and peak discharge (Q p ) computed using the equations are : 

q p - 0.7664 m 3 /sec/Km 2 ; t p = 2.779 hours; and Q p — 16.79 m 3 /sec. 

The values of Q p obtained by the CWC approach is much higher than those 
obtained by the approach suggested in the present study or the Snyder’s approach. 

The set of equations derived by CWC is essentially based on data available 
for larger catchments and these do not seem to apply to the smaller catchments. 
The value of LL C /S for this basin is 14.25. From the curve developed by CWC 
between LL C /S and q p , to get a value of q P of the order of 0.7 m 3 /sec/Km z , the 
value of LLc/S has to be of the order of 5000 which means that the slope S has to 
be very small, which is not realistic. The approach proposed by CWC, thus, seems 
to be valid only for larger basins and not smaller catchments. 



CHAPTER 5 


USGS DISTRIBUTED ROUTING RAINFALL-RUNOFF MODE 

S.l Introduction 

The distributed parameter model can generally take into account the spatial 
and the topographic characteristics of the watershed. Hence it is a more realistic 
representation of the watershed. There are a large number of distributed parameter 
models for the watersheds (Fleming, 1975, Beven,1985). A simple model to simulate 
runoff hydrographs was developed by the Water Resources Division of the United 
States Geological Survey (Dawdy et al ,1978). The model known as the Distributed 
Routing Rainfall Runoff Model (DR3M) is deterministic and it combines the soil 
moisture accounting and rainfall excess components developed by Leclerc and 
Schaake *(1973). Wherever possible, a physical interpretation is placed upon the 
parameters used in the model. Though developed essentially for routing urban flood 
discharges through a branched system of pipes to natural channels, Murthy (1988) 
used this model for simulating flood hydrographs from a semiurban and rural 
catchment and a hilly catchment and concluded that this model is reasonable and 
realistic for such watersheds. Furthermore, it is a relatively simple model with 
physically based parameters. Hence it is proposed to use the USGS model in the 
study and the following description generally follows Dawdy et al (1978). 

Inputs to the model include daily rainfall and rainfall during unit rainfall 
days on which major storm precipitation occurs, the duration being a multiple of 5 
minutes from atmost three raingauge stations; daily pan evaporation, and a physical 
definition of the basin discretised into as many as 50 segments including overland 
flow, channel and reservoir segments. The model maintains a unit duration soil 
moisture accounting during rainy days and daily accounting on unit rainfall days. 
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5.2 Structure of the Model 

The model is divided into four major components: a soil moisture or water 
balance component; rainfall excess component; a routing component and an optimising 
component. The overall flow chart of the model is given in Fig. 5.1. 


5.2.1. Soil moisture accounting component 

The schematic diagram of the soil moisture accounting component is shown in 
Fig. 5.2. This component determines the effect of antecedent conditions on 
infiltration. Soil moisture is modelled as a two layer system, one representing the 
base moisture storage (BMS) and the other upper wetted part caused by infiltration 
into soil moisture storage (SMS). 

During unit rainfall days, moisture is added to SMS based on Philip’s 
infiltration equation (Philip, 1954) 

4i - k ( i + ) ts.i) 

&T 1 


where i 

k 

P 

m Q 

m 


accumulated infiltration volume in wetted soil column since the 
start of the infiltration; 
capillary conductivity of soil; 

capillary potential at wetting front in soil column at the time at 
which the infiltration rate is being computed; 
uniformly distributed moisture content of soil column when 
infiltration started; and 

moisture content, uniformly distributed through wetted column at 
the time at which the infiltration rate is being computed. 




Fig* 5 .I Flow chart for DR3M (Adapted from Dawdy et al*,1978) 












RAINFALL 

INPUT 



Fig. 5.2 Schematic flow chart of soil-moisture accounting component 
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It is assumed that at the beginning of infiltration, the entire soil column is at a 
uniform moisture content m 0 and at any time during infiltration, the moisture 
content m above the wetting front is also uniform. Eventhough these assumptions 
do not give exact depth of the wetting front from the surface, they appear to give 
useful results (Philip, 1954). 

On other days, a specified proportion of daily rainfall (RR) infiltrates into 
the soil. If the daily precipitation is less than the daily infiltration rate, the 
infiltration rate is reset equal to the the daily precipitation rate. Irrigation can be 
accounted for in the daily water balance as an additional input through supplied 
irrigation rates for each month. 

Evapotranspiration takes place from SMS based on availability and otherwise 
from BMS with the rate determined from pan evaporation multiplied by a pan 
coefficient(EVC). Moisture in SMS drains into BMS with a controlling parameter 
(DRN) determining the rate. Storage in BMS has a maximum value (BMSN) equal to 
the field capacity moisture storage of the active soil zone. Zero storage in BMS is 
assumed to correspond to wilting point condition in the active soil zone. When 
storage in BMS exceeds BMSN, the excess is spilled to a deep storage. These spills 
could be the basis for routing interflow and baseflow components, if desired. 
However, this option is not included in the present model. Table 5.1 shows the 
various parameters for soil moisture accounting. Since the original data were 
generally in FPS units and the program was also in FPS units, it was decided not to 
modify the units in the program. 
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Table 5.1 Parameters for Soil-Moisture Accounting and Infiltration 


Soil-Moisture Accounting 


Parameters 

DRN 

EVC 

RR 

BMSN 


A constant drainage rate for redistribution of soil moisture 
between SMS and BMS, in inches per day 

A pan coefficient for converting measured pan evaporation to 
potential evapotranspiration 

The proportion of daily rainfall that infiltrates into the soil 
for the period of simulation excluding unit-rainfall days 
Soil-moisture storage at field capacity, in inches. 


Infiltration 


Parameters 


KSAT 

- The effective saturated value of hydraulic conductivity, in 

inches per hour 

RGF 

- Ratio of suction at the wetting front for soil moisture at 

wilting point to that at field capacity 

PSP 

- Suction at wetted front for soil moisture at field capacity, 

in inches of pressure 


5.2.2 Rainfall excess component 

i) Impervious surface: Two types of impervious surfaces are considered by 
the model, namely, effective impervious surfaces and noneffective impervious 
surfaces. The effective impervious surfaces are those which are directly connected 
to the channel drainage system such as streets and roofs which drain onto 
driveways and streets. The noneffective impervious surfaces are those which drain 


onto the pervious areas. 
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The only abstraction from rainfall on effective impervious area is impervious 
retention, which is to be specified and must be filled before runoff from effective 
impervious areas can occur. Evaporation occurs from impervious retentions during 
periods of no rainfall. 

Rain falling on non-effective impervious areas is assumed to run off to the 
surrounding pervious area. The model assumes that this occurs instantaneously and 
that the volume of runoff is uniformly distributed over the contributing pervious 
area. This volume expressed as inches over the pervious area is added to the rain 
falling on the pervious area prior to the computation of pervious area rainfall 
excess. 

ii) Pervious surfaces: The point potential infiltration (FR) is computed by 
Philip’s equation (Philip, 1954) given by 

FR - KSAT (1 + (5.2) 

where KSAT is the effective saturated soil capillary conductivity and PS is the 
capillary potential at the wetting front. PS varies over the range from the field 
capacity to wilting point linearly as (Fig. 5.3) 

PS - PSP (RGF - ~ RQ g^ BMS ) (5.3) 

where PSP is the suction at wetting front at field capacity, RGF is ratio of suction 
at wilting point to that at field capacity and BMSN is the effective soil moisture 
storage at field capacity. 

Point potential infiltration (FR) is converted to effective infiltration over 
the basin using a scheme for area variability of infiltration, first presented by 
Crawford and Linsley (1966). If SR represents the supply rate of rainfall for 
infiltration and QR represents the rate of generation of rainfall excess, the area 




(BMS-O) (BMS=BMSNi 


SOIL- MOISTURE CONTENT 


Fig. 5.3 Variation of soil-moisture potential 



PERCENTAGE OF AREA WITH INFILTRATION CAPACITY 
EQUAL TO OR LESS THAN INDICATED VALUE 


Fig. 5.4 Area variability model for infiltration 
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variability equations are 

QR - SR 2 /2FR SR < FR (5.4 a.) 

QR - SR - (FR/2) SR > FR (5.4 b.) 

Fig. 5.4 shows the schematic representation of these relationships. The infiltration 
parameters, described above, have been summarised in Table 5.1. Two different soil 
types can be handled by the model with separate soil moisture accounting and 
infiltration parameters for each soil type. 

5.2.3 Routing component 

A drainage basin is represented in this model as a set of segments which 
jointly describe all sub basins in the total basin. There are four basic types of 
segments: (i) overland flow segment, (ii) channel segment, (iii) reservoir segment, and 
(iv) nodal segment. The basin is divided into a number of segments in such a way 
that the essential basin geometry and physiography are taken into account for 
runoff computation. 

(i) Overland flow and channel segments: A channel segment is permitted to 
receive upstream inflow from as many as three other segments, including 
combination of other channel segments, reservoir segments and nodal segments. It 
may also receive lateral inflow from as many as four overland segments. The 
overland flow segments receive uniformly distributed lateral inflow from excess 
precipitation. 

Kinematic wave theory is used for both overland flow and channel routing. 
The partial differential equation to be solved for each channel and overland 
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segment is 


SA . SQ _ 

6t 6x 


q 


(5.5) 


where A is the area of flow, Q is the rate of flow, q is the rate of lateral flow, t is 
time and x denotes the distance along the segment increasing in the downstream 
direction. The dependent variables A and Q are functions of two independent 
variables x and t. The relationship between A and Q is given by 

Q - a A m (5.6) 

The length of segment is L and the outflow hydrograph from the segment is Q(L,T). 
The inflow hydrograph to the upstream end of the segment is Q(0,t) and these are 
the boundary conditions needed to solve Eqn.5.5. In case there is an upstream 
inflow, the solution will also depend on that inflow. Therefore, the solution must be 
some function of x, t, q and Q(0,t). The outflow hydrograph is given by 

Q (L,t) - f ( Q (0,t), q(t) ) (5.7) 

The two parameters a and m of the kinematic wave model can be determined by 
analysis of the physical characteristics of the basin segments. For example, if 
Chezy’s equation is employed to quantify friction, then 

a - C (Sq) 0 6 and m - 3/2 (5.8 a) 

where C is the Chezy’s friction factor and S 0 is the slope. If Manning’s relationship 
is used, 

a - i (So) 0 5 (5.8 b) 

where n is the Manning’s friction factor. Functional relationships of cx. and m for 
channel and overland flow segments are listed in Table 5.2. 

Numerical techniques are used to approximate Q(x,t) at discrete locations in 
the (x,t) plane. In this case, the rectangular grid of points was selected at intervals 



Table 5.2 Relationships for estimating kinematic wave parameters a and m 


Type of 
segeent 

ITYPE(l) 

<*: 

ft 

Def i ni tioi 
of PARAMd, 

Rectangular 

conduit 

1 

1*49 JfLOPEm 

F RN7 1 7 (PARAMd, ir > 

1.67 

width of 
conduit 

Circular pipe 

2 K49 

FRN(I) 

2/3 

(PARAM ( I , 1 ) /4) SLOPE ( I ) 

1.0 

diameter of 

pipe 

Triangular 
cross-secti on 

3 

K49 JSLOPfjrr 

FRN ( 17 lPARAMd,n ) 

1.33 

width at 
1-foot depth 

External specifi- 
cation of and a 

4 

PARAM (1,1) 

PARAMd , 

>2) 

Overland flow 
(turbulent) 

5 or 15 

1. 49-JsLOPE ( I ) 

FRN ( I ) 

1.67 


Overland flow 
(laeinar ) 

6 or 16 

66. 4*JSL0PE ( I ) 

0. 0000 141*FRN ( I ) 

3.0 



FRN(I) * Friction coefficient for segaent 1 
SLOPE(I) * Slope of segaent I 
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of time At and distance Ax. The value of Ax varies from segment to segment but 
the value of At is constant for all segments. Four points of a finite difference mesh 
are shown in Fig.5.5. The purpose of the finite difference scheme is to solve for A 
and Q at point d, given the values of Q and A at points a,b, and c. 

The model contains two different finite difference equations and selects the 
appropriate one at each point in the solution. This is done in an attempt to keep the 
solution errors small while maintaining an unconditionally stable solution. 


t-FAt - 


t + 



X- AX X 

Fig. 5 A Four point finite - difference mesh 
The decision depends upon the parameter 

0 - m ^ ^ . am A 1 "- 1 (^) 

Ax Ab Ax 


-►X 


(5.9) 


If 0 > 1, the equations used are 

Qd “ Q c + q Ax - Ax (Ac - AJ/ At , and 
A« - (Qd/a) 1/m 


(5.10) 

(5.11) 


This involves only mesh points a, c and d. It was derived by substituting (A c -A a )/At 
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for 6 A/6t and (Q d -Q c )/Ax for 6 Q/6x. 

If 0 < 1, the equations used are 

A d « A b + qAt + At (Q. + Q b ) , and (5.12) 

Q d - a A d m (5.13) 

Initial values of A and Q are given along the entire X-axis. At t - 0, the 
model sets A -= 0 and Q - 0 everywhere. During solution, the upstream inflows are 

given and Eqn 5.6 is used to compute the upstream boundary condition for A. 

Equations 5.10 to 5.13 are solved by the model beginning with x = 0 and proceeding 
downstream at intervals of Ax to x = L. 

5.2.4 Selection of Ax and At 

Dawdy et al (1978) use the following approach to select Ax and At. Two 
factors are important for At. One is the frequency of rainfall input to the 
catchment and the other is the frequency reponse characteristics of the catchment. 
For a basin having overland flow length L 0 , overland kinematic wave parameter a 0 
and m 0 , channel length Lc and channel flow parameters a c and m c , 

At x 0.1 (to + t c ) (5.14) 

where 

to - [ — 7 ] 1/m ° seconds (5.15) 

oto (ie / 43200) m °" 1 

and t c - [ ^ jl/iDo seconds (5. 16 ) 

<x c ( N U L o /43200) nic " i 

where i e is the maximum rainfall intensity in inches/hour, N is the number of sides 

of channel with overland flow segments contributing lateral inflow (1 or 2) and L c 
and L 0 are in feet. 

After At is selected, Ax is chosen for each segment. This is done by 
specifying the number NDX of Ax segments in each segment. NDX may vary from 


segment to segment. 



The finite difference solution will be an exact solution if Ax and At are 


selected so that characteristics passing through point a (Fig.5.5) also passes through 

point d. In such a case, solution A d and Q d depends only on A a and Q a and lateral 

inflow along the characteristics curve between a and d. In the special linear case 

when m - 1, an exact solution may be obtained by selecting the ratio ^ - a*- 0 

am A™' 1 . Hence Ax =* a*At. a* is the effective value of a that would be used for 

a linear approximation to a given kinematic flow segment. For overland flow, 

a 0 * = L 0 /t 0 ; For channel flow, a* - L c /t c . Thus, Ax 0 =^ L 0 At/Ax ; 

Ax c a- L c At/Ax; and NDX 0 = and NDX C 

Ax 0 At Ax c At 

In general, the smallest subarea of interest and the highest intensity of rainfall 
should be used to estimate At. 

5.2.5 Optimisation component 

An option is included in the model to calibrate the soil moisture and 
infiltration parameters for drainage basins having rainfall runoff data. The method 
of determining optimum parameter values is based on the technique devised by 
Rosenbrock (1960). The objective function is the sum of the squared deviations of 
the logarithms of computed and measured storm runoff volumes. To start the fitting 
process, the model is assigned an initial set of parameter values and upper and lower 
bounds of each parameter. The suggested values by Dawdy at al are shown in Table 
5.3. The estimation of all the parameters simultaneously is complicated and 
unnecessary. The procedure involves estimation of parameters in a number of 
phases. Initially some of the less sensitive parameters are assumed constant and the 
more sensitive parameters are optimised. The process itself may be iterated in a 
number of phases. For example, the last three parameters (Table 5.3) were found to 
be less sensitive and assuming them as constant, the remaining parameters are 
optimised. While optimising the four parameters, each of them can be 
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Table 5.3 Initial values of soil-moisture and infiltration parameters and their bounds 


Parameter 

Initial 

value 

Lower 

limit 

Upper 

limit 

Units 

PSP 

5.00 

1.00 

15.00 

Inches 

KSAT 

0.10 

0.01 

1.00 

Inches/hour 

RGF 

10.00 

5.00 

20.00 

Dimensionless 

BMSN 

5.00 

1.00 

15.00 

Inches 

EVC 

0.70 

0.50 

1.00 

Dimensionless 

RR 

0.90 

0.65 

1.00 

Dimensionless 

DRN 

0.50 

0.10 

1.00 

Dimensionless 


optimised separately keeping the rest constant till an acceptable result is obtained. 
The optimisation procedure revises the parameter magnitudes and recomputes the 
objective function using the revised set of parameter magnitudes. If the result is an 
improvement, the revised set is accepted; if not, the method returns to the previous 
best set of parameters. 

It is found during the studies that the results depend very much upon the 
initial values of the parameter set and hence, the initial values should be chosen 
with care or alternate sets of initial values may be tried. 


5.3 Implementation of the Computer Program 

The original program developed by. USGS Water Resources Division is written 
in FORTRAN IV and is available from the User’s Manual for the model (Dawdy et 
al,1978). This program was implemented in the microVAX II computer system of the 
Centre for Water Resources Engineering and Management at IIT Kanpur by Murthy 
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(1988) and was available for the study. Some minor modifications were to be made 
while implementing the program, since the original program was developed for CD C 
7600 machine. 

5.4 Analysis Using DR3M 

DR3M requires the basin to be divided into a consistent set of overland flow 
and channel segments. A main channel runs through the entire length of the basin 
with its tributaries. For convenience, the main channel is divided into a certain 
number of segments and the tributaries are treated as separate channel segments. 
While delineating these segments the number of segments is to be balanced; 
uniformity within segments is to be maintained as far as possible; segmentation of 
overland flow is to be done with balanced length of overland flow to the designated 
channel segments; and compatibility between overland flow segments and the 
designated channel segments is also to be maintained. The area of each segment is 
planimetered and slopes are determined from contours in a topo sheet of the basin. 
The length of overland flow is computed by dividing the area of each segment by 
the length of the channel segment into which it contributes as a lateral flow. 
Crosssections of the river channel in different reaches are required for the 
analysis. The crosssection data are generally available only at the bridge sites for 
the basins selected for the study. Two options are available for the roughness 
coefficient for the overland flow segments. A single roughness coefficient can be 
assigned to an overland flow segment or else the segment can be assigned two 
roughness coefficients, one for the pervious and the other for the impervious 
surfaces. 

Direct impervious areas in the basin are to be estimated with care because 
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they contribute directly to the channel without any loss (Rao and Ramaseshan, 1978). 
In the absence of data, particularly of rocky and hilly impervious areas in these 
basins, the impervious area is initially assumed to be negligible and ignored. An 
attempt is made to parameterise this factor in the final stages of the study. In case 
data were available, it will be desirable to identify the spatial distribution of such 
areas and use them appropriately in the model. 

5.4.1 Analysis of a basin using DR3M 

To illustrate the details of the analysis of a basin using DR3M, the details of 

the Talma river basin is presented. This basin (Fig.5.6) has a catchment area of 

* „ 

42.12 Km and has one raingauge station. The data for this basin are available 
generally from 15 May to 15 October for the years 1965 to 1969. Since evaporation 
data were not available for any station within the basin, pan evaporation was 
calculated for Guwahati which is the nearest station for which data were available 
and considered reasonably representative of the basin. Christiansen’s equation was 
used in the estimation of monthly pan evaporation with monthly normal values of 
the meteorological characteristics of Guwahati. A pan coefficient of 0.7 was used to 
estimate the monthly potential evaporation. Daily normal potential evaporation was 
used in the model. 

The basin is divided into a limited number of overland flow and channel 
segments (Fig.5.6). A main channel runs through the entire basin and there is a 
major tributary. The main channel is divided into two channel segments CH20 and 
CH21 respectively and the major tributary is designated as channel segment CH22. 
On either side of the channel segments, there are overland flow segments and they 
are designated as OFOl to OF06. Direct lateral inflow into downstream channel 
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section is ignored. Thus the entire basin is divided into 9 segments consisting of 3 
channel segments and 6 overland flow segments. The details of these segments, for 
example, the length and slope of the channel segments and the length, area and 
average slope of the overland flow segments are determined from topo sheet. The 
length of overland flow was computed by dividing the area by the length of the 
channel segment into which it contributes as lateral inflow. 

The crosssection of the river channel in different reaches are required. 
However generally details of only one crosssection, viz., that at gauging station was 
available. While it may be possible to measure the crosssectional details at other 
upstream sites, it was not possible to do so in the study because of limitations of 
time, resources and remoteness of the region. From a comparison of a large number 
of crosssections of small rivers in the regions, approximate cross- sections for 
upstream channel segments were adopted in the study. The channel and overland 
flow segment characteristics of the basin are shown in Table 5.4. 

Since rainfall and runoff data are available from 15 May, and May is a hot 
month with very limited rainfall, the simulation is started every year from the 
beginning date of record with the initial condition of the basin corresponding to a 
very dry condition. A warming up period of about 15 days in the month of May 
brings the soil moisture, baseflow etc., to the realistic conditions and the model fits 
only important floods generally between July and September. In the case of Talma 
basin, initially 12 storm flood events were identified for use in the study (Sec.3.5). 



117 


Table 5.4 Segment characteristics of Talma basin 


Channel characteristics 

Segment 

Upstream 

Length 

Slope 

Width 

Manning’s n 

No. 

segments 

feet 

ft/ft 

ft 


CH20 


29040 

0.007 

80 

0.035 

CH21 


27456 

0.007 

90 

0.(85 

CH22 

CH20.CH21 

12672 

0.009 

no 

0.040 

Overland flow characteristics 

Segment 

Channel segment 

Length 

Slope 


No. 

for drainage 

feet 

ft/ft 


OFOl 

CH20 


4439 

.0047 


OF02 

CH20 


3351 

.0063 


OF03 

CH21 


3223 

.0047 


OF04 

CH21 


2486 

.0060 


OF05 

CH22 


2127 

.0056 


OF06 

CH22 


3291 

.0037 


It 

should be noted 

that in 

this model, the simulation is 

done not for a storm 


event but for all the storms in a year sampled continuously over time. The initial 
simulation used the parameters and bounds (Table 5.3) suggested by Dawdy et al, 
viz., PSP- 10; KSAT - 0.10; RGF-10; BMSN - 5.0; EVC - 0.70; RR-0.90; and DRN - 
0.50. The results of the analysis are shown in columns 5 and 6 in Table 5.5. A 
comparison of observed and simulated floods indicate that the simulated volume of 
runoff are generally much higher and the simulated peaks are also higher except for 
small flood events 2, 3, and 5. The time to peak of the simulated flood events is 
comparable to that of observed data except for storm events 3, 4, 7, 9, 11 and 12. 







TABLB 5,5 RESULTS OF ANALYSIS FOR TALKA BASIN 
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Events 3 and 12 are small floods and flood event 9 has multiple peaks. However, 
time to peak is not a major parameter in this study. 

There are seven major parameters in the model that can be fitted. Some of 
these can be treated as regional constants or may be evaluated separately from 
hydrograph analysis, seasonal water balance etc. Since the model is relatively 
insensitive to the three soil moisture accounting parameters, DRN, EVC, and RR, 
they are considered to be constants. Based on the soil characteristics, DRN, the 
constant drainage rate for redistribution of soil moisture between SMS and BMS of 
0.5 inches/day; EVC, the pan coefficient for converting measured pan evaporation to 
potential evapotranspiration of 0.70; and RR, the proportion of rainfall that 
infiltrates into the soil for the period of simulation excluding unit rainfall days of 
0.90, are adopted in the study. It is proposed to consider in greater detail the 
other four parameters in the modelling of the small watersheds. 

Initially all the 12 storm flood events were considered to be equally 
important and the sum of the squared errors between the observed and simulated 
floods for all the flood events was considered to be a measure of goodness of fit of 
the model. Using the optimisation option, the four parameters were optimised from 
the initial values. A summary of the results of this analysis is shown in column 8 
and 9 of Table 5.5. The results are poor for the relatively small floods 2, 3, and 5 
as in the case of initial parameters and are otherwise much better for the larger 
floods. The values of the parameters and the initial and final values of the 
objective function are also shown in the Table. The objective function decreases 
slightly from around 8.7 to around 5.8. 

A number of studies were carried out to understand the relative importance 
and sensitivity of each of the four model parameters. They included the 
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optimisation of various combinations of parameters and for different combination of 
floods, say all large floods and all small floods. The results indicate that the 
volume of runoff is significantly affected by the parameter KSAT and the results 
were generally less sensitive to PSP, RGF, and BMSN. Further studies are reported 
later. 

5.4.2 Analysis of all the basins using DR3M 

For all the basins used in the study, a few assumptions were made. Since no 
evaporation data were available, the information available at Guwahati was used 
(Sec. 5.4.1) and the monthly potential evaporation was estimated for each of the 
basins using the monthly normal values of the meteorological characteristics of 
Guwahati. Daily potential evaporation values were used in the model. The cross- 
sectional details were available only at the bridge site. Approximate crosssections 
for the various upstream segments were adopted in the study from a comparison of 
a large number of crosssections of small rivers in the region. 

i) Bat jhora Basin : This basin (Fig.5.7) has a catchment area of 7.33 Km i) 2 and 
the data were available for the years 1964 and 1965. The basin has one raingauge 
station. The main channel is divided into 4 channel segments , CH20, CH21, CH22, 
and CH23 and there are 8 overland flow segments, one either side of each of the 
channel segments and designated as OFOl to OF08. 

There are 9 storm flood events in the basin. Assuming the initial values of 
the parameters as suggested by Dawdy et al, simulation was carried out and the 
results are shown in Table 5.6. A comparison of observed and simulated floods 
indicate that the simulated volumes of runoff are generally much higher and the 




TABLE 5.6 RESULTS OF ANALYSIS FOR BATJHORA BASIN 


jOOD OBSERVED VALUES NO OPTIMISATION 
4BER 


PEAK EFF.RAIN PEAK EFF.RAIN 


(1) 

(2) 

(3) 

(4! 

(5) 

1 

160 

0.356 

365 

0.500 

2 

399 

1.322 

2601 

2.888 

3 

310 

0.848 

609 

1.153 

4 

233 

0.590 

466 

0.804 

5 

310 

1.399 

1390 

1.324 

6 

813 

3.779 

1402 

3.137 

7 

813 

2.470 

891 

2.332 

8 

163 

0.567 

912 

1.424 

9 

194 

0.385 

745 

0.824 




PSP 





KSAT 





RGF 





BMSN 



OBJECTIVE 


INITIAL 

2.380 


FUNCTION VALUES 

FINAL 




S I M U L 

A T E D 

V A L U 

OPTIMISATION WITH 

OPT. WITH PSP=10 

INITIAL VALUES 

RGF 

= 18 

ALL FLOODS 

ALL FLOODS 

PEAK 

EFF.RAIN 

PEAK 

EFF.RAIN 

(6) 

(7) 

(8) 

(9) 

223 

0.315 

229 

0.321 

2443 

2.574 

2464 

2.614 

424 

0.781 

469 

0.848 

311 

0.541 

273 

0.478 

1194 

1.043 

1116 

0.955 

1286 

2.550 

1260 

2.440 

805 

2.044 

790 

1.997 

631 

1.046 

567 

0.939 

248 

0.395 

250 

0.398 


7.586 


10.000 


0.112 


0.102 


16.460 


17.774 


1.890 


2.400 


2.380 


1.222 


1.126 


1.157 


S 

OPT. KITH PSP=10 OPT. WITH PSP=10 
R6F=18 ;BMSN=10 

ALL FLOODS ALL FLOODS 


PEAK 

EFF.RAIN 

PEAK 

EFF.RAIN 

(10) 

(ID 

(12) 

(13) 

270 

0.37 

229 

0.321 

2543 

2.673 

2464 

2.614 

531 

0.963 

469 

0.848 

349 

0.632 

278 

0.486 

1130 

0.986 

1116 

0.956 

1288 

2.568 

1260 

2.44 

812 

2.064 

790 

1.998 

574 

1.001 

566 

0.94 

269 

0.432 

361 

0.498 


10.000 


10.000 


0.096 


0.102 


18.000 


18.000 


1.180 


10.000 


0.300 


1.226 


0.220 


1.215 


TPEAK = TIME TO PEAK IN HOURS 
EFF.RAIN s EFFECTIVE RAINFALL 


PEAK = PEAK DISCHARGE IN CUSECS 
IN INCHES 
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simulated peaks are also high for all the floods indicating that the system 
parameters are different from the initial values. Keeping the three soil moisture 
accounting parameters (EVC, RR, DRN) constant, the remaining four parameters, viz., 
the infiltration parameters PSP, KSAT and RGF and the soil moisture accounting 
parameter BMSN were optimised from the initial values. The results are shown in 
Table 5.6. The objective function value decreased from 2.380 to 1.126. It is found 
that except for flood events 2,5,6, and 8, the fit with optimal parameters is much 
better than that with initial parameters. The flood events 2, 5, 6 and 8 indicate 
very large differences that may be due to very high abstractions or data errors. 

ii) Sankhini Basin: This basin (Fig. 5.8) has a catchment area of 7.76 Km 2 . 
The data for this basin are available only for the period 1964 and 1965. This basin 
has one raingauge station. The main channel is divided into 4 channel segments - 
CH20, CH21, CH24, and CH27. There are four tributaries to the main channel and 
these form 4 more channel segments, viz., CH22, CH23, CH25, and CH26. Thus there 
are 8 channel segments and correspondingly 16 overland flow segments, OFOl to 
OF16. 

There are 5 storm flood events in the basin. The results of the analysis are 
shown in Table 5.7. The results of the simulation run assuming the initial 
parameters for the watershed show that both the simulated peak and volume of 
runoff are much higher than the observed ones. The results of the optimisation run 
are shown in columns 6 and 7 of the Table. The objective function has decreased 
significantly from 4.850 to 0.632 and there is a much better agreement between the 
peak values and the effective rainfall values for the observed and simulated floods. 
For floods 1 and 3, the simulated volume of runoff are higher while for the other 3 
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TABLE 5.7 RESULTS OF ANALYSIS FOR SANKHINI BASIN 


FLOOD OBSERVED VALUES NO OPTIMISATION OPTIMISATION WITH OPT. WITH PSP=10 
NUMBER INITIAL VALUES 


ALL FLOODS ALL FLOODS 



PEAK 

EFF.RAIN 

PEAK 

EFF.RAIN 

PEAK 

EFF.RAIN 

PEAK 

EFF.RAIN 

(1) 

(2! 

(3) 

Hi 

(5) 

(6) 

(7) 

(8) 

(9! 

1 

338 

1.666 

734 

4.526 

376 

2.187 

475 

2.789 

2 

263 

1.869 

444 

3.359 

155 

1.142 

134 

1.007 

3 

279 

2.529 

676 

7.855 

450 

3.776 

481 

3.990 

4 

168 

1.149 

288 

2.570 

124 

1.052 

170 

1.443 

5 

132 

0.642 

336 

2.251 

64 

0.437 

54 

0.373 


PSP 

KSAT 

RGF 

BMSN 

OBJECTIVE INITIAL 

FUNCTION VALUES FINAL 


3.052 

10.000 

0.366 

0.174 

17.404 

6.290 

12.591 

8.020 

4.850 

4.120 

0.632 

1.200 






SIMULATED 

V A L U 

E S 







OPT. WITH PSP=10 

OPT. WITH PSP=10 

OPT. WITH PSP=10 

OPT. WITH PSP=10 

OPT. WITH PSP=10 



RGF=18 

RGF= 1 8 

RGF=18 

RGF=18; 

;BMSN=10 

RGP=18 ;BMSN=10 



ALL FLOODS 

FLOODS 2,4 

FLOODS 1,3 

ALL FLOODS 

FLOODS 2,4 



PEAK 

EFF.RAIN 

PEAK 

EFF.RAIN 

PEAK 

EFF.RAIN 

PEAK 

EFF.RAIN 




PSP 

(10) 

(ID 

(12) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 


KSAT 












RGF 

485 

2.848 

641 

3.827 

268 

1.576 

402 

2.346 

651 

3.894 


BMSN 

136 

1.023 

230 

1.680 

88 

0.675 

114 

0.867 

249 

1.821 

OBJECTIVE 

INITIAL 

483 

4.028 

551 

5.321 

387 

2,977 

449 

3.602 

559 

5.487 

FUNCTION VALUES 

FINAL 

173 

1.464 

206 

1.759 

111 

0.947 

146 

1.236 

210 

1.799 



51 

0,351 

63 

0.431 

31 

0.221 

43 

0.297 

65 

0.44 




10.000 


10.000 


10.000 


10.000 


10.000 




0.171 


0.139 


0.277 


0.208 


0.135 




18.000 


18.000 


18.000 


18.000 


18.000 




8.740 


10.395 


8.045 


10.000 


10.000 




4.120 


0.662 


1.999 


2.660 


1.629 




1.290 


0.192 


0.296 


0.837 


0.201 


TPEAK = TIME TO PEAK IN HOURS PEAK = PEAK DISCHARGE IN CUSECS 
EFF.RAIN = EFFECTIVE RAINFALL IN INCHES 
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floods, they are lower. It may be seen that higher effective rainfall is associated 
with higher floods peaks (flood events 1 and 3) and lower effective rainfall is 
associated with lower flood peaks. Apparently, the infiltration conditions for these 
set of floods are different. 

iii) Mujnai Basin: Mujnai river basin (Fig.5.9) has a catchment area of 38.40 
Km 2 . The data for this basin were available for the years 1964 and 1965 and this 
basin has a raingauge station. A main channel runs through the basin and it has one 
major tributary. The main channel is divided into 4 channel segments and the 
tributary forms the fifth channel segment. Thus this basin has 5 channel segments 
CH20, CH21, CH22, CH23, and CH24 and correspondingly 10 overland flow segments, 
OFOl to OFIO. 

This basin has 10 storm flood events during the period of record. The 
results of the analysis are shown in Table 5.8. The simulation run with the initial 
values of the parameters (columns 4 and 5) indicates that all but one of the 
simulated peaks and all the volumes of runoff are higher than the observed ones. 
The optimisation run with the initial parameters was carried out; the objective 
function decreases significantly from 12.081 to 1.563 and the results appear to be 
much better. Flood events 1„3, and 9 show a higher simulated values; 4, 6, 7 and 8 
have lower simulated values while flood events 2, 5 and 10 have comparable values. 
Thus as in the case of Sankhini river basin, this basin also seem to have at least 
two distinct set of flood events, one having a higher infiltration condition and the 
other lower infiltration condition. 

iv) Khaga Basin : This basin ( Fig.5.10 ) has a catchment area of 



SL .HO OBSERVED VALDES 


TABLE 5.8 RESULTS OP ANALYSIS FOR KUJKAI BASIN 


SIMULATED VALUES 

NO OPTIMISATION OPTIMISATION WITH OPT. WITH PSP=10 OPT. WITH PSP=10 OPT. WITH PSP=10 



PEAK 

EPF. RAIN 

PEAK 

EFF.RAIN 

INITIAL VALUES 

ALL FLOODS 

PEAK EFF.RAIN 

ALL FLOODS 

PEAK EFF.RAIN 

RGF=18 

ALL FLOODS 

PEAK EFF.RAIN 

RGP=18;BMSN=10 
FLOODS 5,6,7,8,10 
PEAK EFF.RAIN 

(1) 

(2) 

(3) 

(4) 

(5) 

(6) 

(T) 

(8) 

0! 

(10) 

(ID 

(12) 

(13) 

1 

1359 

1.313 

4704 

3.722 

3283 

2.091 

3478 

2.278 

4637 

3.627 

3731 

2.531 

2 

851 

0.932 

3145 

3.105 

967 

0.898 

1265 

1.098 

3041 

2.989 

1463 

1.263 

3 

1851 

1.778 

6977 

5.227 

5981 

3.997 

5970 

3.983 

6897 

5.113 

6167 

4.196 

4 

1372 

0.938 

1612 

1.049 

641 

0.519 

634 

0.514 

1474 

0.981 

726 

0.587 

5 

1048 

0.775 

1818 

1.366 

990 

0.759 

980 

0.752 

1767 

1.332 

1142 

0.867 

6 

579 

0.490 

878 

1.178 

408 

0.566 

404 

0.560 

846 

1.141 

474 

0.670 

7 

3028 

2.181 

2990 

3.364 

2220 

2.393 

2209 

2.380 

2962 

3.326 

2421 

2.622 

8 

341 

0.537 

1901 

2.052 

314 

0.310 

302 

0.327 

1766 

1.859 

361 

0.382 

9 

387 

0.658 

1925 

2.172 

696 

0.62 

683 

0.609 

1795 

1.998 

815 

0.730 

10 

1062 

1.078 

2135 

2.600 

1149 

1.088 

1138 

1.075 

2044 

2.447 

1282 

1.261 



PSP 

9.909 

10.000 

10.000 

10.000 


KSAT 

0.198 

0.199 

0.103 

0.173 


RGF 

19.416 

14.234 

18.000 

18.000 


BMSN 

1.025 

4.690 

5.167 

10.000 

OBJECTIVE 

INITIAL 12.081 

12.081 

9.364 

9.364 

1.450 

FUNCTION VALUES 

FINAL 

1.563 

1.620 

8.493 

0.505 


PEAK = PEAK DISCHARGE IN CUSECS EPF. RAIN EFFECTIVE RAINFALL IN INCHES 




FIG-5-10 SEGMENTATION OF KHAGA BASIN 
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66.00 Km 2 and has two raingauge stations. The data for this basin are available 
for the years 1981 and 1984 and from this, 5 flood events were identified for 
analysis. The main channel was divided into 5 channel segments and there are four 
tributaries which form another 4 channel segments. Thus there are 9 channel 
segments, CH20 to CH28 and 18 corresponding overland flow segments, OFOl to OF18. 
The results of the analysis are shown in Table 5.9. 

The simulation run with the initial parameters showed higher values for 
flood peaks and volumes compared to the observed ones for flood events 2, 3, 4, and 
5 while in the case of flood event 1, the simulated values were very much lower. 
Optimisation with the initial values of the parameters showed that for flood events 
1 and 3, the simulated values were very much lower than the observed ones and for 
flood events 4 and 5, they were very much higher while for the flood event 2, the 
simulated and observed values were comparable (Table 5.9). The objective function 
value has decreased from 5.255 to 3.265 and yet is large. Further analysis with 
several alternate parameters indicated that in case flood event 1 were to be 
modelled well, there were very large differences in the case of all other floods. 
This indicates data errors in flood event 1. Simulation with a number of other 
parameter combinations indicated higher peaks. A comparison of times to peak 
indicated the need for introducing additional attenuation and time delay through 
storage elements in the river and this is considered in Subsec. 5.4.3. 

v) Dharsi basin : This basin (Fig. 5.11) has a catchment area of 91.40 Km 2 . The data 
for this basin are available for the years 1968 to 1970 and 12 storm flood events 
were identified. This basin has two raingauge stations. The main channel is divided 
into 5 channel segments, CH20, CH23, CH24, CH24, and CH27. There are 3 major 
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TABLB 5.9 RESULTS OF ANALYSIS FOR KHAGA BASIN 

FLOOD OBSERVED VALUES NO OPTIKISATION 

NUMBER 


1 = 0 Is 2H0URS 



PEAK 

EFF.RAIN 

TPBAK 

PEAK 

EFF.RAIN 

TPEAK 

PEAK 

BFF.RAIN 

TPEAK 

(1) 

(21 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 

1 

13365 

5.626 

7 

6863 

3.443 

5 

9080 

3.443 

7 

2 

5371 

2.106 

12 

13605 

4.692 

2 

12859 

4.692 

11 

3 

1466 

0.615 

9 

1960 

0.957 

6 

2754 

0.95? 

9 

4 

1999 

1.369 

13 

6798 

4.731 

10 

7501 

4.731 

12 

5 

1319 

0.840 

8 

6854 

4.262 

6 

7575 

4.262 

8 


SIMULATED VALUES 



OPTIMISATION WITH 

OPT. WITH PSP=10 

OPT. WITH PSP=10 

OPT. WITH PSP=10 

OPT. WITH PSP=10 

OPT. SITH PSP=10 


INITIAL VALUES 



RGF=18 

RGF=18;BMSN=10 

RGF=18 

RGF=18;BMSN=10 


ALL FLOODS 

ALL FLOODS 

ALL FLOODS 

ALL FLOODS 

FLOODS 2,3,4, 5 

FLOODS 2,3 


PEAK 

BFF.RAIN 

PEAK 

EFF.RAIN 

PEAK 

EFF.RAIN 

PEAK 

BFF.RAIN 






(11) 

(12) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

(20) 

(21) 

(22) 


4906 

2.198 

4627 

2.089 

4724 

2.134 

5561 

2.341 

3615 

1.632 

10253 

4.655 


5340 

2.464 

5119 

2.376 

5208 

2.415 

7564 

3.475 

3896 

1.841 

12241 

5.643 


794 

0.421 

773 

0.422 

676 

0.393 

162 

0.128 

442 

0.277 

1085 

0.563 


3914 

2.894 

4132 

2.937 

4166 

2.959 

4643 

3.488 

3559 

2.349 

8921 

7.142 


4652 

3.048 

4924 

3.375 

4962 

3.422 

5190 

3.724 

4588 

3.000 

6241 

5.411 


PSP 

3.147 


10.000 


10.000 


10.000 


10.000 


10.000 


IS AT 

0.369 


0.172 


0.165 


0.123 


0.240 


0.024 


RGF 

18.286 


6.250 


18.000 


18.000 


18.000 


18.000 


BMSN 

6.310 


6.874 


6.000 


10.000 


6.087 


10.000 

JBCTIVE 

INITIAL 

5.255 


4.465 


4.466 


5.860 


3.962 


2.339 

NCTION VALUES 

FINAL 

3.265 


3.655 


3.726 


5.821 


6.087 


0.979 


TPEAK = TIME TO PEAK IN HOURS PEAK = PEAK DISCHARGE IN CUSECS 
EFF.RAIN = EFFECTIVE RAINFALL IN INCHES 
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FIG- 5-11 SEGMENTATION OF DHARSI BASIN 
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tributaries joining the main channel at various locations and they are designated as 
channel segments CH21, CH22, and CH25. Thus there are 8 channel segments and 
correspondingly 16 overland flow segments, OFOl to OF16. 

The results of analysis are shown in Table 5.10. The initial simulation run 
showed that all the simulated values were very much higher than the observed 
ones. A comparison of times to peak, as in the case of Khaga basin, indicated the 
need for introducing additional attenuation and time delay through storage elements 
in the river. The optimisation run (columns 11 and 12) reduced the objective 
function from 31.897 to 6.707. Also for flood events 1, 2, 6, and 11, the simulated 
values were lower while for the others, it was higher. Only for flood event 7, they 
are comparable. Thus for this basin also, two distinct groups of flood events with 
different infiltration conditions were identified. 

vi) Mansari basin: This is the largest among the basins considered in th* 
study. The catchment area of the basin (Fig. 5.12) is 213.00 Km 2 and it has ' 
raingauge stations. The data for this basin were available for the years 1970 an 
1971 and 5 storm flood events were identified. The main channel is divided into 
channel segments, CH21, CH22, CH24, CH26, CH28, CH30. 5 major tributaries drai 
into the main channel at various locations and each of these form a channel segmer 
and they are designated CH21, CH23, CH25, CH27, and CH29. Thus there are 1 
channel segments and 22 corresponding overland flow segments in this basin. 

The results of the analysis are shown in Table 5.11. Simulation with initi 
values of parameters showed that for all the floods, the simulated values of tl 
peak were very much higher compared to the observed peaks. The volumes 
runoff were also high except for flood number 1. The simulation run wi 



TABLE 5.10 RESULTS OF ANALYSIS FOR DBARSI BASIN 


FLOOD OBSERVED VALUES 

NUMBER NO OPTIMISATION 



PEAK 

EFF.RAIN 

TPEAK 

PEAK 

EFF.RAIN 

TPEAK 

PEAK 

EFF.RAIN 

TPEAK 





(i; 

1 (2) 

(3) 

(4) 

(5) 

(6) 

(7) 

(8) 

(9) 

(10) 





i 

743 

0.499 

13 

5678 

2.814 

6 

3421 

2.814 

12 





2 

3037 

3.548 

35 

7456 

7.414 

7 

3766 

7.414 

34 





3 

2522 

3.082 

39 

10445 

8.047 

15 

4270 

8.047 

27 





4 

1379 

1.291 

15 

5879 

7.601 

6 

3613 

7,601 

11 





5 

1245 

0.781 

13 

10287 

7.047 

4 

3841 

7.047 

21 





6 

1713 

1.178 

17 

5645 

2.653 

13 

3398 

2.653 

18 





7 

880 

0.717 

10 

5056 

1.941 

5 

2930 

1.941 

9 





8 

1773 

1.059 

14 

9138 

7.508 

30 

3810 

7.508 

33 





9 

1658 

1.001 

11 

8375 

4.378 

4 

3692 

4.378 

11 





10 

582 

0.434 

9 

6492 

4.449 

5 

3621 

4.449 

10 





11 

688 

0.532 

16 

5505 

2.357 

8 

3201 

2.357 

14 





12 

702 

0.484 

16 

7617 

3.717 

9 

3640 

3.717 

14 






OBJECTIVE 


INITIAL 


31.897 



31.897 







FUNCTION VALUES 

FINAL 

















SIMULATED 

VALUES 







OPTIMISATION KITH 

OPT. WITH PSP=10 

OPT. WITH PSP=10 

OPT. WITH PSP=10 

OPT. WITH PSP=10 

OPT. NITH PSP=lfi 



INITIAL VALUES 



RGF=1 8 

RGF=18 FLOODS 

RGF=18 

RGF=18;BNSN=10 



ALL FLOODS 

ALL FLOODS 

ALL FLOODS 

2,3,4 »5,6,8,9 

FLOODS 2,4,6 

FLOODS 2,4,6 



PEAK 

EFF.RAIN 

PEAK 

EFF.RAIN 

PEAK 

EFF.RAIN 

PEAK 

EFF.RAIN 

PEAK 

EFF.RAIN 

PEAK 

EFF.RAIN 



(ID 

(12) 

(13) 

(14) 

(15) 

(16) 

(17) 

(18) 

(19) 

(20) 

(21) 

(22) 



552 

0.429 

536 

0.411 

549 

0.420 

374 

0.305 

1531 

1.095 

2142 

1.342 



2254 

2.430 

2232 

2.440 

2257 

2.485 

1600 

1.730 

2848 

3.579 

3183 

4.047 



3576 

4.665 

3612 

4.754 

3629 

4.812 

3429 

3.892 

3710 

5.428 

3723 

5.479 



1338 

3.556 

1428 

3.561 

1461 

3.619 

1037 

2.704 

1839 

4.225 

1982 

4.36 



3725 

4.143 

3745 

4.248 

3749 

4.300 

3629 

3.473 

3789 

4.862 

3796 

4.964 



843 

0.617 

653 

0.490 

704 

0.525 

328 

0.274 

1322 

0.921 

1448 

0.999 



833 

0.557 

737 

0.502 

781 

0.527 

251 

0.210 

1322 

0.825 

1442 

0.885 



2181 

1.732 

2668 

2.180 

2272 

1.778 

1447 

1.147 

2840 

2.442 

2344 

1.899 



2542 

1.612 

2619 

1.619 

2664 

1,722 

2219 

1.316 

2883 

2.019 

2943 

2.074 



957 

0.873 

1058 

0.922 

1083 

0,943 

752 

0.648 

1377 

1.197 

1443 

1.251 



465 

0.376 

496 

0.398 

508 

0,407 

350 

0.297 

688 

0.526 

744 

0.563 



2217 

1.312 

2111 

1.245 

2172 

1.282 

1783 

1.054 

2365 

1.443 

2889 

1.948 



PSP 

13.890 

10.000 

10,000 

10.000 

10.000 

10.000 


KSAT 

0.205 

0.245 

0.240 

0.342 

0.188 

0.180 


RGF 

14.680 

16.010 

18.000 

18.000 

18.000 

18.000 


BMSN 

7.730 

6.350 

7.316 

8.970 

6.299 

10.000 

OBJECTIVE 

INITIAL 

31.897 

26.910 

25,029 

12.730 

3.934 

2.726 

FUNCTION VALUES 

FINAL 

6.707 

7.500 

7,260 

5.550 

1.465 

0.049 


TPEAK = TIME TO PEAK IN HOURS PEAK = PEAK DISCHARGE IN CUSECS 
EFF.RAIN * EFFECTIVE RAINFALL IN INCHES 




TABLE 5.11 RESULTS OF ANALYSIS FOR MANSARI BASIN 


FLOOD OBSERVED VALDES NO OPTIMISATION 

NUMBER 



PEAK 

EFF.RAIN 

TPEAK 

PEAK 

EFF.RAIN 

TPEAK 

(1) 

(21 

(3) 

U) 

(5! 

(6) 

(7! 

1 

368 

0.025 

13 

47 

0.025 

5 

2 

913 

0.582 

25 

8601 

0.922 

16 

3 

523 

0.242 

22 

5281 

0.376 

15 

4 

1414 

1.102 

25 

6503 

1.408 

17 

5 

1583 

0.706 

21 

9851 

1.127 

15 




SIMULATED 

V A L U 

E S 




OPTIMISATION WITH 

OPT. WITH PSP=10 

OPT. WITH PSP=10 

OPT. WITH PSP=10 

OPT. WITH PSP=10 

INITIAL VALUES 


RGF=18 

RGF=18 

RGF=18;BMSN=10 

ALL FLOODS 

ALL FLOODS 

ALL FLOODS 

FLOODS 2, 3,4, 5 

FLOODS 2, 3, 4, 5 

PEAK 

EFF.RAIN 

PEAK EFF.RAIN 

PEAK 

EFF.RAIN 

PEAK 

EFF.RAIN 

PEAK 

EFF.RAIN 

(8) 

(9) 

(10) (11) 

(12) 

(13) 

(14) 

(15) 

(16) 

(17) 

31 

0.026 

29 0.024 

27 

0.023 

13 

0.015 

0.07 

0.0002 

2704 

0.922 

2397 0.826 

2185 

0.766 

1717 

0.616 

1611 

0.622 

794 

0.376 

739 0.409 

799 

0.44 

478 

0.256 

433 

0.243 

3303 

1.408 

3055 1.323 

2314 

1.066 

2131 

0.938 

2117 

0.976 

3082 

1.127 

2788 1.026 

2016 

0.78 

1928 

0.735 

1950 

0.786 



PSP 

6.000 

10.000 

10.000 

10.000 

10.000 


KSAT 

0.127 

0.095 

0.092 

0.128 

0.133 


RGF 

9.023 

6.970 

18.000 

18.000 

18.000 


BMSN 

4.511 

4.250 

1.122 

3.522 

10.000 

OBJECTIVE 

INITIAL 

4.517 

4.749 

6.060 

0.396 

0.230 

FUNCTION VALUES 

FINAL 

2.779 

2.765 

2.797 

0.034 

0.031 


TPEAI = TIME TO PEAK IN HOURS PEAK = PEAK DISCHARGE IN CUSBCS 
EFF.RAIN = EFFECTIVE RAINFALL IN INCHES 
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optimisation using the initial set of parameter values also showed that the simulated 
peaks were very high. The objective function decreased slightly from 4.517 to 
2.779, showing a large deviation between observed and simulated hydrographs during 
the flood events. Simulation with a number of combination of parameters also 
indicated high peaks. A comparison of times to peak indicated the need for 
additional storage elements in the river and this is considered in Sec .5 .4. 3. 

5.4.3 Simulation with storage 

The three basins, Khaga, Dharsi and Mansari, are larger in catchment area in 
comparison to the other basins. The necessity of introducing a storage component in 
the watershed model was indicated in the earlier phase of the study (Sec.5.4.2). 
Linear storage elements with storage discharge relationship S - KQ with storage 
coefficient K can be used in the model at appropriate locations. For simplicity, a 
single lumped linear reservoir at the gauging site was adopted in this study. Trial 
values of storage coefficients around the indicated values were tested without 
optimisation and with optimisation to decide on the storage coefficient K and 
appropriate parameters of the distributed parameter model. For example, results of 
analysis for Dharsi basin are shown in Table 5.12. This indicates a storage 
coefficient of K - 6 hours and the corresponding set of parameters. Similar studies 
for Khaga and Mansari showed a delay time of 2 hours and 9 hours respectively. 

5.4.4 Different infiltration conditions 


For four basins, viz., Sankhini, Mujnai, Khaga and Dharsi, there are two 



TABLE 5.12 RESULTS OF ANALYSIS WITH STORAGE FOR DHARSI BASIN 




OBSBRVED 




SIMULATED 

(WITHOUT OPTIMISATION) 


FLOOD 

PEAK EFF.RAIN 

TPEAK 

PEAK EFF.RAIN 

TPEAK 

PEAK 

TPEAK 

PEAK 

TPEAK 

NUMBER 




K = 1 

HOUR 


K = 

3 HOURS 

K = 

6 HOURS 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

I 

743 

0.499 

13 

4409 

2.814 

7 

3821 

13 

3421 

15 

2 

3037 

3.548 

35 

5187 

7.414 

8 

4406 

9 

3766 

34 

3 

2522 

3.082 

39 

5998 

8.047 

16 

4469 

33 

4270 

36 

4 

1379 

1.291 

15 

4568 

7.601 

9 

3864 

15 

3613 

16 

5 

1245 

0.781 

13 

5857 

7.047 

8 

4975 

8 

3841 

11 

6 

1713 

1.178 

17 

4506 

2.653 

13 

3979 

16 

3398 

17 

7 

880 

0.717 

10 

404 1 

1.941 

6 

3597 

10 

2930 

11 

8 

1773 

1.059 

14 

5463 

7.508 

31 

4253 

32 

3810 

33 

9 

1658 

1.001 

11 

5280 

4.378 

5 

3837 

20 

3692 

21 

10 

582 

0.434 

9 

4665 

4.449 

6 

3896 

9 

3621 

10 

11 

688 

0.532 

16 

4453 

2.357 

8 

3852 

11 

3201 

15 

12 

702 

0.484 

14 

4862 

3.717 

10 

3839 

22 

3640 

22 




OBSERVED 




SIMULATED 

(WITH OPTIMISATION) 



FLOOD 

PEAK EFF.RAIN TPEAK 

PEAK EFF.RAIN 

TPEAK 

PEAK 

TPEAK 

PEAK 

TPEAK 

NUMBER 




K = 

1 HOUR 


K = 

3 HOURS 

K = 

6 HOURS 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

1 

743 

0.499 

13 

1044 

0.424 

9 

764 

9 

552 

12 

2 

3037 

3.548 

35 

3803 

2.43 

32 

3092 

33 

2254 

34 

3 

2522 

3.082 

39 

4677 

4.605 

15 

3938 

20 

3576 

27 

4 

1379 

1.291 

15 

2242 

3.556 

9 

1808 

11 

1338 

11 

5 

1245 

0.781 

13 

5210 

4.143 

8 

4127 

9 

3725 

21 

6 

1713 

1.178 

17 

1589 

0.617 

15 

1177 

17 

843 

18 

7 

880 

0.717 

10 

1768 

0.557 

7 

1212 

8 

833 

9 

8 

1773 

1.059 

14 

3855 

1.732 

32 

3069 

33 

2181 

33 

9 

1658 

1.001 

11 

3839 

1.612 

10 

3355 

11 

2542 

11 

10 

582 

0.434 

9 

2001 

0.873 

8 

1371 

9 

957 

10 

11 

688 

0.532 

16 

841 

0.376 

12 

638 

13 

465 

14 

12 

702 

0.484 

14 

3841 

1.312 

13 

3044 

13 

2217 

14 



PSP 


13.893 


PEAK = PEAK DISCHARGE IN CUSECS 




KSAT 


0.205 


EFF.RAIN = 

EFFECTIVE RAINFALL IN INCHES 



RGF 


14.68 


TPEAK = TIME TO PEAK IN HOURS 





BMSN 


7.727 








OBJECTIVE 

INITIAL 

31.93 








FUNCTION VALUE FINAL 


6.711 
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distinct groups of floods due to differring infiltration conditions. For example, for 
Sankhini basin (Table 5.7), flood events 2 and 4 pertain to a low infiltration 
condition while flood events 1 and 3 correspond to a high infiltration condition. 
Since the present study is concerned with estimation of design flood, the parameter 
estimates corresponding to those of lower infiltration conditions are considered to 
be appropriate for the purpose of the study. Hence optimising the two parameters 
KSAT and BMSN with PSP - 10.00 and RGF - 18.00 for floods 2 and 4 of Sankhini 
basin, the values of these parameters are obtained as KSAT = 0.139; and BMSN — 
10.395. For the floods under higher infiltration conditions, these values are 0.277 
and 8.045 respectively. These are shown in columns 12 to 15. For the lower 
infiltration conditions, the simulated values of peak and the effective rainfall are 
found to be slightly low for flood event 2 while it is high for the flood event 4. 
The values for the other floods are farther away, being very high in the case of 
floods 1 and 3 and very low in the case of flood 5. The objective function value 
reduces form 0.662 to 0.192 in this case. 

Similarly for Mujnai basin, (Table5.8), flood events 1,2,3 and 9 (Fig.5.12a) 

pertain to a low infiltration condition while flood events 6,7, and 8 correspond to a 
high infiltration condition. Flood events 5 and 10 have comparable values and flood 
event 4 corresponds to a very wet condition. For further analysis flood events 
5, 6, 7,8 and 10 (Fig.5.12b) are only considered since for design flood estimation, 
parameters corresponding to wet infiltration conditions are considered to be 
appropriate. 

Similarly for Dharsi basin (Table 5.10), flood events 2, 3, 4, 5, 6, 8, and 9 are 
the major floods. Keeping the values of PSP = 10.00 and RGF = 18.00, the 
simulation run for these floods yield values of KSAT = 0.342 and BMSN — 8.97. 
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DR3M - MUJNAI FLOOD #1 DR3M - MUJNAI FLOOD # 2 




FIG. 5.12 a. SIMULATED HYDROGRAPHS OF MUJNAI BASIN (Dry initial condition) 

This value of KSAT indicates that there is a dominance of some floods and higher 
infiltration conditions. Flood events 2, 4, and 6 correspond to lower infiltration 
conditions during actual storm as is evident from the lower values of simulated 
peaks and effective rainfall and correspondingly the higher values of historical 
flood peaks and effective rainfall. Simulating these three floods with the constant 
values of PSP and RGF, yield values of KSAT = 0.188 and BMSN - 6.299. The 
objective function value for this condition has an initial value of 3.934 and a final 
value of 1.465. 

For Khaga basin (Table 5.9), when simulation is carried out for all the floods. 
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DR3M — MUJNAI FLOOD # 5 


DR3M - MUJNAI FLOOD # 6 




DR3M - MUJNAI FLOOD # 7 


DR3M - MUJNAI FLOOD # 8 




DR3M - MUJNAI FLOOD # 10 


DR3M - MUJNAI FLOOD # 4 




FIG.5.12 b. SIMULATED HYDROGRAPHS FOR MUJNAI BASIN 
Floods 5,6,7,8 and 10 — > Wet initial condition; Flood 4 --> Very wet initial condition 
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Table 5.13 Model parameters after initial simulation 


Co nstant par ame ters 

: EVC 

- 0.70 

; RR - 

0.90; DRN 

- 0.50 

Bas i n 

Area 

Km 2 

PSP 

inches 

KSAT 
inches 
per hr 

RGF 

BMSN 

inches 

Objective 

Function 

Initial Final 

K 

Hours 

Batjhora 

7.33 

7.586 

0.112 

16.46 

1.89 

2.380 

1.126 


Sankhini 

7.76 

3.052 

0.366 

17.40 

12.59 

4.850 

0.632 


Mujnai 

38.40 

9.909 

0.198 

19.42 

1.025 

12.081 

1363 


Talma 

42.10 

8.179 

0.167 

16.63 

7.731 

8.699 

5.762 


Khaga 

66.00 

3.147 

0369 

18.29 

6.310 

5.255 

3.265 

2.0 

Dharsi 

92.40 

13.890 

0.205 

14.68 

7.730 

31397 

6.707 

6.0 

Mansari 

213.00 

6.00 

0.127 

9.02 

4311 

4316 

2.779 

9.0 


the other three parameters are optimised, the objective function value at the end 
of the optimisation run is 5.999. This differs only slightly from the earlier value 
of 5.762. It is also found that by changing the value of PSP to 10.00, there is only 
marginal change in the values of flood peaks and the effective rainfalls and in a 
few cases there is an improvement. Similarly, comparison of columns 8 and 9 with 
columns 10 and 11 shows that when PSP - 10, changing RGF to 18.00 produces no 
significant change in the final objective function value (6.011) or in the peak and 
the effective rainfall values. Similar studies for all the basins indicate that the two 
parameters PSP and RGF can be taken as constant within the observed range for all 
the basins with values of 10.00 and 18.00 respectively. (Table 5.13). The results of 
the simulation run with optimisation with the values of PSP — 10.00 and RGF — 
18.00 for the model are shown in Tables 5.5 to 5.13 and are summarised in the Table 


5.14. 
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Table 5.14 Model Parameters Values with PSP - 10; and RGF — 18.00 


Constant Parameters: PSP - 

10.00; RGF - 18.00; EVC- 

0.70; RR 

-0.90; DRN-0.S0 

Basin 

KSAT 

BMSN 

Objective Function 

Initial Final 

Batjhora 

0.096 

1.180 

0.30 

0.22 

Sankhini 

0.171 

8.740 

4.12 

1.29 

Mujnai 

0.103 

5.170 

9.36 

8.49 

Talma 

0.163 

8.198 

7.53 

6.01 

Khaga 

0.165 

6.000 

4.47 

3.73 

Dharsi 

0.240 

7.320 

25.03 

7.30 

Mansari 

0.092 

1.122 

6.06 

2.80 


The corresponding values of peak and the effective rainfall to this run are also 
shown in Tables 5.5 through 5.11. 

It is seen from Table 5.14 that the BMSN values range from 1.12 to 8.74 
inches. A number of studies with varying values of BMSN were done and it was 
observed that the model is not very sensitive to this parameter. A value of 10 
inches for BMSN corresponding to an upper limit (Table 5.13) is adopted. Thus 
KSAT is found to be the only sensitive parameter of the model for simulating 
floods in the region. Simulation is carried out to optimise this parameter (KSAT) for 
each of the basins and the results are presented in the last two columns of Tables 
5.5 through 5.11. It may be seen that on the basis of peak and effective rainfall of 
each of the floods and the value of the objective function the results are 
comparable to those of the respective previous columns. Accordingly constant 
parameter values of PSP -10.00 inches; RGF- 18.00; and BMSN- 10.00 inches are 
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indicated as regional parameters values and they are adopted in the study. KSAT, 
however, was one of the parameters to which the results of simulation were 
sensitive. The values of KSAT for various basins, assuming the other parameters as 
constant are shown in Table 5.15 and Fig. 5. 14 


Table 5.15 Variation of KSAT in the region among the basins 


Constant Parameters: PSP -=10.00; RGF“ 18-00; 

RR — 0.90;DRN — 0.50 

BMSN = 10.00; EVC-0.70; 

Basin 

KSAT 

Objective Function 


inches/hour 

Initial 

Final 

Batjhora 

0.102 

1.226 

1315 

Sankhini 

0.135 

1.629 

0301 

Mujnai 

0.173 

1.450 

0305 

Talma 

0.182 

0343 

0303 

Khaga 

0.123 * 

5.860 

5.820 

Dharsi 

0.180 

2.726 

0.449 

Mansari 

0.133 * 

0.230 

0.031 


* Values of the parameter inconsistent. 

It is observed from the Table and Fig. 5.14 that the values of KSAT for the 
smaller basins, viz., Batjhora and Sankhini are very low, apparently due to the 
presence of impervious areas which are not considered in the present study so far. 
The values of KSAT for the other three basins considered are 0.173 inches/hour 
for Mujnai, 0.182 inches/hour for Talma and 0.180 inches/hour for Dharsi. The 
average value of KSAT for these three basins is 0.178 inches/hour (4.53 mm/hour) 
and this compares well with 5.00 mm/hour value of infiltration index used in Sec. 4.6 
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and agrees well with the silty loam or sandy clay loam soil present in the region. It 
was decided to use this value of KSAT for the smaller basins with directly 
contributing impervious percentage of catchment area as a parameter. Simulation was 
carried out with a value of KSAT ** 0.178 inches/hour and a few trial runs showed 
that 20% of the area as impervious for Batjhora and 10% for Sankhini, the results 
of the simulation are comparable to the earlier runs. It was then decided to use this 
value of KSAT and simulate floods in all the basins and the results are shown in 
Table 5.16. These results seem to be realistic and hence are accepted in the study. 

For the region under study, the following parameters are found to be regional 
parameters; 


PSP 

10.000 inches 

KSAT 

0.178 inches/hour 

RGF 

18.000 

BMSN 

10.000 inches 

EVC 

0.700 

RR 

0.900 

DRN 

0.50 


5.6 Application of DR3M to Kurti Basin 

Kurti basin was studied in Subsec. 4.5.1 and rejected for further study 
because of suspected data errors. Results of simulation of floods for this basin 
using DR3M with the model parameters for the region with historical data for 
rainfall are shown in Table 5.17. They show wide variation between the observed 
and Simula '^d flood runoffs and peaks. The floods and peaks are much larger in the 
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TABLE 5.16 RBSULTS OF ANALYSIS FOR ALL THE BASINS USING THE REGIONAL PARAMETERS 
BATJHORA SANKHINI HOJNAI 

OBSERVED SIMULATED OBSERVED SIMULATED OBSERVED SIMULATED 


PEAK 

EFF.RAIN 

PEAK 

EFF.RAIN 

PEAK 

EFF.RAIN 

PEAK 

EFF.RAIN 

PEAK 

EFF.RAIN 

PEAK 

EFF.RAIN 

1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

160 

0.356 

123 

0.251 

338 

1.666 

492 

2.924 

1359 

1.313 

3659 

2,465 

399 

1.322 

1708 

1.503 

263 

1.869 

171 

1.318 

851 

0.932 

1418 

1.220 

310 

0.848 

303 

0.727 

279 

2.529 

506 

4.542 

1851 

1.778 

6102 

4.135 

233 

0.590 

345 

0.696 

168 

1.149 

189 

1.647 

1372 

0.938 

699 

0.565 

310 

1.399 

974 

1.064 

132 

0.642 

83 

0.564 

1048 

0.775 

1094 

0.833 

813 

3.779 

1057 

2.051 





579 

0.490 

453 

0.638 

813 

2.470 

650 

1.688 





3028 

2.181 

2361 

2.556 

163 

0.567 

548 

1.144 





341 

0.53? 

346 

0.36? 

194 

0.385 

486 

0,705 





387 

0.658 

779 

0,694 









1062 

1.078 

1246 

1.207 

ERVIOUS 


20 




10 







:A PERCENTAGE 

TALMA KHAGA DHARSI 

OBSERVED SIMOLATED OBSERVED SIMULATED OBSERVED SIMULATED 


PEAK 

EFF.RAIN 

PEAK 

EFF.RAIN 

PEAK 

EFF.RAIN 

PEAK 

EFF.RAIN 

PEAK 

EFF.RAIN 

PEAK 

EFF.RAIN 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 

23 

24 

179 

0.14? 

121 

0.263 

13365 

5.626 

5450 

2.394 

743 

0.499 

2116 

1.329 

396 

0.732 

82 

0.199 

5371 

2.106 

5645 

2.635 

3037 

3.548 

3168 

4.022 

739 

1.548 

202 

0.441 

1466 

0.615 

750 

0.342 

2522 

3.082 

3122 

5.554 

4704 

8.401 

5004 

8.374 

1999 

1.369 

4783 

3.318 

1379 

1.291 

1979 

4.354 

251 

0.418 

56 

0.145 

1319 

0.840 

5483 

3.587 

1245 

0.781 

3796 

4.959 

501 

0.710 

721 

1.369 





1713 

1.178 

1441 

0.955 

1213 

2.359 

1147 

2.499 





880 

0.717 

1436 

0.882 

625 

0.357 

436 

0.735 





1773 

1.059 

2333 

1.894 

3093 

5.997 

3936 

5.680 





1658 

1.001 

2942 

2.071 

3523 

9.741 

2846 

7.006 





582 

0,434 

1441 

1.248 

2703 

3.126 

2875 

4.820 





688 

0.532 

742 

0.562 

598 

0.650 

517 

0.96? 





702 

0.484 

2886 

1.940 


MANSARI 

OBSERVED SIMULATED 

PEAK EFF.RAIN PEAK EFF.RAIN 

25 26 27 28 

368 0.025 51 0.038 

913 0.582 1932 0.719 

523 0.242 586 0.395 

1414 1.102 2606 1.196 

1583 0.706 1492 0.899 


PARAMETERS: 

PSP 10.000 INCHES 

KSAT 0.178 INCHES/HOUR 

RGF 18.000 

BMSN 10.000 INCHES 

EVC 0.700 

RR 0.900 

DRN 0.500 


PEAK = PEAK DISCHARGE IN CUSECS 
EFF.RAIN = EFFECTIVE RAINFALLIN INCHES 
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case of flood events 1 to 3 and less in other cases. Perhaps as in the case of 
Sankhini, Dharsi and Khaga basins, the results indicate different infiltration 
conditions, viz., a wet infiltration condition for floods 1 to 3 resulting in high 
floods and dry infiltration conditions for floods 4 and 5. The results confirm the 
earlier suspicion of data errors. 

Table 5.17 Results of simulation of floods for Kurti basin 


Flood 


Observed 


Simulated 


No. 

Actual Peak 

Ft 3 /sec M 3 /sec 

Effective 

Rainfall 

Inches 

Peak 

Ft 3 /sec M 3 /sec 

Effective 

Rainfall 

Inches 

1 

1067 

30.22 

2.440 

2601 

73.77 

7.491 

2 

1067 

30222 

1.879 

1912 

54.15 

42270 

3 

2244 

63.55 

5.946 

4099 

116.1 

14.650 

4 

1476 

41.80 

2.608 

745 

21.09 

1.400 

5 

1393 

39.45 

2.855 

661 

18.72 

1.700 

5.6.1 

Estimating design 

flood peak for 

Kurt i 

basin 



To estimate the design flood for the Kurti basin, the antecedent condi tons 
for design storm are to be simulated in terms of specified wet initial soil moisture 
conditions. In case storm precipitation data for the given basin were not available, 
it is possible to use precipitation data from an adjacent raingauge station or basin in 
the simulation. The design storm could have been introduced when the simulated 
discharge corresponds to design baseflow (Fig .4 .11). However, since precipitation data 
are available for the basin, design storm was introduced immediately following the 
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day of the highest simulated flood, viz., flood event 3 and the simulation was 
carried out using the design storm of duration 2 hours which is the time of 
concentration for the basin, and a 12 hour design storm. Since the simulated flood 
peak for flood 3 is of the order of 116.1 m 3 /sec, much larger than the design flood 
peak, to consider other initial wet conditions the design storm is introduced 
separately also immediately following the day of the second highest simulated flood 
(Qpeak " 73.77 m 3 /sec), viz., flood event 1 and the results are also summarised in 
Table 5.18. 


Table 5.18 Results of design flood estimation 


Design storm Peak discharge value Baseflow Peak discharge 

duration Ft 3 /sec M 3 /sec Ft 3 /sec M 3 /sec with baseflow 

hours of 25.44 M 3 /sec 

(1) (2) (3) (4) (5) <6) 


Case (i) Design storm follows highest simulated storm flood event 

2 3037 86.01 1259 35.63 75.82 

12 3157 89.41 1259 35.63 79.22 

Case Design storm follows second highest simulated storm flood event 

2 2722 77.09 835 23.65 78.88 

12 3030 85.81 835 23.65 87.60 


Since the baseflows are different for the two cases studied, as an approximation, a 
correction for difference in baseflow is applied to the simulated peak discharges and 
the corrected peak values are shown in column 6. The results indicate that the 
design peak discharges are comparable to those of Sec. 4.7. Furthermore a large 
initial baseflow perhaps may lead to a smaller peak discharge. A 12 hour design 
storm hyetograph seems to indicate a realistic design peak in both cases. 
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Accordingly the following steps for estimation of design flood hydrograph for 
an ungauged basin is suggested. 

i) From available physiographic and hydrometeorological data regional 
parameters for the model are estimated including their variations, if 
any, in space or as a function of area. 

ii) Physiographic data for the basin are collected and adapted in the model. 

iii) Using precipitation data for the basin, if avialable, or else precipitation 
data for a nearby station or basin, the watershed is simulated to identify 
the larger peaks. Data for a couple of years is considered to be adequate. 

iv) Using IMD data, a 12 hour design storm for the basin is estimated and is 
introduced on the day when the simulated discharge is approximately 
equal to the design baseflow for the basin. Alternatively the design storm 
may be introduced the day next to that of the high peaks. 

v) By simulation the design flood is obtained. In case the baseflow is 
different from the design baseflow, the design flood hydrograph may be 
may be appropriately corrected. 

5.7 Conclusions 

The study indicates that DR3M with constant parameters for the region 
performs satisfactorily for the small basins in terms of prediction of peaks and 
effective rainfall. From consideration of parsimony of parameters, a simple constant 
parameter model is considered to be realistic representation of small basins in the 
region. 

The conclusions have been reached in a heuristic manner. The UH approach 
could be used in case adequate rainfall-runoff data are available or if a regional 
relationship on the basis of UH has already been obtained. But since the 
relationships are essentially derived based on basin area or similar parameters 
empirically, they may be subject to large errors and are not recommended. 

Topography, channel crosssections and soil data may be easily availble for 
field organisations and can be used with physically based regionally derived 
parameter estimates iri a distributed parameter model like DR3M to derive design 
flood. Accordingly a distributed parameter model is recommended for flood 
estimation in small catchments except when at least some rainfall - runoff data are 
available for the basin. 



CHAPTER 6 


SUMMARY, CONCLUSIONS AND SUGGESTIONS 

FOR FUTURE WORK 

6. 1 Summary: 

Estimates of flood are required for the design and economic appraisal of a 
variety of engineering works such as culverts, spillways etc. Limited or non- 
availability of data is a major constraint in flood estimation for watersheds. 

Several approaches are available for estimation of floods. When adequate and 
appropriate data are available for the basin, frequency analysis may be applied or 
else unit hydrograph (UH) and other rainfall -runoff procedures may be applied to 
model the basin system. When this is combined with design rainfall, the design flood 
may be derived. Data available for small watersheds are generally inadequate for 
these approaches. Regional flood frequency analysis is generally not valid for small 
watersheds because of data limitations. Then other empirical approaches for 
regionalisation are to be used to model either the variation of flood peaks or 
rainfall-runoff relationships. 

From available data regression relationships are generally derived to relate 
peak discharge to hydrometeorological features of the storm and physiographic 
features of the watershed. Even when adequate data are available for establishing 
such relationships, they are limited to the regions for which such relationships 
have been derived. Since generally adequate data are available for estimating design 
storms, it seems worthwhile to derive relationships between storm rainfall and 
streamflow, perhaps in terms of unit hydrographs. The parameters of the UH may 
then be related to the characteristics of the storm and the physiographic features 
of the watershed to yield a regional relationship for design purpose between storms 
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and floods. 

A more realistic approach is to model the catchment runoff process according 
to physical principles. Data concerning topography, soil and vegetation are available 
from remote sensing. The physically based models can now be implemented even in a 
PC environment. There has also been a greater availability of technical know-how 
on the applicability of such models. Hence, it seems worthwhile to consider the use 
of physically based distributed parameter models as an alternative to regionalisation 
for flood estimation in small catchments. 

Information regarding physiography, soils, vegetation and the drainage 
pattern of the catchment is used in the physically based models. The catchment is 
subdivided into a network of channel and overland flow segments. Using the physical 
characteristics of each component and appropriate hydraulic equations, the 
hydrologic processes are simulated using the distributed parameter model. Necessary 
values of the coefficients and constants are specified from physical characteristics 
of the watershed or obtained from calibration of the model with observed data. 

For purposes of hydrometeorological studies, India has been divided into 26 
subzones by the Central Water Commission (CWC), each of which is assumed to be 
hydrologically homogeneous. Research Designs and Standards Organisation (RDSO) of 
the Ministry of Railways has collected, and in collaboration with CWC and India 
Meteorology Department (IMD) has analysed a large amount of hydrometeorologicl 
data for a large number of catchments draining into various railway culverts and 
bridges in different subzones. This study is essentially a regional case study for 
flood estimation in small catchments in a part of the north Brahmaputra basin, 
designated as subzone 2(a), to compare particularly in a very limited data 
environment the two approaches for rainfall-runoff relationships, viz.. 
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regionalisation using UH and a physically based model. For a number of small 
watersheds in this region, available hydrometeorological data were very limited and 
inadequate and so data available were yet to be analysed and interpreted. 

The objectives of the present study are: 

i) To develop a regional relationship for UH parameters from available limited 
data for a number of basins to evaluate the validity and applicability of UH 
based regionalisation procedure for flood estimation. 

ii) To apply a physically based model to small watersheds in the region to 
evaluate the validity of such a model to model the regional variation, if any, of 
the parameters of the model and also to judge their validity for flood 
estimation, and 

iii) To compare these two approaches for regional flood estimation for small 
catchments. 

The study is limited to small watersheds in the western part of the north 
Brahmaputra basin, for which data were available. Furthermore, the UH is 
represented by the Nash conceptual model and the physically based model considered 
in the study is the USGS Distributed Routing Rainfall Runoff Model (DR3M). 

Catchment plan, crosssection at the bridge site and hourly rainfall and stage 
data at railway bridges for about 10 catchments in the region were available for the 
study. The catchment areas ranged from about 7.00 Km 2 to 215.00 Km 2 . The data 
were available for a period ranging from two to three years. After preliminary 
analysis for data errors, only data of seven basins were found to be suitable for 
analysis. 

In the UH approach for regionalisation, it is necessary to represent the UH in 
terms of a limited number of parameters. A conceptual model with two parameters N 



and K proposed by Nash was adopted in the study. The various steps involved in 
the regionalisation of the parameters of this model are: 

i) Estimate the UH parameters for a storm in a basin, 

ii) Repeat the step for a number of storms in a basin to explain the variation, if 
any, of the parameters in terms of hydrometeorological characteristics of the 
storm, and 

iii) Repeat the above two steps for a number of basins in the region so that the 
residual variation of parameters can be correlated to the physiographic 
characteristics of the basin. 

Various parameters like the coefficient of runoff, the infiltration index and 
baseflow are found to vary from storm to storm in each basin rather erratically. 
Attempts made to correlate them with storm and runoff parameters were 
unsuccessful. The variation of the infiltration index among the storms was very 
large. Correlation with a number of storm flood parameters also showed wide 
variation. However, a general study of the individual storms indicated that a value 
of 5.00 mm/hour can be considered to be the general minimum value for the 
infiltration index for all the basins, particularly during intense storms generally met 
with in the design. This value corresponds to the silty loam soil or sandy clay loam 
soil of the region. 

The UH parameters N and K of the Nash model vary from storm to storm and 
from basin to basin. A smaller NK generally leads to a quicker flood and a high 
peak. A quasilinear variation of N and K as a function of peak discharge (Q pe , k ) is 
assumed and for each of the basin, the regression equations relating N, K and NK 
with Q p))(lk is derived. Regional regression equations relating N and K corresponding 
to the maximum observed peak in each basin are as follows: 
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N - 1.792 + 0.0088 A, and K - 2.766 + 0.01 6A 
where A is in Km 2 , K is in hours and N is dimensionless. 

The baseflow varies from storm to storm and generally increases with the 
basin area. Baseflow/Km 2 for all critical storms are plotted as a function of basin 
area and an enveloping curve is drawn. Since the higher values of baseflow lead to 
larger floods, the enveloping curve for baseflow/Km 2 is recommended for use with 
the design flood. 

A design flood estimated using the regional relationships for a small 
watershed in the region is derived and is compared with other approaches generally 
used in India. The results are found to be satisfactory. 

Many distributed parameter models are available for simulating rural and 
hilly catchment runoff. The DR3M has been tested for some rural and hilly basins 
in India with satisfactory results and hence this model was adopted for the study. A 
drainage basin is represented in this model as a set of channel and overland flow 
segments in such a way that the essential basin geometry and physiography are 
taken into account for runoff computation. The model has four components - a soil 
moisture or water balance component, the infiltration or rainfall excess component, a 
routing component and an optimisation component. The model has seven physically 
based parameters - four soil moisture accounting parameters, viz., (i) DRN, a 
constant drainage rate for redistribution of soil moisture in inches/day, (ii) EVC, a 
pan coefficient for converting measured pan evaporation to potential 
evapotranspiration; (iii) RR, the proportion of daily rainfall that infiltrates into the 
soil for period of simulation excluding unit rainfall days and (iv) BMSN, the soil- 
moisture storage at field capacity in inches, and three infiltration parameters viz., 
(i) PSP, the suction at wetted front for soil moisture at field capacity in inches of 



pressure, (ii) KSAT, the effective saturated value of the hydraulic conductivity in 
inches/hour and (iii) RGF, the ratio of suction at wetted front for soil moisture at 
wilting point to that at field capacity. These parameters are fitted to the model for 
a basin by Rosenbrock’s algorithm and the iterative optimisation among basins is 
done heuristically. The results of the study indicated that a constant parameter 
model for the small catchments in the region is realistic. The model parameters 
agree with the watershed characteristics. 

It is realised from the study that the directly contributing impervious area 
percentage to be given as input data from the characteristics of the watershed is 
very important in estimating the flood peak. Such data were not available for the 
region. They were parameterised and estimated by simulation and compared with 
general topography and basin area and are considered to be satisfactory. However, 
with the use of remote sensed data (not available for the study but generally 
available for government and other field organisations), it will be possible to use 
this information for realistic modelling of small watersheds. The design flood 
estimated for a small watershed using this model compared favourably with other 
estimates. 

6.2 Conclusions: 

1. For estimation of design flood in a small basin, since design rainfall can be 
estimated separately by standard procedures, it seems preferable to use a 
regional rainfall-runoff relationship with design storm rather than using a purely 
empirical approach. 

2. The UH parameters are not generally constant. Hence it is necessary to adopt 
a quasilinear approach and criteria for estimation and regionalisation of desien 



parameters. Regional UH approaches involve a large amount of empiricism. When 
sufficiently large amount of data are available they may be used for 

regionalisation. 

3. For the design flood for the region under study, a constant infiltration 
parameter during the design storm, a baseflow related to the basin area and a 
regional quasilinear relationship for the UH parameters were found to be 

satisf actory. 

4. DR3M is a relatively simple distributed parameter model which can be used 

for simulation of runoff in small watersheds. 

5. Impervious area is a very important parameter in the DR3M which has to be 

estimated with care. 

6. Even in a very limited data environment, DR3M has been found to work 

satisfactorily. 

7. DR3M may be used to realistically estimate the design flood, particularly 
when data are very limited, since field data concerning topography, soil and 
vegetation may be easily collected by remote sensing and other means or by 

modelling and regionalisation. 

8. Because of easy availability of, and competence with computers and softwares 
in field organisations, it seems preferable to adopt a distributed parameter model 
with relevant field data rather than empirical, statistical relationships of 

questionable validity. 

9. Distributed parameter models, may, eventually replace empirical regression 
models in hydrology not only for large watersheds but also for small watersheds. 



6.3 Suggestions for Further Study 


Based on the results of the study, the following suggestions are made for 

further study: 

1. The study used limited data of questionable reliability. It seems worthwhile 

to use reliable data of sufficiently long length. 

2. Correlation and regionahsation of model parameters with physiographic 
characteristics of the basin need further study. 

3. Comparison of results with other distributed parameter models may identify 
an appropriate model for use in the case of small watersheds. 

4. It is also necessary to test such models in a variety of regions to identify 
variation of parameters in space and as a function of area. The variation, 

if any, may be used for defining homogeneity of a region and identification 
of hydrologically homogeneous regions. 
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